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ABSTRACT 


An  analytical  and  experimental  program  was  conducted  to  develop  acoustic  fatigue  design 
criteria  for  aircraft  structures  subjected  to  intense  noise  in  a  high  temperature  environment. 
Equations  for  the  dynamic  response  of  a  buckled  panel  were  formulated  for  simply  supported 
boundary  conditions  using  large  deflection  plate  theory.  Random  amplitude  acoustic 
fatigue  testing  of  representative  aircraft  structure  was  accomplished  at  temperatures  up  to 
600°F  to  provide  data  for  correlation  with  the  analytical  results.  Empirical  design  criteria 
are  presented  in  the  form  of  design  equations  and  nomographs  for  predicting  the  combined 
thermal  and  dynamic  response  of  aircraft  structures. 
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Kinetic  energy;  temperature  increase  (note  context) 

Temperature  of  box  structure  skins,  see  Appendix  I 
Kinetic  energy  of  plate  element 


T 

r 


Temperature  of  box  structure  rib,  tee  Appendix  I 
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Kinetic  energy  of  substructure 
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Critical  buckling  temperature  of  center  bay  of  nine-bay  panel  -  F  obove 
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Ciritical  buckling  temperature  of  simple  panel  -  °F  above  ambient 
Time 

Thickness  of  box  structure  ribs  and  skin,  see  Appendix  I 
Strain  energy 

Strain  energy  of  plate  element 
Strain  energy  of  substructure  elements 
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panel  -  inch 
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Coordinate  direction 

Coordinate  direction  from  stiffening  member  centroid 
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Centroid  of  data  set,  see  Section  IV 
Coefficient  of  thermal  expansion  -  (in/in)/°F 

Cross-section  warping  constant  of  substructure  about  shear  center  -  in^ 

Cross-section  warping  constants  of  substructure  about  the  skin  attachment 
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Standard  deviation  from  the  mean 
Mean  stress  -  psi  (or  ksi) 

Thermal  expansion  stress  -  psi  (or  ksi) 

Normal  stress  in  the  x  and  y-directions  -  psi  (or  ksi) 

Thermal  mean  stress  in  x  and  y-directions  -  psi  (or  ksi) 

Thermal  buckling  stress  in  x  and  y-directions  -  psi  (or  ksi) 

Dynamic  stress  in  x  and  y-directions  -  psi  (or  ksi)  rms 
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Abbreviations 

dB  Decibel  (Re:  0.0002  microbar) 

F  Fahrenheit 

ksi  1000  psi 

psi  lb/in^ 

rms  Root  mean  square 

RT 


SPL 


Room  temperature 
Sound  pressure  level 


I  -  INTRODUCTION 


The  development  of  acoustic  fatigue  design  criteria  has  historically  followed  on 
empirical  approach,  since  many  of  the  factors  that  affect  high  cyclic  fatigue  are  not  pre¬ 
dictable.  Ecrly  efforts'  were  directed  toward  establishment  of  o  broad  base  of  general 
design  information  in  the  form  of  empirical  design  nomographs,  applicable  to  several 
structural  configurations.  These  design  charts  were  based  on  discrete  frequency  fatigue 
data  that  had  been  converted  to  an  equivalent  random  amplitude  fatigue  curve  through 
the  use  of  Miles^  single  degree -of- freedom  theory  and  the  Miner-Palmgren®  cumulative 
damage  concept.  The  aircraft  designer  found  these  criteria  lo  be  of  considerable  value  and 
he  could,  with  a  certain  amount  of  individual  judgment  and  experience,  effectively 
control  acoustic  fatigue  problems  for  conventional  subsonic  aircroft.  However,  increasing 
aircraft  performance  and  size,  coupled  with  increasingly  stringent  requirements  for  structural 
efficiency,  created  a  need  to  refine  these  conservative  criteria  to  eliminate  unnecessary 
weight  from  structural  designs. 
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Subsequent  programs  were  accomplished  to  refine  and  extend  the  range  of  appl icat ion 
of  the  then-existing  design  criteria.  These  later  programs  included  acoustic  fotigue 
tests  of  structural  panels  using  random  amplitude  excitation.  The  program  of  Reference  4 
included  fatigue  tests  of  flat,  sti ffened-sk in  and  honeycomb  panels.  A  subsequent  program® 
continued  the  refinement  of  stiffened-skin  design  criteria  by  developing  empirical 
design  methods  for  the  structural  support  members. 

The  results  of  these  previous  programs  have  considerably  reduced  the  uncertainties 
involved  in  predicting  dynamic  response  and  life  characteristics  of  conventional 
structures  subjected  to  both  propulsion  system  and  aerodynamic  noise  at  ambient  tempera¬ 
tures.  However,  when  unusual  structural  configurations  or  environmental  conditions 
are  encountered,  the  applicability  of  existing  design  methodology  decreases,  and  the 
judgment  of  the  design  engineer  must  be  relied  on  more  heavily.  Structural  applications 
are  commonly  encountered  today  in  the  near  field  of  an  operating  turbojet  or  turbofan 
engine  where  severe  noise,  high  temperatures,  static  loading,  and  vibratory  buffet 
occur  simultaneously  or  in  conjunction  with  each  other. 

At  least  one  attempt1"*  has  been  made  to  define  the  effects  of  low  frequency  dynamic 
loads  and  elevated  temperatures  on  the  acoustic  fatigue  resistance  of  flat  and  curved 
simple  panels,  and  curved  honeycomb  structures.  One  of  the  significant  conclusions 
from  this  program  was  that  combined  environments  should  be  investigated  in  limited  com¬ 
binations;  the  combined  environment  problem  should  be  approached  gradually  and  system¬ 
atically.  The  primary  cause  of  the  fatigue  failures  experienced  during  this  program  was 
dynamic  response  due  to  acoustic  excitation,  acting  in  combination  with  thermal  mean 
stresses  caused  by  the  structural  heating. 

The  program  described  in  this  report  was  undertaken  to  extend  tk  basic  design  technology 
for  stiffened  structures  at  ambient  temperatures  to  include  she  effects  of  simultaneous 
application  of  thermal  and  acoustic  environments.  An  analytical  development  is 
described  in  Section  II  for  the  dynamic  response  of  heated  structures  before  and  after 
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rhermal  buckling.  The  primary  purpose  of  Hie  analytical  effort  was  to  identify  the  para¬ 
meters  which  describe  the  structural  response;  then  the  data  requirements  for  the  ex¬ 
perimental  program  (Section  111)  were  defined  in  detail.  Measured  data  were  correlated 
with  the  analytical  results  in  Section  IV;  these  empirical  relations  were  then  used  to 
establish  the  design  methods  and  nomographs  presented  in  Section  V.  A  preliminary 
analytical  development  for  the  dynamic  response  of  heated  box  and  curved  structures 
is  presented  in  Appendices  I  and  II. 

The  results  of  this  investigation  are  also  summarized  in  o  design  handbook  for  ease  of 
application  by  the  design  engineer.  That  report,  AFFDL-TR-73-1 55,  Part  II,  contains  the 
empirical  equations,  design  nomographs,  and  computer  programs  described  herein. 


II  -  ANALYTICAL 


The  emphosis  of  the  onalytical  effort  was  placed  upon  developing  techniques  for  under¬ 
standing  the  vibration  of  heated  structures  through  maximum  utilization  of  available  ana¬ 
lytical  results  reported  in  the  literature.  Since  the  beginning  point  of  this  development 
was  tiie  analysis  previously  developed  by  Rudder,  the  primary  task  was  the  modification  of 
the  room  temperature  stiffened  panel  analysis  to  include  the  effect  of  a  spatially  uniform 
temperature  increase  upon  the  dynamic  characteristics  of  the  structure. 

It  was  assumed  that  the  panel  vibration  modes  and  the  panel  buckled  modes  are  such  that 
the  stiffeners  remain  straight  along  the  axial  direction  and  only  rotate  about  the  attach¬ 
ment  line  of  the  stiffener  and  the  cover  sheet.  The  analytical  model  considers  on  array  of 
simply  supported  panels  constrcined  at  the  boundaries  so  that  the  slopes  between  adjacent 
panels  are  compatible  (zero  shear  condition).  This  permits  the  use  of  sinusoidal  mode 
shapes  below  and  above  the  panel  buckling  temperature.  The  validity  of  this  assumption 
was  checked  by  using  the  results  presented  by  Timoshenko^  for  a  stiffened  ponel  uniformly 
compressed  by  inplone  edge  forces. 

The  Rayleigh-Ritz  method  is  used  to  derive  the  equations  of  motion  for  the  structure.  An 
expression  is  developed  for  predicting  the  temperature  increase  required  to  cause  the  ponel 
tc  buckle,  this  temperature  increase  is  defined  as  the  critical  temperoture .  The  tempera¬ 
ture  increase  is  defined  relative  to  a  room  (or  ambient)  temperature  at  which  no  thermal 
mean  stresses  are  present.  Below  the  critical  temperature  the  panel  dynamic  response  is 
calculated  uti  I  izing  linear  small  deflect'on  plate  theory.  Above  the  critical  temperature, 
large  deflection  (von  Karman)  plate  theory  is  used.  The  large  deflection  plate  analysis  is 
linearized  to  obtain  first  order  effects  for  estimating  the  ponel  dynamic  response  above  the 
critical  temperature.  All  of  the  assumptions  are  discussed  os  they  ore  introduced,  The 
effect  of  a  temperature  increase  upon  material  properties  is  also  considered  (i  ,e . ,  the 
temperature  dependence  of  Young's  Modulus  and  the  coefficient  of  expansion). 

Based  on  the  simple  panel  analysis,  expressions  for  the  frequency  and  buckling  amplitude 
of  a  flat ,  stiffened ,  n  ine  -bay  panel  structure  are  developed  .  These  relations  inc  lude  the 
effect  of  substructure  stiffness  and  mass;  however,  the  substructure  is  assumed  to  be 
temperature-independent . 


A.  Simple  Ponel  Response 

The  geometry,  nomenclature,  and  sign  convention  for  the  inplone  forces  is  indicated  in 
Figure  I ,  which  shows  the  simple  panel  configuration  .  Assuming  simple  supports  at  the  edges 
of  the  panel,  the  dynamic  and  buckled  modes  of  the  panel  can  be  described  by  sinusoidal 
functions  in  the  inplone  coordinates.  The  analysis  assumes  a  spatially  uniform  temperoture 
increase,  T,  over  the  surface  of  the  panel.  Below  the  critical  temperature,  the  thin  panel 
is  assumed  to  be  initially  flat,  and  small  deflection  plate  theory  is  used  to  develop 
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the  equations  of  motion  .  Above  the  critical  temperature,  the  thin  panel  is  assumed  to 
have  a  budded  shape,  and  the  equations  of  motion  are  developed  by  application  of  large 
deflection  (von  Kantian)  plate  theory. 


FIGURE  1  .  SIMPLE  PANEL  GEOMETRY  AND  NOTATION 


I .  Pre-Buckled  Response 

In  the  following  discussion,  the  panel  skin  temperature,  T,  refers  to  the  temperature 
Increase  reloti  ve  to  a  predefined  zero  mean  stress  condi  tion  (i  .e  . ,  ambient  temperature) . 

Assuming  o  transverse  displacement  function  of  the  form 

w(*,y,t)=^^  sin  (”"<?” )  S'n  (  )  ^mn(t}  (1) 

m  n 

end  using  small  displacement  plate  theory,  the  strain  energy  of  the  simple  panel  has  the 
form 

a  b  a  b 

U  =  J"  J  (c^w)  dydx  i  ^  h  j"  J'  •»  tt  -m~  dydx  (2) 

0  0  0  0 

Ai.ile  the  kinetic  energy  has  the  form 

a  b 

T  =  ^  Yh  J'  J'  w2  dydx  (3) 

0  0 
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Here,  and  throughout  this  development,  the  spatial  and  time  dependence  functions  have 
been  dropped  from  the  basic  parameter  after  the  initial  introduction  to  simplify  the 
expressions . 


For  small  deflection  plate  theory,  the  mean  stresses  are  defined  in  terms  of  the  tempera¬ 
ture  increase,  T,  os 


a  - 

x 


a 

y 


F.aT 

=  _n^y 


(4) 


Substituting  the  assumed  displacement  function  and  the  obove  mean  stresses  into  the  strain 
and  kinetic  energy  expressions,  the  equation  of  motion  for  the  (m,n)  mode  is  determined 
by  Logronge's  equation,  and  is 


"4DF2 

•  ■  mn  _  .  . 

q  -*  - tt-tt  (1  -C  r  q  -  0 

mn  ,  l.  mn  'mn 

vha  b 


(0  *c  r<:  1) 

mn 


(5) 


where  C  -  F  , ,  ,/F  :  F 

mn  -II.  mn  mn 
temperature,  T  ,  is 


(b/a)m  ■*  (a/b)n  ;  r  -  T/f  ;  and  the  critical  buckling 


T 

c 


2,2_ 

n  h  F 


II 


1 2cLab(  1  -*•  v) 


(5a) 


Assuming  simple  harmonic  motion,  the  expression  for  the  natural  frequency  is 


fW  °  2ob 


m 


(I  -C  r) 

mn 


mn 


(c  r  1) 

mn 


(6) 


It  is  now  advontogeous  to  discuss  certain  features  of  these  results.  First,  the  critical 
temperature,  Tc ,  is  defined  in  terms  of  the  panel  geometry,  with  the  only  material 
properties  being  the  coefficient  of  linear  expansion,  a,  and  Poisson's  ratio,  v.  Noting 
that  the  frequency  expression,  given  by  Equation  (6),  is  a  function  of  the  mode  number 
(m,n),  it  is  obvious  that  the  frequency  squared  for  a  given  mode  decreases  linearly  with  an 
increase  in  temperature.  The  magnitude  of  the  temperature  differential  required  to  yield  a 
zero  natural  frequency  for  a  given  mode  is  broadly  defined  as  the  critical  buckling 
temperature  for  that  mode.  However,  from  Equation  (5)  it  is  seen  that  the  lowest  critical 
temperature  occurs  for  the  fundamental  mode  (m  =  n  -  1),  ond  that  the  critical  temperatures 
for  the  higher  modes  are  related  to  the  fundamental  mode  critical  temperature  by  the  panel 
aspect  ratio  and  mode  numbers. 


Since  the  frequency  expression  given  by  Equation  (6)  is  limited  to  positive  or  zero  values, 
this  analysis  will  be  used  to  estimate  the  panel  fundamental  frequency  response  for  a  range 
of  temperature  such  that  T  ^  Tc .  In  other  words,  this  analysis  will  be  limited  to  tempera  ¬ 
ture  increases  below  the  critical  temperature.  Also,  oil  references  to  the  critical 
temperature  will  imply  the  temperature  increase  required  to  cause  the  fundamental 
frequency  to  equal  zero. 
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Filially,  it  should  ho  not  til  iliut  Hu.-  equations  ot  motion  tui  each  mud-;  aie  uncoupieii 
below  tlie  critical  tempcuilui  ••  .  Above  the  ci  i  r  ica !  iriufinrnti'K;,  M'll.'nin,  1  ^  shows  ti.al  tin; 
buckled  panel  dynamic  i  espouse  is  coupled  for  the  dynamic  modes  different  frum  the 
buckled  mode,  but  the  dynamic  mode  •_ an  expanding  to  die  truckled  mode  is  oil  oupled 
from  all  other  dynamic  modes.  This  result  will  be  applied  to  estimate  the  simple  panel 
response  above  the  critical  temperature. 


2.  Post-Buckled  Response 

The  analysis  presented  here  is  a  simpl  ication  of  the  analytical  approach  devcluped  by 
Shulmanl  ‘  for  calculating  the  dynamic  response  of  a  thin,  rectangular  ,r..nel  subjected 
to  a  spatially  uniform  temperature  increase.  Since  the  dynumic  response  of  the  buckled 
panel  is  required,  large  deflection  plate  theory  is  used  .  Basically,  the  von  Korman 
theory  ^ ^  is  used  w  i  fh  an  assumed  disp  lac  ement  and  Airy  stress  func  t  ion  .  ^  1  The  expres¬ 
sion  for  the  strain  energy  of  the  buckled  panel  is  formulated  from  the  assumed  displace¬ 
ment  made  shape  and  the  str  ess  function  .  The  nonlinear  form  of  the  stress  function  is 
linearized  la  obtain  a  first  order  approximation  for  the  equations  of  motion. 

Since  this  analysis  is  aimed  at  providing  approximate  results,  certain  assumptions  wilt  be 
introduced  based  on  the  pie-buckled  analysis  (below  the  critical  temperature)  .  First,  it 
is  assumed  that  the  px.net  buckles  in  the  luiidaniental  mode  (m  n  -  1)  corresponding  to  the 
lowest  buckling  temperature.  Then,  it  is  assumed  that  the  only  significant  dynamic  moae 
is  the  fundamental  mode  of  lire  panel  . 

The  amplitude  of  the  buckled  ■aural  mg  panel  is  defined  by 

W,n(t)  -  Wq-  q(t) 


where  W  is  the  static  buckled  panel  amplitude,  ond  q(t)  is  the  dynamic  panel  amplitude 
(both  corresponding  to  the  fundamental  mode  of  the  panel).  The  assumed  transverse  dis¬ 
placement  has  the  form 

w(x,;,t)  =-•  sin  |  ~-)  Sin  (-^)w)|(.)  (7) 

For  trie  panel,  the  stress- ;ti uin  relationships  or  ; 


x  (I  -  v)  _  2)  --x  / J 


xy  2(1  -t  \)  ’  xy 


E  J  _  E 

/  ft  ri~ 


(I  ■  u  ) 


and  the  strain  •  displacement  relationships  are 
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(9) 


1  2 

€  -  u ,  •+  —  w .  ;  Y  -  u ,  *  v.  ■+  w.  w , 

x  *  2  x  xyyxxy 

.  1  2 
e  -  v,  +  x-  w, 
y  y  2  y 

For  flie  assumed  mode  given  by  Equation  (7),  a  suifoble  choice  for  the  stress  function  Isas 
the  form  11,12 


F(x,y,t)  =  ^  (°x y2  +  °yx2)  +  f  j(x,y,  t) 


(10) 


vhere 


F^x.y.t)  »Eh  \  N  [(£)  m2s  (g)"2]  cos(^~ ) 


iCOSlb  )fJ° 


m-0  n-0 


and  (m,n)  take  on  integer  volues.  The  significance  of  the  mean  stresses  o  and  will 
appear  '!  ortly.  The  stress  function  (10)  and  the  displacement  function  (7)  ore  related  by 


v4F  -  -7^Nt  +  Eh(w^  -  w,  w,  )  (11) 

T  xy  xx  yy 


/here 


h/2 

NT  =  a£  f  T(*'y-z)dz 
-h/2 

2 

For  a  spatially  uniform  temperature  increase,  V  s  0. 

The  relationships  between  the  stress  function,  F,  and  the  stresses  given  in  Equation  (8)  are 

lv7  =  F ,  ;  ha  =  F,  ;  hT  -  ~F,  (12) 

x  yy  y  xx  xy  xy 

Substituting  the  assumed  stress  function.  Equation  (10),  into  Equations  (12),  it  is  evident 
thut  the  shear  stress  vanishes  along  the  edges  of  the  panel . 

The  coefficients  fmn  are  obtained  by  substituting  the  sumed  mode,  Equation  (7),  into 
Equation  (11),  with  the  result  that 

12  12 
f  =  -  W*  ■  f  =  -  Vv 
02  2  11  '  20  2  11 

while  all  other  coefficients  vanish  identically. 

Substituting  these  coefficients  into  Equation  (10),  the  stress  function  becomes 
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(13) 


F(*,y,t) 


h  2  - 

o  (3  y  ^  °  x 

Ik  y 


^  -*  F,(x,y,t) 


Substitution  •  ihe  stress  function  into  Equotions  (12)  results  in  t fie  static  stress  relations 

1 


?  -  r  F ,  =  C7  -  cos  (  -r~Z- )  W ^ 

X  h  yy  x  „  2  o  ' 


r  F , 

h  xx 


8a- 

2, 


11 


.“E  :  2nX  :  2 

T - *cos - )wn 

y  s?  0  11 


(14) 


T  =  0 

xy 


N 3w ,  substituting  the  strcin-displocement  relationships  foi  and  ey,  Equation  (9),  into 
the  stress-strain  relationships  for  -Tx  and  3y/  Equation  (8),  and  solving  for  the  displace¬ 
ments  u,  arid  v,  gives 
x  'v  a 


u , 

X 


2 

'x 


V 


T  V 

-•  aT  -  T  a  ■+ 
y  ex 


1  2 


(15) 


Integrating  these  displacements  and  assuming  that  the  edges  of  the  panel  are  rigidly 
restrained  from  moving  in  the  (x,y)  plane  yields 


a  a  a  a 

f  u'*d*  ~  =  0  -  afiT  +  ~  J  a^dx  -  ^  J  a^dx  -  ^  J  w^dx 


b  b 

dy  -  v(b)  -  v(0)  -  0  -  baT 

L. 

0 


D 

A 


■  F  /  V*  ’  F  /  y*  -  5  /  w' 


,?  dy 

y 


(16) 


Combining  the  expressions  for  the  stresses  from  Equation  (14)  and  the  assumed  displacement 

fur-- Mon  from  Equation  (7)  with  Equations  (16)  and  solving  for  the  thermally  induced  mean 

stn.-tvos,  n  and  3  ,  produces 
x  y  r 
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rr^y 
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(1  'V)  8ab(l-v2) 
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To  summarize  the  results  thus  far,  it  is  seen  that  for  the  assumed  displacement  and  stress 
functions  the  stresses  are  given  by  Equation  (14)  with  the  meon  stresses  defined  by  Equation 
(17) .  The  shear  stresses  vanish  identically  for  the  assumed  form  of  the  stress  function. 
Figure  2  illustrates  the  thermal  stress  variation  with  temperature  for  a  typical  panel . 

Now,  in  terms  of  the  displacement  function,  w,  and  the  stress  function,  F,  the  expression 
for  the  strain  energy  is  taken  in  the  form  presented  by  Shulmanl  1 

U  r  7  I  f  ^Vdydx  -  ^h  f  f  {■2fV2dydx 
0  0  0  0 
a  b 

*  5  If  -xw2x+VVdyd* 

0  0 

Substituting  the  displacement  and  siress  functions,  together  with  the  mean  stresses,  the 
above  stroin-eneigy  expression  becomes 

.2' 


(18) 


U  = 


m40 

8a  b 


where 


Fll(1  ~  r>;  +  8 


2  2 


hi 


W 


1 1 


J 


jwn 


(19) 


R  -  3[(5  -  v‘)F ‘  j  -  2(5  f  v)()  -v)] 


(19a) 


and  r  and  T^  remain  as  previously  defined. 

Now  the  displacement  function  was  assumed  to  have  an  amplitude  of  the  form  W ^  j (t)  - 
WQ+  q(t),  where  WQ  is  the  static  buckled  panel  amplitude  and  q(t)  is  the  dynamic 
(vibratory)  panel  amplitude.  Then,  expansion  of  the  displacement  function  gives 


2  2  2 
W7 .  =  W  +2W  qtq' 
!  o  o^  ^ 


(20) 


W, ,  -  W4  +  4W^q  t  6W2q2  •+  4W  q2  +  q4 
O  O  0  o 


1 1 


Considering  the  stotic  problem,  the  dynamic  amplitude  q(t)  in  Equations  (20)  is  set  to  zero 
and  the  strain  energy  minimized  with  respect  to  the  buckle  amplitude,  W0,  to  obtain 
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*u 

?w 

o 


»  dw  r  _ 

_ £  1 F2  (I 

4a b  111' 


r)+iR 


-  0 


(r  *  1) 


Then 


W 

o 


2h  FnL(r  -  1)^]1/2 


(21) 


(22) 


describes  the  panel  static  buckling  amplitude  as  a  function  of  the  temperature  increase. 

The  effect  of  increasing  temperature  on  buckling  amplitude  is  shown  in  Figure  3  for  various 
aspec  t  ratios . 


Now,  it  is  clear  that  if  the  expanded  forms  of  the  displacement  function.  Equations  (20) , 
are  used  to  obtain  the  modal  stiffness  of  the  panel,  the  stiffness  will  be  nonlinear  in  the 
coordinate  q(t) .  Assuming  that  the  static  buckled  amplitude  is  much  greater  than  the 
dynamic  amplitude  (WQ  »  |q(t)|),  then  expressions  for  W~  and  can  be  approximated 
by 


9  ^  9 

W7,  %  W*  +  q 
11  o 

4  4  2  2 

W, ,  *=  W  -  6W  q 

t  1  o  o 


(23) 


Substituting  these  approximations  in  Equation  (19),  the  expression  for  the  strain  energy 
becomes 


or 


(24a) 


2  2  21 
-r)W‘+  2F,1(r  -  l)q  J 


(24b) 


where  Equation  (22)  hos  been  used  to  obtain  the  simplified  Equation  (24b) . 


The  kinetic  energy  of  the  buckled  panel  has  the  form 

c  b 


T  -  iyh 


f  I 

0  0 


W2  dydx 


(25) 


Substituting  Equation  (7)  into  Equation  (25)  and  performing  the  indicated  operations  yields 
the  resulr 


T  - 


1  L  U-2 

o  Vhabq 

o 


I  1 


(26) 


FIGURE  3.  ANALYTICAL  BUCKLING  AMPLITUDE  VARIATION  WITH  TEMPERATURE 
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The  equation  of  motion  is  obtained  by  using  Lagrange's  equation  and  expressions  (26)  and 
(24)  for  the  kinetic  and  strain  energy,  respectively. 


J Yhob’  *  Jar  2Fn(r  ■ ,,q  =  0 


(27) 


or 


4  7x1 1 F  > 


(r  *  1) 


(28) 


Assuming  harmonic  motion,  the  expression  for  the  response  frequency  is 

1  /2 

w-ss'iitlSrh'-')] 


(29) 


Comparing  this  result  to  the  fundamental  frequency  below  the  critical  temperature, 
Equation  (6),  it  will  be  noted  that  as  the  panel  temperature  is  increased  from  the  room 
temperature  condition  (zero  thermal  stress),  the  panel  fundamental  mode  decreoses  to  zero 
at  the  critical  temperature.  Above  the  critical  temperature  (after  the  panel  buckles),  the 
fundamental  frequency  increases  at  a  rate  1 .414  times  that  of  the  sub-critical  temperature 
frequency  decrease.  Defining  the  room  temperature  fundamental  mode  frequency  as  f0, 
and  using  Equations  (6)  and  (29),  the  fundamental  frequency  can  be  simplified  to 


,  ,riR\^y2 

o  2ab  LvhJ 


f(r)  =  f  [1  -  r] 
o 


1/2 


(r  =  0) 


(Oir^l) 


(30) 


f(r)  -  fQ[2(r  -  I) j 


,1/2 


(r  ^  1) 


These  results  are  graphically  shown  in  Figure  4;  the  frequency  rotio  decreases  to  zero  at 
buckling  and  then  increases  after  buckling. 


3.  Summary  of  Simple  Panel  Results 

For  ease  of  comparison  and  reference,  the  following  results  are  presented  in  terms  of  the 
panel  temperature  increase,  T,  relative  to  the  critical  temperature  increase,  T  ,  for  the 
panel .  c 


Critical  Temperature  Increase,  T  : 

2U2C 

tt  h  F . 


Tc  =  12aob(l  +  vf  :  F 1 1  "  (b/a)  +  (a/b) 
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f(r)/fn  -  FREQUENCY  RATIO 


FIGURE  ANALYTICAL  FREQUENCY  VARIATION  WITH  TEMPERATURE 
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Panel  Temperature  Ratio: 


r 


Strain  Energy  (Linearized): 


Panel 


(0  <  r  «  1) 

u  "  E3T  ljf  li*1  '  r>Wo 1  2Ftl(r "  '*q2| 

(r  *  1) 

Ampl  i  tude: 

W  =0 

(0  •-  r  ^  1) 

o 

V\  -  2hF.  '(r-l)Al1/2 
o  1  1 

(r  -  1) 

R  -31(5  -v2)F^  -2(5*  v)(1  -v). 


Mean  Stresses: 


(0  s  r  ?  1) 


EaT 

2p 

n  E  ; 

i  (  —  I  +  v  1 

-)|w 

vn^r 

8ab(l  -v2) 

r  a  1 

b  ')  i 

EoT 

A  | 

t 

_Q  1  O 

> 

ig»lw 

7y  (1  -v) 

8ab(l  -v2) 

Thermal  Stresses: 


2, 

-  TT  E  ;  L 
7=0 - COS 

y  y  8a2 


-  *  )w2 

a  o 


(r  *  1) 
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Equation  of  Motion  (Linear  Free  Vibration): 

..  -4< 

q  + - ^ 1  '  r  0 

Yho  b 

q  *  - *-L-2(r  -  l)q  =0 

Yha  b 


Response  Frequency: 


TsrLvhj 


1/2 


2ab 


f(r)  =  f  rJ  -  r],/2 
o 


f(r)  =  f  '//(r  -  1)] 
o 


1/2 


(0  r  1) 

(r  •*  1) 


(r  =  0) 


(0  s  r  s  1) 


(r  *  1) 


6  .  Nine-Boy  Panel  Response 


The  simple  panel  analysis  described  in  the  previous  section  will  now  be  used  to  develop  a 
model  of  a  flat,  nine-bay  stiffened  panel  configuration,  illustrated  in  Figure  5.  This 
anclysis  is  based  upon  the  results  developed  by  Rudder, 8  and  follows  the  basic  assumptions 
that  slope  compatibility  relations  are  preserved  across  the  attachment  line  of  the  stiffener 
to  the  skin  panel  and  that  the  stiffener  remain  straight  (i.e.,  no  bending  or  buckling  of 
the  stiffeners  is  allowed). 

The  nine-bay  panel  analysis  will  be  presented  in  terms  of  the  temperature  increase  of  the 
panel,  similar  to  that  for  the  simple  panel  analysis.  Since  the  simple  panel  analysis 
assumed  a  uniform  temperature  distribution  over  the  surface  oreo  and  the  thickness  of  the 
panel,  it  is  likewise  assumed  that  the  temperature  distribution  over  the  surface  and  through 
the  thickness  of  the  sheet  (see  Figure  5/  of  the  nine-buy  panel  is  uniform.  This  assumption 
implies  that  the  edges  end  the  stiffeners  of  the  nine-bay  panel  are  insulated  from  the  cover 
sheet. 


1_  Substructure  Characteristics 

The  results  presented  in  Reference  8  can  be  used  directly,  since  the  substructure  is  assumed 
to  remain  at  room  temperature. 

Assuming  that  the  stiffeners  only  rotate,  and  imposing  slope  compatibility  for  the  cover  sheet 
across  the  stiffeners,  the  expression  for  the  strain  energy  0f  stiffeners  parallel  to  the  x-axis  is 
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UNIFORM  TEMPERATURE,  T 


b)  SUPPORT  STRUCTURE  GEOMETRY 


FIGURE  5.  NINE-BAY  FLAT  STIFFENED  PANEL  CONFIGURATION 
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The  stiffener  cross-section  geometry  and  nomenclature  is  illustrated  in  Figure  6. 


2.  Skin  Characteristics 


Since  the  kinetic  energy  of  the  cover  sheet  is  independent  of  temperature,  this  result  will 
be  presented  first.  Assuming  a  fundamental  mode  response  for  the  nine-bay  sheet  (i.e., 
each  bay  responds  in  its  fundamental  mode)  and  applying  slope  compatibility  relations. 
Table  I,  across  the  stiffener  attachment  lines,  the  kinetic  energy  is 


TABLE  I 

SLOPE  COMPATIBILITY  RELATIONS  FOR  A  NINE -BAY  STIFFENED  PANEL 

VALUES  OF  C../  W 

i  j  mn 


i  =  1 

;  -  2 

i  -  3 

i  =  1 

(-l)m  n(a]b)/a2b2) 

(-i)n(b,/b2) 

(-l)m+n(o1b|/o2b2) 

i  “  3 

HAo/oj) 

1 

<-lf(0l/a2) 

i--3 

(-l)m+n(o1b)/a2b2) 

H)n(b,A2) 

(-l)m+n(a1b1/a2b2) 

Expressions  for  the  strain  energy  of  the  cover  sheet  os  a  function  of  the  temperature 
increase  con  be  developed  by  using  the  results  from  Equation  (2)  for  temperatures  below 
the  critical  temperature,  and  from  Equations  (19)  and  (23)  for  temperatures  above  the 
critical  temperature. 


For  a  panel  bay  with  dimensions  a  x  b  and  thickness  h,  the  expression  for  the  strain  energy 
of  the  fundamental  mode  as  a  function  of  the  temperature  ratio,  r,  is 


U 


-'.I*?! 

I  , 


(0«rN) 


(34a) 
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a)  STIFFENER  GEOMETRY  AND  NOTATION 
PARALLEL  TO  Y-AXIS 
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2 

|  'I  o  \ 

\  h  /  J 

q!i 

where 


(34b) 


F11(b.,a.)  =  (b./a.)  +  (a./bO 

R..  =  3((5  -  v2)F2  (b.,o.)  -  2(5  r  v)(!  -  v)] 

'I  n  I  ' 

The  slope  compatibility  coefficients  can  now  be  used  to  obtain  the  structural  response  in 
terms  of  the  amplitude  of  the  center  bay.  It  is  also  necessary  to  scale  the  temperoture 
parameter,  r j j ,  since  the  critical  temperoture  is  a  function  of  the  panel  area  (see  Equotion 
(5)) .  The  scale  factors  on  r  are  referenced  to  the  critical  temperature  of  the  center  boy. 

It  is  now  convenient  to  introduce  the  notation 


r..  -  d..r 
'I  '  I 


d.. 

‘I 


°ibif1l(b2'°2) 


..w  =  W  C . .  q..(t)  =  q(t)c.. 

I|  O  o  l|  i|  i| 


(35) 


where  r  is  the  critical  temperature  ratio  of  the  center  bay,  d;t  ore  the  scale  parameters 
for  the  temperature  rise,  C;j  are  the  slope  compatibility  coefficients,  ond  WQ  and  q(t) 
represent  the  static  buckling  and  vibration  amplitudes,  respectively,  of  the  fundamental 

mode  of  the  center  bay.  (The  constants  C..  are  defined  in  Table  I  .) 

•I 

Introducing  this  notation  into  Equotion  (34a)  and  summing  ali  terms  for  each  of  the  nine 
panel  bays  (see  Figure  5),  the  expression  for  the  total  strain  energy  of  the  cover  sheet  for 
temperatures  below  the  critical  temperature  becomes 


U 

P 


h4D 

^?2 


F*(r)q‘ 


(0  £  r  i  r  ) 

o 


(36) 


where 


F*(r)-F^(b2,a2)(1  -  r)  -  2^  jF^b,  ,o2)(i  -  d2Jr) 

+ -v”  4(^)f>  -dn'» 
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The  parameter  rQ  ij  the  -,alue  of  the  temperature  ratio  for  which  the  strain  energy  is  zero, 
and  will  be  evaluated  later.  For  r  =  0  (ambient  temperature)  this  expression  is  identical 
to  Equation  (17)  of  Reference  8. 


Likewise,  introducing  the  notation  of  Equation  (35)  into  Equation  (34b)  ond  summing  ail 
terms  for  each  of  the  nine  panel  bays,  the  expression  for  the  strain  energy  of  the  cover 
sheet  for  temperatures  above  the  critical  temperature  becomes 


R..  =  3l(5  -  v2)F^.(b.,a.)  -  2(5  + v)(l  -v)j 
H  I'  I  i 

Equation  (36)  represents  the  strain  energy  of  the  heated  skin  for  temperotures  below  the 
critical  temperature  while  Equation  (37)  is  the  strain  energy  of  the  heated  skin  for 
temperatures  above  the  critical  temperature. 


Frequency  Response 


The  equation  of  motion  for  the  fundamental  mode  of  the  nine-bay  panel  is  obtained  by 
summing  the  strain  energy  terms  and  the  kinetic  energy  terms  for  the  cover  sheet  and  the 
stiffeners,  and  then  applying  Lagrange's  equation. 


Since  the  kinetic  energy  is  independent  of  the  temperature  increase,  this  result  will  again 
be  presented  first.  The  kinetic  energy  for  the  fundamental  mode  of  the  nine-boy  panel 
structure  is,  from  Equations(32)  and  (33), 


Below  the  critical  temperature  (0  £  r  <  rQ),  the  strain  energy  for  the  fundamental  mode  of 
the  nine-bay  panel  is 
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(39) 


U  CF*(r)  +  K  +  K  ]W2 

80202  sx  s9  22 


(0  £  r  *  r  ) 


where  Equations  (31)  and  (36)  hove  been  used  to  obtain  this  result,  ond  again, 

value  of  r  such  that  U(r  )  =  0. 

o 

The  substructure  stiffness  terms  in  the  preceding  equation  are 


4n2Er* 

K  -  — v~ 

sx  n  2u 

Da2b2 


4r2EP 
K  =  - 

S/  °a2b2 


,tkx2t2(;' 


,+  kV2(r. 


<14  kx|) 


0  +  k- \) 


Equations  (38)  and  (39)  may  be  reduced  to  the  following  conventional  forms 
T  =?M1  ^22  °nd  U  =  IKl(r)  W22  <0*r*r) 


where  the  mass  and  stiffness  are  defined  by 

.3 


V? 


Vha2b2 


,-4)-(S)3 


2/ 


hr 


2  2 
+  4n  vl  4 
P*  2 
b2  u 


]  +  2|  — 


2  2 
+  4n  Yl  -4 
py  2 

a2 


1  +  24 


4 

K,(r)  =  xrrr  +  K  +K  1 
1  4<32b2  s*  sy 

The  equation  of  motion  of  the  nine-bay  panel  below  the  critical  temperature  i< 
free  vibration, 

M  q  +  K  ,(r)q  =  0  (O^rsr  ) 

I  i  o 

and  assuming  simple  harmonic  motion,  the  response  frequency  is  given  by 

f(r)  [K|(r)/M1]1/2  (OsrSrj 


r  is  the 
o 

(40) 

(41) 

(42) 

(43) 

(44) 

then,  for 

(45) 

(46) 
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for  o  temperoture  increose  above  the  critical  temperature,  rhe  k'netic  energy  expression, 
Equation  (39),  is  still  valid  and  Equations  (31)  and  (37)  are  used  to  obtoin  the  expression 
for  the  strain  energy.  As  indicated  in  the  development  of  Equation  (37)  the  following 
approximations  are  assumed 

2  2  2 

W72*Woq 

U  0  (47 

4  4  2  2 

W,_*»W  t  6W  q 
22  o  a^ 

where  W0  is  the  buckling  amplitude  of  the  center  bay  and  q(t)  is  the  dynamic  amplitude 
of  the  center  bay  for  the  fundamental  mode  (see  Equation  23). 


It  is  evident  that  the  kinetic  energy  expression  is  then  simply 

T-iM}q2  (r*ro)  (43) 

where  is  defined  by  Equation  (43)  since  the  moss  is  independent  of  temperature  . 

The  strain  energy  of  the  nine-bay  panel  is  obtained  by  using  Equations  (31)  ond  (37),  and 
is 

I T  r  ol  1 


=  fr'l!  ■F*(r)  +  §R‘(ir)  4K  4K 

^2^2  LL  8\h/J  sx  sy 


r 
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L 

Lf*(r) 

+  K  +  K 

SX  Sy 

Setting  q(t)  =  0  and  minimizing  the  result  with  respect  to  W q, 


n4DW 

cU  o 

?W  8o_b„ 
o  2  2 


*(r)  **  2K  +  2K  -  0 

sx  sy  2  \  h  / 


Then  neglecting  the  trivial  case  of  W  -  0,  the  displacement  ratio  (W  /h)  is  obtained  as 

o  o 

/w  •?  , 

Itc)  V  (,8ro>  <5 

"here  F*(r)  is  defined  by  Equation  (36)  and  R*  is  defined  by  Equation  (37)  . 

Applying  Lagrange's  equation  to  the  kinetic  and  strain  onergy,  the  equation  of  motion  for 
the  fundamental  mode  of  the  buckled  nine-bay  panel  is  derived  for  free  vibrotion  and 
takes  the  conventional  form 
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(52) 


M 


jq  +  K2(r)q  -  0 


where  the  mass  was  identified  previously  and  the  stiffness  is 


y,  (  V  _  ”4P 

K2(r)  "  "2^T2 


[F*(r)  +  K  +  K  ] 

sx  sy 


(53) 


Comparison  of  the  expressions  for  K|(r)  and  l<2(r)  reveals  that  both  vanish  for  the  some 
value  of  r,  which  is  to  be  expected.  Indeed,  the  result  simplifies  to 

«2(r)  -  -  2K  ^(r) 


(54) 


K2(ro>=Kl(ro)=0 


The  response  frequency  is  determined  from  Equation  (50)  to  be 


f(r) 


(r*rj 


2n 


(55) 


The  results  presented  to  this  point  may  be  simplified  by  returning  to  Equations  (35)  and  (36), 
Using  the  defin'tion  for  d;j  and  F|  ] (b: , a;)  the  plote  stiffness  parameter  may  be  expressed  as 


where 
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Finally,  the  value  of  r  at  which  the  strain  energy  of  the  structure  becomes  zero  is  derived 
from  Equation  (39),  and  is 


r 

o 


'F|(fi* 


K 

sx 


(57) 


Using  this  definition  of  r 


fhe  previous  results  for  stiffness  may  be  further  simplified  to 
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Ill  -  EXPERIMENTAL 


This  overall  program  was  primarily  experimental,  with  the  preceding  analysis  serving 
as  the  means  of  establishing  parameters  to  be  measured  during  the  test  program  for 
later  correlation  with  the  analytical  results.  The  experimental  program  was  conducted  in 
two  phases:  (1)  vibratory  fatigue  tests  of  coupon  specimens,  and  (2)  acoustic  fatigue 
tests  of  multi-bay  stiffened  panels.  Both  series  of  tests  were  conducted  at  room  and 
elevated  temperatures  with  specimens  fabricated  of  7075-T6  aluminum  and  6AI-4V  titanium 
(mill-annealed)  alloys. 

Coupon  vibratory  fatigue  tests  were  conducted  to  generate  basic  material  fatigue  properties 
for  each  alloy  ar  room  and  elevated  temperatures  in  the  absence  of  mean  (thermal) 
stresses.  These  data  were  used  in  the  design  of  stiffened-skin  specimens  which  were  fatigue 
tested  under  random  amplitude  acoustic  loading.  The  stiffened  panel  test  program  provided 
static  (thermal)  and  dynamic  test  data  for  correlation  with  the  analytical  results, 

A.  Coupon  Fatigue  Tests 

Fatigue  tests  were  conducted  on  cantilever  beam  specimens  to  develop  random  loading 
fatigue  properties  at  room  and  elevated  temperatures.  The  two  alloys  were  tested  at 
the  following  temperatures: 

7075-T6  Aluminum  -  Room  temperature  and  300°F 
6AI-4V  Titanium  -  Room  temperature  and  600°F 

Cantilever  beam  specimens  were  used  to  develop  zero  mean  stress  fatigue  data  and  to 
evaluate  the  thermal  degradation  effects  on  each  alloy.  The  plain  (unnotched)  and 
riveted  specimens  are  detailed  in  Figure  7.  The  riveted  specimens  incorporated  two 
rivets  to  represent  the  stress  concentration  in  the  stiffened-skin  test  specimens. 

All  specimens  were  fabricated  from  close  tolerance  sheet  using  standard  manufacturing  and 
quality  control  processes  and  procedures.  Each  specimen  was  sheared  oversize  and  milled 
to  the  final  size  to  remove  deformed  or  compressed  material  from  the  edges  caused  by  the 
shearing  operation.  All  holes  were  located  relative  to  the  free  end  of  the  beam  and  jig- 
drilled  to  provide  close  control  over  specimen  cantilever  length.  Rivets  were  installed  by 
a  Drivematic  riveting  machine  to  give  uniform  rivet  installation. 

1 .  Test  Procedure 


Three  specimens  were  fatigue  tested  simultaneously  in  each  set-up,  using  random  amplitude 
vibratory  excitation.  The  room  temperature  test  set-up  is  shown  in  Figure  8.  The  specimens  were 
supported  by  steel  clamp  blocks  (detailed  in  Appendix  III)  ottached  to  the  table  of  on  MB 
E lectronics  C-1 OE  electro-mechanical  shaker.  Phenolic  inserts  were  used  for  room  temperature 
tests,  while  steel  inserts,  with  asbestos  insulation  on  each  side  of  the  specimen,  were  used  for 
all  elevated  temperature  tests  to  minimize  heat  flow  from  the  specimen  to  the  clamp  blocks. 
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1 .00 


a)  PLAIN  SPECIMEN 


A  FLUSH  HEAD  RIVET,  D  =  .  125 


b)  RIVETED  SPECIMEN 


FIGURE  7.  COUPON  SPECIMEN  CONFIGURATION 
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The  room  temperature  set-up  was  duplicated  for  the  elevated  temperature  tests,  with 
tungsten-filament  quartz  lamps  and  a  thermal  enclosure  added.  Two  quartz  lamps  were 
located  above  the  specimens  to  provide  heat,  as  shown  in  Figure  9.  The  specimens  and 
lamps  were  contained  in  an  enclosure,  shown  in  Figure  10,  to  maintain  a  uniform 
thermal  environment  around  the  specimens.  A  slot  in  the  front  of  fhe  enclosure  allowed 
visual  observation  ^f  the  specimen  response  at  the  elevated  temperatures. 

Specimen  temperatures  were  monitored  by  a  thermocouple  attached  to  a  1/2  inch  wide 
strip  of  the  specimen  alloy  rigidly  fixed  between  adjacent  specimens  at  the  clamp  line. 

Strain  levels  were  measured  with  uniaxial  strain  gages  bonded  to  the  specimen  at  the 
locations  shown  in  Figures  7  and  11.  Generally,  only  the  center  specimen  yyos  instrumented 
since  specimen  responses  were  essentially  identical.  Room  temperature  curing  adhesives 
were  used  to  bond  all  gages. 

Excitation  levels  were  monhored  hy  means  of  an  accelerometer  attached  directly  to  the 
shaker  table. 

a.  Strain/Acce leration  Calibration  -  The  strain  gage  fatigue  life  was  limited  by  the 
high  acceleration  forces  experienced  during  resonance  testing.  Also,  the  gages  could  only 
be  used  at  ambient  temperatures,  since  a  room  temperoture  curing  adhesive  was  used  to  bond 
the  gages  to  the  specimens.  This  required  establishment  of  a  strain/acceleration  transfer 
function  for  use  in  maintaining  specimen  strain  levels  after  gage  failure. 

Narrow-band  random  input  spectra  were  used  for  all  tests,  with  a  typical  bandwidth  of 
30  Hz.  The  specimen  resonances  were  located  just  below  the  upper  cut-off  frequency  to 
maintain  uniform  inpuf  levels  as  the  resonant  frequency  decreased  during  failure. 

Simultaneous  overall  strain  response  and  input  acceleration  levels  were  measured  at 
various  input  levels,  using  the  test  spectrum  described  above.  These  data  were  plotted 
as  shown  in  Figure  12  to  produce  an  excitation  versus  response  curve  for  each  set  of 
specimens.  The  acceleration  level  then  become  the  control  parameter  for  maintaining 
constant  strain  levels  for  the  duration  of  the  test. 

2. _ Fatigue  Tests 


The  specimens  were  fatigue  tested  at  the  appropriate  test  temperature  with  the  excitation 
spectrum  used  for  the  strain/acceleration  calibration.  The  specimen  temperature  was 
allowed  to  stabilize  prior  to  application  of  vibration  for  the  elevated  temperature  tests. 
Periodic  inspections  were  performed  during  testing  to  detect  initiation  of  fatigue  failure. 
During  these  inspections,  the  test  was  stopped,  and  the  random  excitation  source  was  replaced 
with  sinusoidal  motion  in  order  to  accurately  determine  the  icsonant  frequency.  The  speci¬ 
men  temperature  was  maintained  during  these  inspections.  A  visual  inspection  was  also 
performed  on  the  room  temperature  specimens  at  this  time. 

Fatigue  failure  was  defined  as  the  time  at  which  the  resonent  frequency  decreased  2%  below 
its  initial  value.  This  was  determined  from  a  plot  of  frequency  versus  test  time  for  each 
specimen,  os  shown  in  Figure  13. 
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FIGURE  10.  ELEVATED  TEMPERATURE  TEST  THERMAL  ENCLOSURE 


FIGURE  12.  STRAIN/ACCELERATION  CALIBRATION  CURVE 


FAILURE 


*H  -  ADN30O3yd  3SNOdS3d 


-  ADN3n03dd  3SNOdS3d 


FIGURE  13.  TYPICAL  FREQUENCY  TIME  HISTORIES  FOR  COUPON  SPECIMENS 


Figures  l4and  15  are  typicai  probability  density  ::na!y>c3  oi  t 're  acceleration  excitation  and 
strain  response  signals  for  two  response  amplitudes.  As  indicated  in  the  figures,  these 
data  approximate  a  Gaussian  distribution. 

3. _ Test  Results 

The  coupon  fatigue  data  are  tabulcted  in  Tables  II  and  III  while  Figures  16  and  17  present 
the  fatigue  curves  developed  from  these  data.  The  median  life  is  shown  by  a  least  squares 
regression  line  drown  through  the  data  points  for  each  configuration  and  test  condition. 

Tire  scatter  present  in  the  data  is  indicated  by  the  95%  confidence  bonds  shown  on  the 
plots.  The  regression  line  and  confidence  bands  were  computed  by  the  methods  delineated 
in  Reference  4,  Appendix  II.  The  results  of  the  statistical  analyses  conducted  on  these 
data  are  shown  in  Table  IV  tor  reference.  The  fatigue  curve  for  riveted  aluminum  at 
room  temperature  was  found  to  agree  very  closely  with  similar  data  published  in  Reference  4. 

B.  Stiffened  Pane1  Fatigue  Tests 
1.  Test  Specimen  Design 

The  fatigue  data  from  the  coupon  tests  were  used  to  modify  the  stiffened  structure  skin 
design  nomograph,  Figure  73  of  AFFDL-TR-67-1  56^.  The  sti  ffened-sk  in  test  panels  were 
designed  using  this  nomograph, with  temperature  effects  included,  and  the  following 
guidelines: 

q 

o  Overall  size  identical  to  the  test  panels  described  in  AFFDL-TR-7)  - 107 
to  use  existing  frames. 

o  Range  of  structural  and  test  parameters  consistent  with  those  of  Ar FDL- TR-67- 1  56  and 

AFFDL-TR-71 -1  07  to  provide  a  logical  transition  from  ambient  to  High  temperature 
design  criteria . 

o  Total  quantity  of  27  aluminum  and  21  titanium  specimens. 

Previous  test  programs  '  °  Have  shown  that  fatigue  failures  could  be  expected  in  the  outer 
boys  ir  all  bays  were  of  the  same  approximate  size.  Therefore,  the  center  bays  of  all  spec¬ 
imens  were  designed  to  have  a  greater  surface  area  than  that  of  any  of  the  adjacent  outer 
bays  since  the  center  bay  was  considered  to  be  the  primary  test  area. 

Eigr.t  aluminum  and  four  titanium  specimens  were  designed  to  be  tested  at  room 
temperature  for  direct  comparison  with  the  data  of  Reference  4.  An  identical  number 
biese  same  specimen  designs  were  to  be  tested  at  eievoted  temperatures.  Another  eight 
c1  the  aluminum  specimens  were  also  designed  with  faying  surface  sealant  to  evaluate  the 
et'tect  of  the  sealant  on  fatigue  endurance;  these  specimens  were  identical  to  eight  of  the 
specimens  without  sealant.  The  specimens  with  sealant  were  included  In  this  program  to 
ascertain  trie  validity  of  previous  empirical  design  cri  teria ,  whichare  based  on  tests 
conducted  on  specimens  without  sealant,  for  use  on  structures  which  use  sealant. 
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FIGURE  14.  AMPLTUDE  DISTRIBUTION  OF  EXCITATION 
Til  ANIUM  COUPON-  SPECIMEN  (RIVETED) 


35 


2 


-2-1  0  1 

RATIO  OF  INSTANTANEOUS  TO  RMS  AMPLITUDE 


b)  OVERALL  STRAIN  LEVEL:  480  m  IN/IN  RMS 

URL  15.  AMPLITUDE  DISTRIBUTION  OF  SPECIMEN  RESPONSE 
TITANIUM  COUPON  SPECIMEN  (RIVETED) 
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TABLE  II 

ALUMINUM  COUPON  SPECIMENS 
SUMMARY  OF  FATIGUE  TEST  RESULTS 

Ploin  Specimens 


Room  Temperature  ] 

1  300°F  I 

Specimen 

No. 

Stress 
ksi  rms 

m 

es 

Specimen 

No. 

Stress 
ksi  rms 

IBH 

1 

22.4* 

1.08  x 

104 

1 

18.5 

1 . 14  x  I05 

2 

22.4* 

6.07  x 

103 

d. 

18.5 

9. 75  x  104  . 

3 

22.4* 

9.00  x 

103 

3 

18.5 

7.42  x  104 

4 

17.6 

2.96  x 

105 

4 

17.0 

2.87  x  105 

1  5 

16.8* 

2.31  x 

105 

5 

17.0 

2.87  x  I05 

6 

16.8* 

1.43  x 

105 

6 

17.0 

2.88  x  105 

7 

16.8* 

1.01  x 

105 

7 

10.6 

1 . 80  x  1 06 

8 

16.8 

6.14  x 

10f 

8 

10.6 

1 .46  x  106 

9 

14.6 

6.80  x 

105 

9 

10.6 

1 . 93  x  1 06 

10 

14.0 

8.26  x 

105 

10 

10.2 

2. 00  x  106 

11 

14.0 

1.53  x 

106 

11 

10.2 

2. 00  x  106 

12 

14.0 

8.80  x 

to5 

12 

8.5 

2.60  x  106 

Riveted  Specimens 


Room  Temperature  1 

300°F 

Specirn* 

No. 

jn  Stress 

ksi  rms 

1  miBsm 

1 

14.0 

1 . 27  x  1 05 

1 

14.7 

4.  40  x  104 

2 

14.0 

1.31  x  105 

2 

14.7 

4.30  x  104 

3 

14.0 

1 .30  x  105 

3 

14.7 

4.08  x  104 

4 

12.0 

4.97  x  105 

4 

9.2 

2.55  x  105 

5 

12.0 

5.12  x  105 

5 

9.2 

3.00  x  105 

6 

12.0 

1.15  x  105 

6 

9.2 

2. 90  x  105 

7 

10.5 

8. 20  x  105 

7 

6.2 

1.22  x  106 

8 

10.5 

7.92  x  105 

8 

6.2 

1.51  x  106 

9 

10.5 

'.20  x  105 

9 

6.2 

1.52  x  106 

10 

6.6 

6.75  x  106 

10 

4.1 

4.  44  x  10^ 

11 

5.9 

4.11  x  106 

11 

4.1 

5.22  x  106 

12 

5.7 

2. 70  x  106 

12 

4.1 

7.82  x  106 

*  Invalid  data  -  ends  of  beam  weighted  arid  response  highly  nonlinear.  Not  used  in 
fatigue  curve . 

NOTE:  Fatigue  life  (in  cycles)  is  the  product  of  resonarn  frequency  and  test  time  to  fail 


TABLE  III 

TITANIUM  COUPON  SPECIMENS 
SUMMARY  OF  FATIGUE  TEST  RESULTS 


Plain  Specimens 


Room  Temperature 

600°F 

Specimen 

■ 

Spec imen 

Stress 

Life 

No. 

■m 

No. 

ksi  rms 

Cycles 

1 

35.5 

3.28  x  105 

1 

26.9 

5. 40  x  105 

2 

35.5 

3. 24  x  105 

2 

26.9 

2.74  x  105 

3 

35.5 

3.55  x  10-j 

3 

26.9 

3.38  x  105 

4 

28.7 

5. 15  x  105 

4 

25.3 

2.74  x  105 

5 

28.7 

9. 45  x  105 

5 

24.8 

2.50  x  105 

6 

28.7 

8. 65  x  105 

6 

24.6 

2.26  x  105 

7 

21.4 

5.18  x  10* 

7 

22.9 

3.38  x  105 

8 

21.4 

3.97  x  106 

8 

22.2 

2.26  x  105 

9 

21.4 

2.57  x  106 

9 

21.5 

2.82  x  105 

10 

21.2 

1 .89  x  106 

10 

19.7 

8.90  x  105 

11 

21.2 

1 .67  x  106 

11 

19.7 

8.38  x  105 

12 

21.2 

5.52  x  106 

12 

19.7 

1.14x106 

13 

17.9 

1.53  x  107 

13 

14.9 

7.18  x  106 

14 

17.0 

3.05  •:  107 

14 

14.9 

7.32  x  106 

15 

14.9 

5.00  x  106 

Riveted  Specimens 


Room  Temperature 

600°F  i 

Stress 

''  j  l' 

Stress 

flEEHi 

ksi  rms 

mzssmm 

ksi  rms 

i 

20.9 

5. 00  x  105 

i 

19.2 

3.48  x  104 

2 

19.4 

4. 1 8  x  105 

2 

19.2 

5.34  x  104 

3 

17.8 

3.82  x  105 

3 

19.2 

5.93  x  104 

4 

16.0 

1.11  x  1 06 

4 

13,6 

4.57  x  105 

5 

16.0 

9.60  x  105 

5 

» 3.6 

2.52  x  105 

6 

16.0 

4.50 x  105 

6 

13.6 

5.02  x  105 

7 

13.9 

2.23  x  106 

7 

10.  1 

8.55  x  105 

8 

13.9 

1 .34  x  106 

8 

10.1 

6.24  x  105 

9 

13.0 

1 .46  x  106 

9 

10.1 

1 . 04  x  1 06 

10 

12.8 

4.10  x  106 

10 

5.5 

1 .65  x  107 

11 

12.0 

3.02  x  10j 

11 

5.2 

8.02  x  106 

12 

11.2 

5.22  x  106 

12 

3.5* 

1 . 87  x  1 07 

*Did  not  fail. 

NOTE:  Fatigue  life  (in  cycle?)  is  the  product  of  resonant  frequency  and  test  time  to  failure. 
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RANDOM  AMPLITUDE  REVERSED  BENDING  -  ZERO  MEAN  STRESS 


o  PLAIN  SHClMINS  - - if  A\f  SQUAMS  ■tGMS'SlON  UNI 

£  i; .  i  I  Ur  '.PiCJMf  - -  --  liMir.  Cf  *»S  CONflOl  K  I  MUi 


o)  ROOM  TEMPERATURE 


FIGURE  17.  COUPON  FATIGUE  TEST  RESULTS 
TITANIUM  SPECIMENS 
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TABLE  IV 

SUMMARY  OF  STATISTICAL  PROPERTIES 
FOR  COUPON  FATIGUE  DATA 


Test 

Sample 

Size 

Regression  Line* 

°c 

Rc 

Al  loy 

Temperature 

°F 

A 

Intercept 

B 

Slope 

Standard 

Deviation 

Correlation 

Coefficient 

Aluminum- 

Plain 

RT 

6 

2.424 

-.215 

** 

** 

Aluminum- 

Plain 

300°  F 

12 

2.362 

-.215 

0.034 

-0.969 

Aluminum- 

Riveted 

RT 

12 

2.261 

-.217 

0.059 

-0.921 

Aluminum- 

Riveted 

300PF 

12 

2.376 

-.260 

0.020 

-0. 996 

Titanium- 

Plain 

RT 

14 

2.438 

-.167 

0.039 

-0. 942 

Titanium- 

Plain 

60tf>F 

15 

2.234 

-.155 

0.039 

-0.915 

Titanium- 

Riveted 

RT 

12 

2.335 

-.190 

0.035 

-0.915 

Titanium- 

Riveted 

600°  F 

1 1 

2.377 

-.230 

0.045 

-0.976 

NOTE:  ^Regression  Line:  Log  c1  =  A  +  B  Log  N  (a  in  ksi) 

**lnsuffic  ient  data  points  -  regression  line  slope  token  same  as 
300^F  condition. 
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The  general  specimen  configuration  is  shown  in  Figure  18,  with  the  detailed  dimensions 
for  each  of  the  panels  listed  in  Tables  V  and  VI  for  the  aluminum  and  titanium  alloys, 
respectively.  More  detailed  drawings  are  presented  in  Appendix  III.  As  noted  in  Table 
V,  two  of  the  aluminum  specimen  designs  had  12  panel  bays  because  of  the  small  center 
bay  in  relation  to  the  overall  specimen  size.  All  other  specimens  had  nine  panel  bays. 

Two  replicates  of  each  specimen  design  were  fabricated  to  provide  statistical  accuracy 
consistent  with  previous  test  programs.  One  additional  titanium  specimen  was  fabricated 
for  an  investigation  of  thermal  buckling  effects  and  was  subsequently  fatigue  tested. 

An  additional  aluminum  specimen  was  fabricated  to  evaluate  the  effect  of  intermediate 
temperature  on  life. 

Aii  specimen  manufacturing  was  accomplished  in  conformance  with  standard  aircraft 
manufacturing  processes  and  procedures,  with  qualify  control  surveillance  at  all  steps. 

All  aluminum  zee-section  stiffeners  were  hot  formed  from  sheet  stock,  while  extrusions 
were  used  for  the  frames  and  clips.  All  titanium  ports  were  chemically  cleaned  during 
manufacture  to  prevent  degradation  by  foreign  materials.  The  titanium  stiffeners  and 
frames  were  fabricated  from  sheet  stock  by  welding  strips  together  to  form  the  desired  sec¬ 
tions.  An  automatic  Tungsten-Inert-Gas  (TIG)  welding  process  was  used  for  these 
longitudinal  welds.  These  parts  were  subsequently  straightened  by  stress  relieving.  This 
approach  was  taken  because  of  the  prohibitive  cost  of  procuring  extrusions  or  hot-forming 
these  titanium  sections  in  small  quantities.  Cold-forming  was  deemed  undesirable  in 
most  instances  because  of  the  large  bend  radii  that  would  result  (at  least  six  times  the 
thickness).  Automatic  riveting,  with  a  Drivmatic  riveting  machine,  was  used  whenever 
possible  for  uniformity.  Monel  countersunk  rivets  (MS20427M4)  were  used  on  the  titanium 
specimen  skins,  while  aluminum  countersunk  rivets  (LS10795  MP4)  were  used  on  the  skins 
of  all  aluminum  specimens.  Typical  aluminum  and  fitonium  specimens  are  shown  in 
Figure  19. 

2h _ Test  Set-up 

The  test  panels  were  installed  in  a  steel  frame,  with  two  identical  specimens  mounted 
side-by-side,  as  shown  in  Figure  20.  Two  steel  frames  were  used  to  minimize  down  time 
between  tests. 

All  tests  were  conducted  in  the  high  intensity,  grazing  incidence,acoustic  tert  facility  shown 
in  Figure  21  .  The  test  frame  was  installed  in  the  facility  wall  with  the  panel  substructure 
exposed  to  the  noise,  opposite  to  the  normal  installation,  to  enable  application  of  heat 
directly  to  the  skin  surface  without'  subjecting  the  lamps  to  the  noise  environment 
(see  figure  22).  Response  tests  conducted  at  the  beginning  of  the  test  program  indicated 
no  discernible  difference  in  panel  response  due  to  panel  orientation  (i.e.,  substructure 
inside  or  outside  the  test  chamber). 

Elevated  temperatures  were  attained  with  o  bank  of  2500  watt  tungsten-filament  quartz 
lamps  mounted  vertically  in  front  of  the  test  panels.  The  3/8-inch  diameter  quartz  lamp 
elements  were  encased  in  a  1 -inch  quartz  tube,  with  the  lamp  ends  isolated  from  the  tube 
for  protection  from  the  noise  environment.  The  lamp  installation  is  shown  in  Figure  22. 

The  lamps  were  powered  by  c  multi-channel  power  supply,  with  two  lamps  per  channel  to 
provide  uniform  temperature  control  over  the  surface  of  the  panel .  Temperatures  were 
contiolled  by  varying  power  input  through  a  silicone  controlled  rectifier  on  each 
channel . 
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PLAN  VIEW 

NOTE:  EDGE  MEMBER  OMITTED  FOR  CLARITY 
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FIGURE  18.  STIFFENED-SKiN  SPECIMEN  CONFIGURATION 
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TABLE  V 

ALUMINUM  STIFFENED-SKIN  SPECIMEN  DETAILS 


Specimen  1  No.  of 

No.  of 

1  Specimen  Dimensions 

J 

Test 

No. 

Specimens 

Bays 

b 

IB9 

h 

Temperature 

AL-1 

2 

12 

6 

9 

1.5 

0.032 

0.040 

300°  F 

-2 

2 

12 

5 

10 

2.0 

0.032 

0.040 

300°  F 

-3 

2 

9 

6 

12 

2.0 

0.032 

0.040 

RT 

-4 

2 

9 

6 

12 

2.0 

0.032 

0.040 

300°  F 

-5 

2 

9 

6 

12 

2.0 

0.050 

0.063 

300°  F 

-6 

2 

9 

9 

18 

2.0 

0 .040 

0.050 

RT 

-7 

2 

9 

9 

18 

2.0 

0.040 

0.050 

300°F 

-8 

2 

9 

9 

18 

2.0 

0.063 

0.071 

300°F 

-9 

2 

9 

6 

18 

3.0 

0 .032 

0.040 

300°F 

-10 

/h  2 

9 

6 

12 

2.0 

0.032 

0.040 

RT 

-11 

A  2 

9 

6 

12 

2.0 

0.032 

0.040 

300°F 

-12 

A  2 

9 

9 

18 

2.0 

0.040 

0.050 

RT 

-13 

A  2 

9 

9 

18 

2.0 

0.040 

0.050 

300°F 

-14 

1 

9 

6 

12 

2.0 

0.032 

C  .040 

150°F 

TOTAL: 


NOTES: 


1 .  h  -  skin  thickness  -  inch 

2.  hr  =  stiffener  thickness  -  inch 

3.  Overall  panel  dimensions  21 . 5  x  31.5  inches. 
/ff\  Seal  per  LAC  G230  using  STM40-1 1 1  sealant. 


TABLE  VI 

TITANIUM  STIFFENED-SKIN  SPECIMEN  DETAILS 


Specimen 

No,  of 

Specimen  Dimensions 

Test 

No. 

Specimens 

a 

L± 

II2QHI 

n 

hr 

Temperature 

1 

1 

6 

12 

2.0 

0.024 

0.036 

600°  F 

2 

2 

6 

12 

2.0 

0.032 

0.044 

RT 

3 

2 

6 

12 

2.0 

0.032 

0.044 

400°  F 

4 

2 

6 

12 

2.0 

0.032 

0.044 

600°F 

5 

2 

6 

18 

3.0 

0.032 

0.044 

400°  F 

6 

2 

6 

18 

3.0 

0.044 

0.056 

RT 

7 

2 

6 

18 

3.0 

0.044 

0.056 

400°  F 

8 

2 

6 

18 

3.0 

0.044 

0.056 

600°  F 

r\ 

7 

2 

8 

16 

2.0 

C.  032 

0.044 

400°  F 

10 

2 

8 

16 

2.0 

0.044 

0.056 

400°  F 

11 

2 

8 

16 

2.0 

0.056 

0.068 

400°  F 

TOTAL: 


NOTES: 


1.  h  -  skin  thickness  -  inch 

2.  hr  =  stiffener  thickness  -  inch 

3.  Overall  panel  dimensions  21 . 5  x  31.5  inches. 

4.  All  panels  have  nine  bays. 


b)  TITANIUM  SPECIMEN 

FIGURE  19.  TYPICAL  STI FFE NED-SK I N  TEST  SPECIMENS 
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A  modol  investigation  ai  room  and  elevated  temperature  was  conducted  on  the  first  titanium 
specimen  to  determine  placement  of  the  strain  gages.  Figure  23  shows  the  .truin  gage 
locations  tor  the  aluminum  and  titanium  specimens.  Micro-Measurements  V/K-Q5- 1  25AD-350 
strain  gages  were  bonded  to  the  titanium  specimens  with  BLH  Electronic  EPY  600  adhesive 
cured  at  650°F.  The  1-inch  integral  gage  leads  were  silver  soldered  to  teflon  covered 
wires.  BLH  FAE-12-12S13ET  strain  gages  were  bonded  to  the  aluminum  pane*s  using  EPY 
600  adhesive  cured  at  35GPF.  Lead  wires  for  these  gages  were  30-gage  enameled  wire 
bonded  to  the  panel  surface.  Teflon  tape  was  used  over  the  lead  wires  and  gages  to  pre¬ 
vent  motion  of  the  wire.  Solder  tabs  at  the  edge  of  the  panel,  where  the  response  was  greatly 
reduced,  were  used  to  connect  the  enameled  wire  to  nn  insulated  lead  wire.  A  typical 
gage  installation  tor  the  aluminum  specimens  is  shown  in  Figure  24. 

Chrome l-Alume I  foil  thermocouples,  RdF  Corp.  20112,  wer;  bended  to  the  specimens  with 
EPY  600  adhesive  at  the  locations  shown  in  Figure  23  and  cured  simultaneously  with  the 
strain  gages.  The  thermocouples  were  located  an  the  basis  of  an  initial  thermal  survey  on 
the  first  titenium  specimen,  using  14  thermocouples.  The  locations  selected  for  each 
specimen  provided  skin  temperatures  in  the  vicinity  of  the  strain  gages.  Typical  thermo¬ 
couple  installations  a"  a. so  shown  in  Figure  24. 

3.  Test  Procedure 


Several  tests  were  conducted  cn  each  test  specimen  prior  to  fatigue  testing.  The  static  and 
dynamic  parameters  measured  were  those  indicated  by  the  analysis  to  be  requisite  to  deri¬ 
vation  of  useful  empirical  relations.  Instrumentation  schematics  for  these  tests  are  de- 
ta;lcd  in  Appiendix  III.  The  tests  were  generally  conducted  in  the  sequence  described  in 
tallowing  subsections. 

a.  Room  Temperature  Frequency  -  The  test  specimens,  installed  in  the  test  frame, 
were  mounted  over  due  ctro-mechanical  speakers  as  shown  in  Figure  25.  Low  level 
sinusoidal  noise  excitation  was  applied  to  the  specimen  and  strain  amplitude  plotted  versus 
frequency,  at  ambient  temperature,  for  eacl  strain  gage.  The  frequency  scale  near  the 
fundamental  mode  was  then  expanded  to  provide  betrer  resolution  of  this  mode,  and  the 
sine  sweep  was  repeated.  The  resultan*  ,train-h-<  quency  plot  was  used  to  compute  damping 
ratios  by  the  half-power  method^. 

0.  Temperature  Effects  or.  Frequency  -  The  rest  set-up  described  in  the  preceding 
sub  -St  cHort  was  also  'used  to  determine  the  changes  in  fundamental  frequency  with  in- 
e-  Ang  temperature.  Six  heat  lamps  were  mounted  above  the  specimens  as  shown  in 
Figure  26.  The  .ocm  tempe.ature  Chiadni  pattern  for  the  fundamental  mode  was 
obtained  vsii-  j  cor!-.  \  rtinlcs;  i!m  heat  lamps  we'e  rhe.i  turned  on  and  the  fundamental 
mode  frequency  was  tracker  by  manually  varying  *!'c  speaker  irput  freJUoncy  to  maintain 
c  fundamental  mode  response.  Skin  temperature  and  -otponse  frequency  vie  re  recorded 
on  data  seems  for  later  analysis 

The  critical  buckling  temperature  o'  th’  center  bay  was  a>so  determined  djmng  fnis  series  of 
♦ests.  Tni  ruckling  temperature  was  defined  as  the  tempe.  ‘  ji  incr  ase  ui  which  >Se  fre¬ 
quency  wv  •  a  minimum;  this  was  determined  tron-  a  plot  of  frequr  ncy  \  ■  sus  > emperofure 
fo-  each  •  uecimen. 
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JJJ — ’■  .V.  1  r.  : 
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TEMPERATURE 

SPECIMENS 
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•  THERMOCOUPLE 
-  STRAIN  GAGE 


<k 


BAY 


o)  ALUMINUM  SPECIMENS 


£  BAY 


b)  TITANIUM  SPECIMENS 


FIGURE  23.  STRAIN  GAGE  AND  THERMOCOUPLE  LOCATIONS 
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c.  _ Thermal  Strain  and  Deflection  -  Thermal  strain  and  deflection  of  the  panel 

center  bay  were  measured  on  all  elevated  temperature  specimens  in  the  absence  of 
acoustic  excitation.  The  test  panels,  in  the  test  frame,  were  installed  in  the  wall  of  the 
progressive  wave  test  chamber  with  the  heat  lamps  in  place.  Thermal  strains  were 
measured  with  a  BLH  SR4  Strain  Indicator  while  slowly  increasing  the  panel  temperature; 
the  set-up  for  these  measurements  is  shown  in  Figure  27.  A  compensating  strain  gage, 
identical  to  the  specimen  gages,  was  bonded  to  a  narrow  strip  of  the  specimen  alloy  and 
placed  on  the  face  of  the  panel  near  the  edge.  The  compensating  gage  location  was 
experimentally  determined  to  maintain  the  same  temperature  as  that  at  the  panel  gages. 

The  center  bay  displacement  due  to  thermal  buckling,  at  the  midpoint  of  the  bay,  was 
measured  concurrently  with  the  thermal  strains.  A  contacting  displacement  gage  was 
temporarily  clamped  to  the  frames  of  the  specimen  as  shown  in  Figure  28.  The  thermal 
strain  ond  displacement  data  were  recorded  on  data  sheets  versus  increasing  temperature 
for  later  analysis. 

d.  Excitation  Spectrum  Shaping  -  With  the  specimens  mounted  in  the  test  chamber, 
heat  was  gradually  applied  to  bring  the  specimens  to  the  desired  test  tempeiature. 
Sinusoidal  strain  response  plots  were  then  made  of  each  strain  gage  at  a  constant  sound 
pressure  level.  The  fundamental  response  frequency  at  the  test  temperature,  as 
determined  from  these  plots,  was  used  to  locate  the  frequency  limits  for  shaping  of  the 
input  noise  spectrum. 

The  input  noise  spectra  were  determined  by  the  type  of  response  obtained  during  these 
frequency  scans.  Where  single  mode  response  was  obtained  a  narrow-band  random 
input  with  a  100  Hz  bandwidth  was  used;  this  did  not  usually  require  additional  spectrum 
shaping.  Where  multi-mode  response  was  obtained,  o  broader  bandwidth  was  used 
(normally  300  Hz  bandwidth)  which  was  shaped  to  a  flat  respc  se  by  using  a  1/3  octave 
bond  spectrum  shaper. 

e.  _ Fatigue  Tests  -  Fatigue  testing  was  initiated  following  completion  of  spectrum 

shaping.  The  desired  temperature  was  stabilized  before  application  of  the  noise 
excitation.  The  test  excitation  level  was  established  by  using  the  modified  design 
chart  described  previously  (See  Section  III.B.l).  Temperatures  on  the  specimen 

were  recorded  on  a  multi-channel  strip  chart  recorder,  while  noise  and  strain  levels 
were  recorded  directly  onto  a  14-track  FM  magnetic  tape  recorder.  Parallel  outputs 
on  a  monitor  panel  allowed  direct  observation  of  the  noise  ond  strain  signals.  Test  levels 
were  controlled  by  mointaining  a  constant  overall  noise  level  at  one  of  the  microphone 
locations . 

Visual  inspection  of  each  specimen  wos  generally  accomplished  per  the  following 
schedule: 


Test  Time 
0  to  1  Hr. 

1  to  4  Hrs. 

4  Hrs.  to  End 


Inspection  Interval 
1  5  minutes 
30  minutes 
60  minutes 
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These  inspection  intervals  were  selected  to  minimize  down-time,  while  maintaining 
the  maximum  possible  resolution  on  life  (failure  detection),  since  a  cooling-heating 
cycle  was  necessary  for  each  inspection.  The  actual  inspection  intervals  varied  from  those 
shown  above  in  some  instances;  the  intervals  were  usually  shorter  (particularly  at  the 
beginning  of  the  test)  for  light  structure,  while  the  intervals  for  the  thickest  panels 
started  out  at  30  minutes. 

4.  Test  Results 


The  data  collected  during  the  various  tests  conducted  on  each  specimen  are  much  too 
extensive  to  be  included  in  their  entirety;  instead,  typical  examples  and  summaries  of  the 
test  data  are  presented,  in  ‘he  same  sequence  as  discussed  in  the  preceding  section. 

cr. _ Room  Temperature  frequency  -  Room  temperature  fundamental  mode  frequencies  are 

tabulated  in  Tables  VII  and  VIII  for  the  aluminum  and  titanium  specimens.  The  frequencies 
shown  were  recorded  at  the  ambient  temperature  (generally  65-70°F).  The  specimens  were 
allowed  to  stabilize  at  this  temperature  for  at  least  12  hours  prior  to  measurement  of  these 
frequencies;  hence,  the  panel  skin  and  substructure,  as  well  as  the  test  frame,  should 
have  been  in  a  neor-equiiibrium  thermal  state. 

Typical  examples  of  the  panel  strain  response  are  shown  in  Figures  29  (a)  and  30  (a) . 

These  show  the  dominant  fundamental  mode  response  for  the  selected  strain  gage  loca¬ 
tions,  which  is  normal  for  room  temperature  response.  Figure  31  presents  the  fundamental 
mode  strain  response  curves  for  this  specimen  with  an  expanded  frequency  scale  to  illus¬ 
trate  the  damping  ratio  measurement.  As  mentioned  previously  the  damping  ratios 
were  determined  by  the  half-power  method,  and  the  results  are  tabulated  in  Tables 
VII  and  VIII.  The  values  listed  are  averages  for  all  the  strain  gages  on  each  specimen. 

An  indication  of  the  damping  variation  between  configurations  is  obtained  from 
Figure  32,  which  shows  the  effect  uf  frequency  on  damping  ratio.  Data  from  four  other 
sources^'  8,  15,  16  are  also  shown  in  this  figure.  All  of  the  data,  except  for  that 
from  Reference  15,  are  for  fundamental  mode  response.  These  data  are  typicol  of  all 
damping  data  in  the  degree  of  scatter  present;  however,  they  may  prove  useful  in 
establishing  realistic  limits  of  assumed  damping  ratios  for  structural  design,  A  regres¬ 
sion  line  was  plotted  through  the  data  centroid,  as  shown  in  Figure  32,  with  a  slope 
taken  as  an  average  of  the  Reference  15  and  16  curves.  This  regression  line  can  be 
used  to  estimate  damping  for  all  modes,  including  the  fundamental.  As  an  alternative, 
the  mean  damping  value  of  0.016  (for  cl!  data  points  plotted)  may  be  used  for  the 
fundamental  mode. 

b.  Temperature  Effects  on  Frequency  -  The  change  in  the  fundamental  frequency 
with  increasing  and  decreasing  temperature  was  measured  for  empirical  correlation  and 
for  determination  of  the  panel  buckling  temperature.  Several  se's  of  data  were  accumulated 
for  most  of  the  specimens,  all  in  tabular  form  of  temperature  versus  frequency.  It 
was  found  that,  after  the  initial  heating  and  buckling,  during  cooling  the  unbuckling 
occurred  at  a  higher  temperature  than  the  initial  buckling,  as  shown  in  Figure  33.  This 
hysteresis  effect  was  attributed  to  the  gradual  heating  of  the  panel  substructure,  causing 
the  panel  to  reach  its  neutral  position  (unbuckling)  at  a  higher  skin  temperature  during  the 
decreasing  temperature  test.  Alternatively,  holding  the  panel  at  a  constant  temperature,  above 
buckling,  will  result  in  gradually  decreasing  buckling  amplitudes  due  to  expansion  of 


55 


TABLE  VII 

ALUMINUM  SPECIMENS 

SUMMARY  OF  FUNDAMENTAL  MODE  FREQUENCY  AND  DAMPING 
PANEL  CENTER  BAY  -  ROOM  TEMPERATURE 


Spec  imen 

Number 

Fo 

Fundamental  Mode 

Frequency  -  Hz 

i 

Fundamental  Mode 
Damping 

AL-1A 

217 

0.0145 

-IB 

225 

0.0080 

-2A 

233 

0.0055 

-2B 

244 

0.0097 

-3A 

166 

0.0135 

-3B 

184 

0.0160 

-4A 

194 

0.0188 

-4B 

171 

0.01 93 

-5A 

295 

0.0109 

-5B 

282 

0.0125 

-6A 

94 

0.0153 

-6B 

105 

0.01 55 

-7A 

76 

0.0203 

-7B 

80 

0.0187 

-8A 

144 

0.0095 

~8B 

141 

0.0192 

-9A 

146 

0.0167 

-9B 

163 

0.0252 

-10A 

215 

0.0139 

-1  OB 

187 

0.0186 

-1 1 A 

173 

0.0152 

-11 B 

148 

0.0192 

-12A 

94 

0.0177 

-12B 

78 

0.0179 

-13A 

86 

0.0195 

-I3B 

95 

0.0153 

-14A 

192 

0,02)2 
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TABLE  VIII 

TITANIUM  SPECIMENS 

SUMMARY  OF  FUNDAMENTAL  MODE  FREQUENCY  AND  DAMPING 
PANEL  CENTER  BAY  -  ROOM  TEMPERATURE 


Spec  imen 

Number 

f 

o 

Fundamental  Mode 
Frequency  -  Hz 

{ 

Fundamental  Mode 
Damping 

TI-1A 

172 

0.0084 

-2A 

210 

0.0176 

-2B 

227 

0.0275 

-3A 

211 

0.0094 

-3B 

191 

0.0156 

-4A 

201 

0.0167 

-4B 

200 

0.0125 

-5A 

171 

0.0172 

-5B 

166 

0.0204 

-6A 

205 

0.0169 

-6B 

178 

0.0154 

-7A 

205 

0.0155 

-7B 

200 

0.0115 

-8A 

195 

0.01 14 

-8B 

203 

0.0217 

-9A 

141 

0.0238 

-9B 

116 

0.0234 

-10A 

169 

0.0152 

-10B 

153 

0.0123 

-1 1 A 

190 

0.0213 

-1 1  B 

191 

0.0184 
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FREQUENCY  -  Hi 


FRfcQUENCY  -  H 7 


RELATIVE  STRAIN  AMPLITUDE  RELATIVE  STRAIN  AMPLITUDE 


DAMPING  RATI 


FUNDAMENTAL  MODE  FREQUENCY  ~  Hz 


the  substnicure.  It  should  be  noted  tfv.it,  if  the  entire  structure  (skin  and  lifft-ners) 
were  heated  uniformly  with  no  externa!  restraints  at  the  structure  boundaries,  the  skin 
wr  'd  not  buckle  because  ti.e  in-plane  residual  stresses  would  be  avoided. 

Panel  skin  buckling  temperatures  were  determined  from  the  above  frequency  •  .emperat  •  .'e 
plots  us  the  temperature  at  '-/Inch  the  frequency  was  a  minimum  (see  figure  3d).  Because 
of  the  h/steresi:  eff “C*,  only  the  initial  heating  cycle  was  used  to  determine  the  ui, skiing 
temperatu  es  since  heating  rates  and  substructure  temperatures  were  not  controlled 
Since  the  ana'yfical  buckling  temperatures  are  defined  as  a  temperature  increase  relat.ve 
to  ami  ien»,  the  measured  initial  temperature  was  subtracted  from  the  buckling  ter.iperc  ure. 
The  critical  buckling  temperatures  for  each  of  the  aluminum  and  titanium  specimen;  ate 
I i sherd  in  Yable  IX. 

c.  Thermal  Strain  and  Deflection  -  Thermal  strains  unci  cent'  ;  bay  buckling  ampli¬ 
tude  were  measured  orly  on  the  elevated  temperature  specimens,  since  rooi.i  temperature 
curing  udhas'ves  were  used  tor  the  strain  gages  ori  the  room  tpirpvrjiure  specimens,  The 
thernal  strain'  were  recorded  for  twe  purposes:  (1)  Determination  of  the  thermal  (moon) 
it'c.s  at  rhe  test  temperature,  and  (2)  correlation  w"  *h  analyficol  data  to  provide  em¬ 
pirical  ..quotions  for  predicting  thermal  strein.  Kigu'e  34  depicts  typical  thermal  strain 

■  jr'ction  with  temperature;  th ::  entire  ‘hermr,  I  strain  data  set  is  not  included  because 
of  its  .iz:  and  the  f aci  that  the  per'inci't  into.r.vjtion  is  contained  in  the  later  empirical 
noire  lotion,  Section  IV, 

Center  buy  buckling  amplitudes  wee  meosured  solely  for  correlation  with  the  cnoly- 
tiusl  results,  io  pros  ide  an  interim  correlation  in  the  thermal  stress  calculation.  Max¬ 
imum  buckling  ,«mpl itudes  at  the  individual  specimen  test  temperatures  are  itemized  in 
i able  X . 

d,  _ Excitation  Spec’rum  Shaping  -  The  results  of  the  sinusoidal,  high  intensity, 

noise  frequency  sweeps  ore  hypli fied  by  Figures  29  (b)  or.d  30  'b).  These  response  plots, 
made  at  the  individual  specimen  test  temperatures,  served  to  locate  the  significant  panul 
response  for  shaping  n T  tne  excitation  spectrum. 

The  response  frequencies',  obtained  during  these  frequency  sweeps  at  room  temperature 
were  always  different  ftom  those  obtained  in  the  low  level  noise  frequency  sweeps. 

The  low  level  sweeps  were  conducted  in  a  room  with  u  controlled  temperature  and 
humidity  t  n .  ..onm?n',  whereas  the  high  intensity  sweeps  v/ere  mode  in  the  grazing 
incidence  test  facility,  with  a  significant  aii  flow  over  the  panel,  where  no  temperature 
contiol  was  possible.  Aside  from  tne  obvious  temperature  difference,  which  has  beep 
noted  in  past  programs^'  ®,  the  pcnel  response  nonlinearities  served  to  drastically  alter 
the  response  frequency .  Many  of  the  specimens,  particular!)'  those  tested  at  room  temper¬ 
ature,  exhibited  highly  nonlinear  response  wheie  increasing  sound  pressure  levels  coused 
the  response  freejutne  its  io  increase;  only  slight  increases  occurred  in  the  ;t.  ain  amplitudes. 
The  combinations  of  these  two  events  were  in  some  cases  offset:  the  frequency  de  ease 
due  to  a  decrease  in  the  ambient  temperature  was  counteracted  by  an  increase  in  the 
panel  respon»e  freq^ncy  due  to  nonlintority . 
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TABLE  IX 


ALUMINUM  AND  TITANIUM  SPECIMENS 
SUMMARY  OF  SKIN  BUCKLING  TEMPERATURES 
PANEL  CENTER  BAY 


Aluminum  Specimens 

Titanium  Specimens  \ 

Spec imen 

Buckli-j* 

Spec  imen 

Buckling* 

Number 

Temperature  -  °F 

Number 

Temperature  -  °F 

AL-1A 

8 

TI-1A 

38 

-IB 

8 

-2A 

_ ★  ★ 

-2A 

10 

-2B 

* 

-26 

n 

-3A 

57 

-3A 

12 

-3B 

44 

-3B 

8 

-4A 

40 

-4A 

14 

-4B 

59 

-4B 

8 

-5A 

28 

-5A 

34 

-5B 

28 

-5B 

24 

-6A 

_ *  * 

-6A 

7 

-6B 

_ ** 

-6B 

11 

< 

hN 

i 

51 

-7A 

6 

-7B 

50 

-7B 

8 

< 

CO 

1 

12 

-8A 

18 

-8B 

16 

-8B 

15 

-9A 

19 

-9A 

7 

-9B 

17 

-9B 

6 

-10A 

31 

-1 OA 

10 

-10B 

19 

-  !  ► 

10 

-1 1 A 

39 

-11, 

9 

-1 1 B 

49 

•IIP, 

6 

-12A 

6  NOTES:  "Buckling  Temperature  in  °F 

-1  ?B 

5 

above  room  temperature  . 

-13A 

4 

""Room  temperature  specimens. 

-13B 

4 

-1 4  A 

18 

TABLE  X 


ALUMINUM  AND  TITANIUM  SPECIMENS 
SUMMARY  OF  THERMAL  BUCKLING  AMPLITUDES 


Aluminum  Specimens  I  Titanium  Specimens 


Specimen 

No. 

Test 

Temperature 

OF 

Maximum 

Displacement 

In. 

Specimen 

No. 

Test 

Temperature 

°F 

Maximum 

Displacement 

In. 

AL-1A 

300 

0.17 

TI-1A 

400 

- 

-IB 

300 

0.17 

-2A 

RT 

0 

-2A 

300 

0.16 

-2B 

RT 

0 

-2B 

300 

0.20 

-3A 

600 

0.24 

-3A 

RT 

0 

-3B 

600 

0.24 

-3B 

RT 

0 

-4A 

400 

0.17 

-4A 

300 

0.24 

-4B 

400 

0.17 

-4B 

300 

0.24 

-5A 

400 

0.26 

-5A 

300 

0.25 

-5B 

400 

0.23 

-5B 

300 

0.25 

-6A 

RT 

0 

-6A 

RT 

0 

-6B 

RT 

0 

-6B 

RT 

0 

-7A 

400 

0.16 

-7A 

300 

0.35 

-7B 

400 

0.20 

-7B 

300 

0.31 

-8A 

600 

0.27 

-8A 

300 

0.35 

-8B 

600 

0.30 

-88 

300 

0.36 

-9A 

400 

0.12 

-9A 

300 

0.26 

-9B 

400 

0.27 

-9B 

300 

0.25 

-10A 

400 

0.12 

-10A 

RT 

0 

-10B 

400 

0.26 

-1  OB 

RT 

0 

-1 1 A 

400 

0.23 

- 1 1 A 

300 

0.20 

-11B 

400 

0.23 

-1IB 

300 

0.24 

-12A 

RT 

0 

-12B 

RT 

0 

NOTE: 

Maximum  buckling 

1  amplitude 

-13A 

300 

0.26 

measured  at  midpoint  of 

-13B 

300 

0.35 

center  bay . 

-14A 

150 

0.12 

The  excitation  spectrum  for  the  majority  of  the  specimens  had  a  relatively  narrow  bandwidth, 
usually  1 00  Hz,  to  concentrate  the  available  acoustic  energy  of  the  specimen  response 
frequency.  Typical  narrow-band  excitation  spectra  are  shown  in  Figure  35.  Several  of  the 
elevated  temperature  specimens,  especially  the  600°F  titanium  panels,  exhibited  a  relative¬ 
ly  flat  response,  with  no  single  significant  mode.  In  these  coses,  broad-band  excitation 
was  used  to  envelope  the  first  two  oi  three  response  peaks.  The  bandwidth  in  this  case 
was  typically  200  to  300  Hz,  as  shown  in  Figure  36.  Although  a  bend-pass  filter  was 
used  to  cut  out  the  unwanted  low  and  high  frequencies,  the  roll-off  of  the  noise  signul 
in  most  cases  was  very  gradual,  due  to  harmonic  distortion  generated  in  the  electropneu- 
niatic  transducers. 

Probability  density  analyses  were  made  of  the  input  spectra  shown  in  Figures  35  and  36  to 
determine  the  distribution  of  instantaneous  peaks.  These  analyses,  presented  in  Figures 
37  and  38,  closely  approximate  a  Gaussian  distribution  as  noted  on  the  figures. 

e.  Fatigue  Tests  -  The  noise  and  strain  signals,  recorded  on  FM  magnetic  tape, 
were  re-recorded  onto  a  loop  recorder  (10-Second  sample)  for  narrow-band  spectrum 
analyses.  Spectrum  analyses,  with  a  5  Hz  nominal  bandwidth,  of  each  strain  gage 
response  were  made  at  the  beginning  of  the  test  and  periodically  throughout  the  test. 

Sample  strain  response  spectra  are  shown  in  Figures  39  through  44.  The  strain  response 
spectra  provided  the  specimen  response  frequency,  while  the  response  strain  level  was 
taken  as  the  overall  rrm  strain.  Representative  excitation  spectra  were  discussed  in  the 
preceding  section,  and  presented  in  Figures  35  and  36. 

Probability  density  analyses  were  made  of  the  strain  response  spectro  shown  in  Figures 
39  through  44  to  determine  the  distribution  of  response  peaks.  These  analyses  are  shown 
in  Figures  45  through  50  and  approximate  a  Gaussian  distribution. 

Table  XI  summarizes  the  fatigue  test  results  for  the  27  aluminum  specimens  while  Table  XII 
summarizes  the  test  results  for  the  21  titanium  specimens.  These  data  are  plotted  in  the  form 
of  fatigue  curves  in  Figures  51  and  52  for  the  aluminum  and  titanium  specimens,  respectively 
The  data  points  were  plotted  and  a  least  squares  regression  analysis  was  performed  on  the 
dota  for  each  test  temperature.  Since  most  of  the  test  conditions  did  not  contain  sufficient 
data  points  on  which  to  base  a  valid  statistical  analysis,  only  the  centroid  of  the  data  set 
was  used  from  these  analysts.  The  slope  of  the  appropriate  coupon  fatigue  curve  (from 
Figure  16  or  17)  was  then  plotted  through  the  data  centroid  to  establish  the  curves  shown. 
The  applii-.-ble  statistical  properties  ore  itemized  in  Table  XIII  for  reference.  The  fatigue 
curves  of  F  lures  51  and  52  show  the  decrease  in  fatigue  life  due  to  the  combined  effect 
of  thermal  s  ress  and  degradation  of  alloy  properties  caused  by  the  elevated  temperatures. 

The  room  temperature  curves  are  lower  (stresswise  for  equal  life)  than  the  comparable 
coupon  fatigue  curve;  this  is  attributed  to  the  difference  in  response  between  the  two  con¬ 
figurations,  The  stiffened  panel  response  was  very  highly  nonlinear  in  most  cases,  with 
relatively  high  damping,  while  the  coupon  response  was  linear  with  low  damping. 

Fatigue  failures  generally  occurred  in  the  rivet  row  at  the  midpoints  o'  the  long  and  short 
sides  of  the  center  bay.  The  distribution  of  failures  at  these  locations  is  shown  in  Table 
XIV. 
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SPECTRUM  SOUND  PRESSURE  LEVEL  -  dB  SPECTRUM  SOUND  PRESSURE  LEVEL  -  dB 


PROBABILITY  DENSITY  PROBABILITY  DENSITY 


RATIO  OF  INSTANTANEOUS  TO  RMS  AMPLITUDE 


FIGURE  37.  AMPLITUDE  DISTRIBUTION  -  NARROW-BAND 
NOISE  SPECTRUM  OF  FIGURE  35. 


-3-2-1  0  1  2  3 

RATIO  OF  INSTANTANEOUS  TO  RMS  AMPLITUDE 


FIGURE  38.  AMPLT'JDE  DISTRIBUTION  -  BROAD-BAND 
NOISE  SPECTRUM  OF  FIGURE  36. 
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•  :Tt  : 


STRAIN  -  a  IN, /IN  RMS  STRAIN  -  a  IN/IN  R 


SPECTRUM  LEVEL:  119  dB 
OVERALL  STRAIN:  350  a  IN/IN 
RMS_  , 

FILTER  BANDWIDTH:  5  Hz  1 
SCAN  RATE:  1.0  Hz/SEC  I 
SAMPLE  TIME:  10  SEC  f 

TIME  CONSTANT:  3  SEC  - f 


STRAIN  GAGE  LOCATION 


F  RtGUt  Ni_  Y  -  H: 

FIGURE  39.  NARROW-BAND  ANALYSIS 

ALUMINUM  SPECIMEN  AL-1A 


200 

r  -  t Out  NCY  -  n i 


FIGURE  40.  NARROW-BAND  ANALYSIS- 

ALUMINUM  SPECIMEN  AL-3A 


70 


STRAIN 


STRAIN  -  u  IN/IN  RMS 


SPECTRUM  LEVEL:  134  dB 
OVERALL  STRAIN:  350  u  IN/IN  RMS 
FILTER  BANDWIDTH:  5  Hz  _  ' 

SCAN  RATE:  1 .0  Hz/$EC 
SAMPLE  TIME:  10  SEC 
TIME  CONSTANT:  3  SEC 


STRAIN  GAGE  LOCATION 


FREQUENCY  -  H; 

FIGURE  43.  NARROW-BAND  ANALYSIS 
TITANIUM  SPECIMEN  TI-3B 


SPECTRUM  LEVEL:  1 36  dB 
OVERALL  STRAIN:  250  u  IN/IN  RMS 
FILTER  BANDWIDTH:  5  Hz 
SCAN  RATE:  1.0  HZ/SEC 
SAMPLE  TIME:  10  SEC 
TIME  CONSTANT:  3  SEC 


STRAIN  GAGE 
LOCATION 


nn 


FREQUENCY  -  Hz 


FIGURE  44.  NARROW-BAND  ANALYSIS 
TITANIUM  SPECIMEN  TI-10A 


PROBABILITY  DENSITY  PROBABILITY  DENSITY 


FIGURE  46.  AMPLITUDE  DISTRIBUTION  OF  STRAIN  RESPONSE 


SPECIMEN  AL-3A  (REF.  FIGURE  40) 
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RATIO  OF  INSTANTANEOUS  TO  RMS  AMPLITUDE 

FIGURE  47.  AMPLITUDE  DISTRIBUTION  OF  STRAIN  RESPONSE 
SPECIMEN  AL-4A  (REF.  FIGURE  41) 


RATIO  OF  INSTANTANEOUS  TO  RMS  AMPLITUDE 


FIGURE  48.  AMPLITUDE  DISTRIBUTION  OF  STRAIN  RESPONSE 
SPECIMEN  AL  -4B  (REF.  FIGURE  42) 


PATIO  OF  INSTANTANEOUS  TO  RMS  AMPLITUDE 


FIGURE  49  AMPLITUDE  DISTRIBUTION  OF  STRAIN  RESPONSE 
SPECIMEN  TI-3B  (REF.  FIGURE  43) 


-2  -1  0  1 

RATIO  OF  INSTANTANEOUS  TO  RMS  AMPLITUDE 


FIGURE  50.  AMPLITUDE  DISTRIBUTION  OF  STRAIN  RESPONSE 
SPECIMEN  TI-IOA  (REF.  FIGURE  44) 


TABLE  XI 

ALUMINUM  SPECIMENS 
SUMMARY  OF  FATIGUE  TEST  RESULTS 


Spec  imen 

No . 

Test 

Temperature 

°F 

Strain 

Spectrum 

Level 

dB 

Response 

Frequency 

Hr 

Li  fe 

Cyc.  les 

Thermal  | 

/j.in/in  | 

t 

Dy  namic 
u  i  n/i  n 

rms 

AL-1A 

300 

2400 

210 

119 

300 

1  . 08  X  1 06 

-IB 

300 

2400 

•  150 

122 

270 

4.60  x  106 

-2A 

300 

2250 

530 

135 

220 

3.95  x  1C5 

-28 

300 

22.50 

390 

131 

490 

1 .32  x  10* 

-3A 

RT 

0 

270 

134 

255 

3.67  x  KY> 

-3B 

RT 

0 

300 

T  31 

275 

4.46  x  106 

-4A 

300 

260C 

230 

135 

405 

8.46  x  1 

-48 

300 

2600 

220 

1 35 

428 

1.16  x  1  06 

-5A 

300 

2750 

140 

128 

258 

2.55  x  !06  1 

-58 

300 

2750 

90 

126 

280 

3.30x  10‘S  1 

-6A 

RT 

0 

230 

132 

200 

2 . 36  x  1 0)7* 

-6B 

RT 

0 

230 

130 

132 

1 .01  x  KV 

-7A 

3C0 

2600 

250 

128 

235 

8 . 63  x  1  C'5 

-7B 

3C0 

2600 

380 

130 

no 

4.03  «  105 

-8A 

300 

2700 

200 

134 

255 

3.44  x  106 

-SB 

300 

2700 

120 

132 

250 

2.  70  x  1C6 

-94 

300 

2600 

290 

127 

370 

6.66  x  10-' 

-9B 

300 

2600 

200 

127 

345 

9.32  x  105 

- 1 OA 

RT 

0 

390 

137 

205 

3.51  x  106 

-I0B 

RT 

0 

400 

132 

215 

2.40  x  106 

-1 1 A 

300 

2600 

140 

133 

295 

7.12  x  105 

- 11  B 

300 

2600 

140 

131 

280 

9. 27  x  105 

-12A 

RT 

0 

350 

132 

135 

3.53  x  106 

-12B 

RT 

0 

310 

134 

135 

5 . 70  x  1C6 

-13A 

300 

2600 

100 

102 

185 

1  .  7  8  x  1 06 

-13B 

300 

2200 

360 

122 

220 

1  .2!  x  106 

-14A 

150 

750 

370 

123 

325 

9.71  x  105 

*No  skin  fa i  iure. 

|  NOTE:  Fat 

igue  life  (in  cycles)  is  t'ne  product  of  res| 

oonse  frequency  and  test  time  to  failure. 
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TABLE  XI! 

TITANIUM  SPECIMENS 
SUMMARY  OF  FATIGUE  TEST  RESULTS 


Spec inen 

Test 

Strain 

Spectrum 

Response 

Life 

No. 

Temperature 

°F 

Thermal 
p  i  n/  i  n 

Dynamic 

/. 

it  .in 

rms 

Level 

dB 

Frequency 

Hz 

C  yc les 

TI-IA 

400 

1750 

190 

134 

355 

4.48  x  106 

-2A 

RT 

0 

600 

137 

290 

2.14  x  106 

-2B 

RT 

0 

465 

139 

270 

1 . 29  x  1 06 

-3A 

600 

2500 

160 

135 

255 

2.52  x  106 

-3B 

600 

2500 

180 

134 

260 

2.57  x  106 

-4A 

400 

1600 

270 

129 

440 

1 .58  x  106 

-4B 

400 

1600 

230 

128 

460 

2.07  x  106 

-5A 

400 

1000 

200 

126 

390 

3.18  x  106 

-5B 

400 

1000 

260 

113 

495 

6.71  x  106 

-6A 

RT 

0 

180 

133 

200 

1 . 76  x  1 07* 

-6B 

RT 

0 

230 

132 

195 

1 .69  x  107* 

-7A 

400 

2000 

80 

129 

440 

1.11  x  107 

-7B 

400 

1 150 

130 

127 

460 

2.16  x  107 

-8A 

600 

2900 

100 

136 

400 

2.16  x  106 

-8B 

600 

2900 

60 

130 

500 

3.60  x  106 

-9A 

400 

1000 

340 

135 

180 

1 . 94  x  1 06 

-9B 

400 

1550 

340 

134 

180 

4.21  x  106 

-10A 

400 

1600 

190 

136 

225 

7.92  x  106 

-10B 

400 

1600 

200 

140 

200 

1 .01  x  1C7 

-1 1  A 

400 

1680 

260 

137 

230 

1.05  x  I07 

-11 B 

400 

1680 

230 

136 

300 

1.17k  107 

*  No  sk  i  n  fa  i  1  ure . 

NOTE:  Fatigue  life  (in  eye 

les)  is  the  product  of  response  frequency  and  test  time  to  failure. 
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□  ROOM  TEMPERATURE 


CYCLES  TO  FAILURE 


FIGURE  51.  FATIGUE  TEST  DATA 

ALUMINUM  SPECIMENS 


□  ROOM  TEMPERATURE 


A  400  F 


O  60  C  F 


CYCLES  TO  FAILURE 


FIGURE  52.  FATIGUE  TEST  DATA 

TITANIUM  SPECIMENS 


TABLE  XIII 

SUMMARY  OF  STATISTICAL  PROPERTIES 
FOR  STIFFENED  PANEL  FATIGUE  DATA 


A 1  loy 

Test 

Temperature 

°F 

No.  Data 
Points 

Standard 

Deviation 

*c 

Correlation 
Coeffic  ient 

A 1  uminum 

RT 

8 

0.056 

-0.833 

A  luminum 

150 

1 

- 

- 

A 1  uminum 

300 

18 

0,165 

-0.673 

Titanium 

RT 

4 

0.111 

-0.931 

Tilanium 

400 

13 

0.146 

-0.564 

Titanium 

600 

4 

0.209 

-0.622 

TABLE  XIV 

STIFFENED  PANEL  FAILURE  DISTRIBUTION 


Test 

Condition 

A 1  loy 

Stiffener  Rivet  Row 
%  of  Total  Failures 

Frame  Rivet  Row 
%  of  Total  Failures 

Room  Temperature: 

|  Aluminum 

100 

0 

|  Titanium 

100 

0 

Elevated 

j  Aluminum 

21 

79 

Temperature: 

|  Titanium 

77 

23 

80 


Figures  53  and  54  show  some  of  the  failures  obtained  for  eoch  of  these  locations  on  the 
aluminum  and  titanium  specimens.  The  room  temperature  failures  were,  in  all  cases, 
identical  la  those  obtained  on  previous  test  programs**'  ^  in  that  failure  occurred  at  the 
location  of  highest  dynamic  stress  (i.e.  at  the  midpoint  of  the  long  side  of  the  panel  bay). 

It  is  not  clear  what  caused  the  failures  during  the  elevated  temperature  tests  to  occur  at  the 
frame  rivet  iow,  on  the  short  side  of  the  panel,  because  of  the  difference  of  the  failure 
distribution  between  the  aluminum  and  titanium  panels.  Il  is  apparent  from  the  measured 
thermal  strains  that  the  highest  thermal  strains  occur  at  the  short  side  of  the  panel,  whereas 
the  highest  dynamic  strains  occur  at  the  midpoint  of  the  long  side  (for  a  fundamental  mode 
response).  The  combination  of  thermal  and  dynamic  strains  at  the  midpoint  of  the  short 
side  may,  in  some  cases,  prove  more  damaging  than  those  ot  the  other  side.  The  frame 
rivet  row  failures  may  in  some  cases  be  attributed  to  dynamic  response  in  a  higher  order 
mode,  most  probably  the  ( I  ,  2)  mode  .  It  was  not  possible  to  visually  observe  the  specimen 
(esponse  during  the  elevated  temperature  tests  because  of  the  close  piOximity  of  the  lamp 
fixture  to  the  skin  surface. 

Regardless  of  the  failure  location,  the  desired  end  product  (a  fatigue  curve  at  the  various 
rest  temperatures)  is  achievable  by  using  only  the  dynamic  strain  component  at  the  failure 
location.  These  curves  were  discussed  previously  ond  shown  in  Figures  51  ond  52. 

Three  other  types  of  failure  were  observed  during  these  tests:  rivet  failure,  skin  failures 
in  the  outer  bays,  and  substructure  failures.  All  of  these  occurred  predominately  in  the 
titanium  specimens,  with  only  two  rivet  failures,  two  outer  boy  skin  failures,  and  three  sub¬ 
structure  failures  noted  during  the  aluminum  panel  tests. 

The  rivet  failures  generally  were  evidenced  as  fatigue  of  the  rivet  heods  at  the  beginning 
of  the  countersink  toper.  Only  one  rivet  failure,  prior  to  skin  failure,  was  experienced 
during  the  aluminum  panel  test  program.  This  occurred  in  the  0.0o3  skin  AL-8A  specimen 
ot  the  stiffener-to-frame  intersection  .  The  rivets  on  the  stiffener-to-frame  clip  of 
Specimen  AL-5B  failed  at  approximately  the  same  time  that  a  skin  crack  occurred  over 
the  frame  . 

Fifteen  of  the  titanium  panels  experienced  rivet  failures.  The  rivet  failures  occurred  os 
fatigue  of  the  heads,  identical  to  that  described  for  the  aluminum  specimens,  and  as 
deformation  of  the  rivet  heads.  This  deformation,  shown  in  Figure  55,  was  extreme  enough 
in  some  instances  to  allow  the  rivet  head  to  pass  through  the  hole  in  the  skin.  These 
failures  were  usually  randomly  scattered  over  the  surfoce  of  the  panel.  In  a  fpw  instances, 
the  rivet  fa i lures  occurred  in  adjoining  rivets,  allowing  the  skin  to  work  loose  from  the 
substructure;  however,  this  occurred  after  the  initial  skin  failure  in  the  center  bay  on 
a  1 1  spec  i  me  ns  but  T  I  -  1 1  A  and  B.  These  spec  imens,  with  the  thickest  skins  (0.  056), 
proved  stronger  than  the  rivets  and  substructure.  (The  rivets  ond  rivet  spacing  remained 
constant  for  all  specimens).  It  was  concluded  that  the  rivet  size  and  substructure 
thicknes:  should  have  been  increased  in  this  case. 
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PLAN  VIEW  OF  SPECIMEN 


SPECIMEN  AI-12A 


FIGURE  53.  TYPICAL  SKIN  FAILURES 
ALUMINUM  SPECIMENS 


FIGURE  55.  TYPICAL  RIVET  FAILURE 
TITANIUM  SPECIMENS 


The  only  aluminum  specimens  to  incur  cracks  in  the  outer  boy  skins  were  AL-5A  and  B. 
These  failures  were  judged  to  be  small  enough  and  sufficiently  far  removed  from  the 
initial  outer  bay  crock  so  as  to  have  negligible  effect  on  the  center  boy  response  and 
thermal  strain.  Approximately  50%  of  the  titanium  specimens  experienced  skin  cracks 
in  the  outei  bay  skins  prior  to  or  concurrent  with  the  center  bay  skin  cracks.  In  all 
cases,  these  were  either  sufficiently  far  removed  from  the  center  bay  (i.e.,  at  the  end 
rivet  on  o  frame)  or  small  enough  that  their  presence  did  not  significantly  affect  the 
center  bay  response. 

f.  Effect  of  Sealant  on  Fatigue  -  As  mentioned  previously,  eight  of  the  aluminum 

specimens,  of  two  different  designs,  were  fabricated  using  faying  surface  sealant  at  all 
joints  where  metal-to-metal  contact  occurred.  These  specimens  were  identical  in  all 
respects,  except  for  the  sealant,  to  eight  of  the  other  specimens,  as  identified  in 
Table  XV. 


TABLE  XV 

SPECIMEN  DESIGNATION  FOR  SEALANT  EVALUATION 


Spec imen 

Size 

Test 

Temperature 

Specimen  Idem 

i  fication 

Without 

Sealant 

With 

Sealant 

6  x  12  x  0.032 

RT 

AL-3 

AL-10 

6  x  12  x  0.032 

300°  F 

AL-4 

AL-11 

9  x  18  x  0.040 

RT 

AL-6 

AL-12 

9  x  18  x  0.040 

300PF 

Al-7 

AL-13 

The  measured  dynamic  data  for  these  specimens  formed  the  basis  for  comparison 
of  the  effects  of  sealant  on  response  and  life.  Room  temperature  fundamental  mode  fre¬ 
quencies  and  damping  rctios  were  tabulated  in  Table  VI!  for  these  specimens. 

A  comparison  of  both  of  these  parameters,  individually,  shows  that  any  effect  that  the 
sealant  has  is  lost  in  the  scatter  in  the  data.  Figure  56,  which  shows  the  damping 
ratios  plotted  versus  the  fundamental  frequency,  clearly  shows  that  faying  surface  sealant 
produces  no  discernable  difference  on  the  dynamic  response.  A  further  evaluation  of  the 
fatigue  data  of  Table  XI  and  Figure  51  reinforced  this  conclusion.  The  results  of  both 
designs  are  intermingled  in  the  scatter  band  of  the  data. 

Since  faying  surface  sealant  has  no  net  effect  on  dynomic  response  or  fatigue  life  for 
the  type  of  structures  considered,  it  can  be  concluded  that  the  design  criteria  generated 
on  previous  programs,  without  sealant,  are  volid  for  structures  where  sealant  is  used.  (This 
includes  most  of  the  secondary  structure  on  current  generation  aircraft,  ) 
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IV  -  CORRELATION  OF  ANALYTICAL  AND  EXPERIMENTAL  RESULTS 


The  principal  objective  of  this  program  was  the  establishment  of  empirical  design  criteria 
for  aircraft  structures  subjected  to  combined  acoustic,  and  thermal  environments.  The 
analytical  development  provided  the  interrelation  omong  the  applicable  parameters,  while 
the  experimental  program  provided  measured  values  of  the  individual  parameters .  Corre¬ 
lation  of  these  data  then  can  provide  the  necessary  empirical  constants,  cr  factors,  for 
derivation  of  the  design  methods. 

The  derivation  of  empirical  parameters  follows  the  sequence  in  which  each  parameter  is 
required  in  the  thermal  and  dynamic  stress  relations,  respective ly  .  Generally,  the  meas¬ 
ured  data  were  plotted  versus  calculated  values  (using  the  resuits  of  the  simple  panel 
analysis),  and  a  least  squares  linear  regression  analysis  was  performed  to  determine  empir¬ 
ical  constants.  The  least  squares  regression  analysis  follows  the  method  of  Reference  4, 
Appendix  II  .  In  some  instances  an  exponential  regression  analysis  was  performed  to  deter¬ 
mine  whether  exponentia  I  factors  should  be  considered;  however,  in  eoch  cose  the  linear 
relationship  was  judged  to  provide  the  best  correlation.  Most  ot  the  correlation  plots 
exhibited  an  off-set  '  i  "he  y-intercept  (measured  data  axis),  because  of  non-linearities 
and  other  variations  n  hie  r  easurea  data.  In  all  instances  this  off-set  was  juaged  to  be 
small,  in  relation  to  h  -  vor:  :nce  in  the  data,  so  that  the  final  equation  could  be  estab¬ 
lished  by  the  apprcxi  mu  r-  .ation 

(  )  -  m(  ) 

e  c 


where  () ^  -  empirical  value 

()  ^calculated  value,  from  analytical  results 

m  =  slope  of  approximate  regression  line  through  the  origin  and  the  data  centroid  (x,  y) . 

The  nine-bay  analysis  results  were  also  correlated  with  the  test  data  to  establish  empirical 
relations  for  use  in  a  digital  computer  program.  This  correlation  effort  is  described  in 
Appendix  IV  . 


A,  Thermal  Stress 


As  the  thermal  stress  is  dependent  on  the  skin  buckling  remperoture  and  buckling  amplitude, 
these  parameters  must  be  established  first. 

1  .  Skin  Buck  I  i  ng  T  emperature 

Skin  critical  buckling  temperatures  were  listed  in  Table  X  and  were  found  to  be  consis¬ 
tently  higher  than  the  values  calculated  from  the  analysis  by  Equation  (5o).  The  correlation 
plot  of  the  buckling  temperature,  Figure  57,  shows  that  considerable  scatter  was  present  in 
the  measured  data,  as  evidenced  by  a  correlation  coefficient,  Rc,  of  0.718.  (Note  that 
the  correlation  coefficient  between  two  variables  is  zero  when  no  correlation  exists,  and 
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FIGURE  57.  BUCKLING  TEMPERATURE  CORRELATION 


unity  when  their  re  lationship  can  be  represented  exactly  by  a  straight  line  .)  Repeated 
measurements  on  the  same  specimen  revealed  that  the  critical  buckling  temperature  could 
vary  sianif  icantly ,  depending  on  the  heating  rate  and  the  temperature  differential  between 
the  skm  and  substructure.  Pre-stresses  introduced  during  specimen  fabrication  (i  .e . ,  mis¬ 
alignment  of  frames  and  stiffeners,  and  other  factors)  olso  influence  the  buckling  temperature. 


The  equation  of  the  regression  line,  (A),  through  the  data  points  of  Figure  57  has  a  2.59°F 
off-set  on  the  measured  buckling  temperature  axis.  A  line,  (B),  plotted  through  the  origin 
and  the  centroid  (x,  y)  of  the  data  set  is  also  shown  on  the  figure .  Comparison  of  the  two 
equations  shows  very  little  difference  in  the  results,  particular!  ly  in  view  ot  the  scatter 
present  in  the  data.  Hence,  the  simpler  approxi motion.  Equation  (B),  was  selected  to 
represent  the  corn  lotion  between  measured  and  analytical  volues,  which  results  in  the 
following  empirical  relation  for  multi-bay  structures: 


-  5.25  h2  F , , 
a  ab(  1  +  v) 


(60) 


where  T  is  the  temperature  increment  above  ambient.  The  temperature  ratio,  r,  remains 
ce 

identical  to  the  analytical  definition,  or 


r  -  T/T 

c 


e 


It  should  again  be  pointed  out  that  the  skin  temperature,  T,  is  the  temperature  increase 
above  ambient . 

2.  Buckling  Amplitude 

The  panel  skin  buckling  amplitude,  W  ,  is  related  to  the  critical  buckling  temperature 
by  Equation  (22).  For  a  specific  configuration,  the  buckling  amplitude  varies  directly  with 
the  square  root  of  the  temperature  increase,  T.  Therefore,  the  measured  temperature  in¬ 
crease,  corresponding  to  a  measured  buckling  amplitude,  was  used  to  compute  analytical 
duckling  amplitudes  by  Equation  (22). 

Because  of  the  multitude  of  data  points  (temperature-displacement  data  at  approximately 
25°F  increments  for  each  elevated  temperature  specimen  to  the  maximum  test  temperature), 
a  digital  computer  program  was  written  to  accomplish  the  individual  computations .  The 
empirical  buckling  temperature  relation,  (60),  was  used  instead  of  the  measured  buckling 
temperatures  to  minimize  the  effects  of  cumulative  errors.  The  specimen  dimensions  and 
temperature-displacement  data  were  entered  into  the  program  individually  for  each  speci¬ 
men  and  tnen  grouped  into  equal  aspect  ratios.  The  data  were  then  plotted  on  a  digital 
plotter  in  one  of  the  following  forms: 

o  Measured  versus  calculated  displacement  ratio 

o  Displacement  ratio  versus  temperature  rotio. 
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The  first-  of  these  plot  formats  is  illustrated  in  Figures  58  through  60  for  aspect  ratios  of 
I  .5,  2  and  3,  respectively.  As  on  the  buckling  temperature  correlation  plot,  two  lines 
are  shown  on  each  figure.  The  first,  line  (A),  represents  the  least  squares  regression  line 
while  line  (8)  represents  the  zero-origin  approximation.  Because  of  rhe  slight  difference, 
the  approximate  regression  line  was  used  in  all  cases.  The  slopes  of  these  lines  are  tabu¬ 
lated  in  Table  XVI  to  illue*rate  the  variation  in  slope  with  aspect  ratio. 

TABLE  XVI 


SLOPE  OF  BUCKLING  AMPLITUDE  CORRELATION  PLOTS 


!  ASPECT  RATIO 

SLOPE 

b/  a 

F.l 

1  .5 

2.17 

2.167 

2.0 

2.50 

2.571 

3.0 

3.33 

3.046 

A  simple  average  of  these  slopes,  when  used  to  derive  an  empirical  relation,  yields 
empirical  displacements  higher  than  measured  for  aspect  ratios  less  than  2,  while  the 
reverse  is  true  for  aspect  ratios  greater  than  2.  The  effect  of  aspect  rotio  on  the  slope 
is  shown  in  Figure  61 .  Variation  of  the  power  of  F]  j  resulted  in  the  conclusion  thof  an 
additional  power  of  0.75  provided  the  minimum  error  in  the  slopes  for  each  of  the  aspect 
ratios.  This  gives  an  empirical  buckling  amplitude  relation  of 


2.50hF 


1.75 

1 1 


f  ( r  —  1 )/ R 


1/2 


(61) 


where  R  remains  as  defined  in  Equation  (19a).  This  equation  shows  a  greater  dependency 
on  aspect  ratio  than  that  indicated  by  the  analytical  equation.  This  additional  dependence 
on  aspect  ratio  results  in  o  reversal  of  the  trends  given  by  the  analysis  and  shown  in  Figure  3. 
The  effect  of  this  aspect  ratio  correction  will  be  mote  clearly  defined  in  the  later  Design 
Method  Section. 

The  measured  displacement  data  were  plotted  versus  temperature  rotio  as  shown  in  Figures  62 
through  64  for  the  three  aspect  ratios.  The  empirical  equation  is  also  shown  for  each  aspect 
ratio,  and  illustrates  the  agreement  between  the  empirical  trend  and  the  dato  with  varying 
aspect  ratio.  Since  most  of  the  test  specimens  had  on  aspect  ratio  of  2,  the  corresponding 
displacement  plot  of  Figure  63  shows  the  magnitude  of  scatter  that  can  be  expected  in  any 
given  set  of  buckling  amplitude  measurements.  To  quantify  the  variance,  all  data  were 
plotted  in  the  correlation  plot  of  Figure  65,  using  Equation  (61).  The  correlation  coeffi¬ 
cient  and  standard  deviation  from  the  mean  are  indicated  in  the  figure,  as  are  the  limits  of 
the  60%  confidence  bands.  Hence,  for  any  given  displacement  estimate,  the  maximum 
error  expected,  for  a  confidence  level  of  60%,  should  be  less  than  *0.59  times  the  skin 
thi  ckness . 
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MEASURED  DISPLACEMENT  RATIO 
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f  FIGURE  62.  COMPARISON  OF  MEASURED  AND  EMPIRICAL  BUCKLING  AMPLITUDE 
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FIGURE  63.  COMPARISON  OF  MEASURED  AND  EMPIRICAL  BUCKLING  AMPLITUDE 
ASPECT  RATIO  =  2 


FIGURE  64,  COMPARISON  OF  MEASURED  AND  EMPIRICAL  BUCKLING  AMPLITUDE 
ASPECT  RATIO  =  3 
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(W  /h)  -  EMPIRICAL  DISPLACEMENT  RATIO 
o'  e 


FIGURE  65.  BUCKLING  AMPLITUDE  CORRE LATION 
ALL  ASPECT  RATIOS 


3. 


Thermal  Strain 


Thermal  strairts  were  measured  during  the  test  program  on  each  of  the  elevated  temperature 
specimens  at  increasing  temperature  increments  of  approximate!/  25°F.  The  buckling 
amplitude  computer  program  was  adapted  to  accept  the  thermal  strain  input  and  compute 
analytical  strains,  as  well  as  to  plot  the  data.  The  analytical  stresses  at  the  midpoints  of 
the  sides  were  converted  to  strains  for  direct  comparison  with  the  measured  data  and  to 
simplify  the  computation,  since  the  elastic  modulus  is  temperature-dependent. 


The  analytical  buckled  panel  strain  is  thus  gisen  by: 


e  (x,b/2)  =  - 

r  aT 

+  rt2^2 

b_  (2-v^)  +  vo 

i>vj 

8ab(  1  -V*) 

a  b  _ 

(62  a) 


r  i 

|2l 

1*  Avo2 

_o_(2-v2)  +  .,  b 

Li-vJ 

8ab(l-v^) 

.  b  o  . 

(62b) 


The  first  component  of  these  stresses  is  the  expansion  strain  aT/(l-v)  ,  which  is  linear 
and  applies  for  all  values  of  r.  It  was,  therefore,  advantageous  to  delete  this  strain  com¬ 
ponent  from  both  the  analytical  value  and  the  measured  strains.  This  was  accomplished 
by  reducing  the  magnitude  of  the  measured  strains  by  the  value  of  the  expansion  strain  at 
the  temperature  for  which  the  data  were  taken.  The  remaining  strain  components,  due  only 
to  the  panel  buckling,  were  correlated  directly  against  one  of  the  following  calculated 
strains  due  to  buckling: 


X 

c 


-Av2 

o 

8ab(l  -v  2) 


b  (2-v  )  + 
a 


(63a) 


-_1 i  r.e^.vb 

8ab(!-v2)  'Lb  a 


(63b) 


The  empirical  equation  for  buckling  amplitude,  (6l),  wosused  in  these  computations  to 
reduce  the  cumulative  error. 


The  test  data  were  measured  at  three  locations,  in  the  direction  indicated  below: 

1  .  Center  of  panel  (x  =  a/2,  y  -  b,  2) 

x 

2 .  Midpoi  nt  of  pane  I  long  side  (x  -  0  or  a ,  y  =  b/2)  ~  e 

x 

3.  Midpoint  of  panel  short  side  (x  =  a-  2,  y  =  0  or  b)  ~  c 

y 

Measured  strains  at  the  three  locations  were  then  plotted  versus  calculated  strains,  with  the 
result  shown  in  Figures  66  and  67.  Absolute  values  of  both  measured  and  calculated  dato 
were  used,  even  though  the  sign  of  the  buckling  stress  is  positive.  The  negative  measured 
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measured  strain 


c 


-  ANALYTICAL  STRAIN  -  a  IN/IN 


FIGURE  67.  THERMAL  BUCKLING  STRAIN  CORRELATION 

LOCATION  3  -  PANEL  SHORT  LENGTH  MIDPOINT 


strains  on  the  plots  are  the  result  of  subtraction  of  the  expansion  strain.  The  least  squares 
regression  lines  and  the  simpler  zero-origin  approximation  are  again  indicated,  with  the 
difference  between  the  two  neg  I  ig  ible  in  comparison  to  the  scatter  in  the  data  .  Therefore  , 
the  simpler  approximation  was  selected  in  each  case,  resulting  in  the  following  empirical 
relations  for  thermal  stress: 


o  Midpoint  of  panel  long  side 

c  =  -[Eo!]  +0.82EWo2  |b(2-v2)  +  va 


1 


6  L’-vJ  ab(l-,2)  La 


bj 


o  Midpoint  of  panel  short  side 

c  *  -  j  EoT  !  +  1.66  EW^  [  a  (2-v2)  +  vb  ] 

!  b  a  j 


e  L 1  J  ab(l-v2)  L1 


(64a) 


(64d) 


These  equations  are  applicable  only  for  the  rnid-plane  of  the  panel  in  the  direction  indicated 
Strains  for  other  locations  can  be  readily  derived  from  these  relations  by  use  of  the  analytica 
equations  involving  spatial  location. 


The  results  of  the  strain  corre lation  are  illustrated  in  Figures  68  through  70,  which  show  the 
measured  data  and  the  empirical  equation  plotted  versus  temperature  ratio  for  several  of  the 
test  specimens.  The  agreement  in  the  trends  of  the  empirical  relation  and  the  measured  data 
is  considered  good  in  view  of  the  considerable  scatter  which  exist  in  the  data  .  . 


The  influence  of  temperature  on  the  coefficient  of  thermal  expansion,  a ,  is  reflected  in  the 
empirical  curves  of  Figures  68  through  70  by  the  curvature  in  the  empirical  line.  The  value 
of  a  for  aluminum  increases  (with  increasing  temperature)  for  ail  temperatures  below  300°F, 
wh  i  le  o  for  titanium  reaches  an  upper  limit  at  a  temperature  of  400 °F . 


B.  Dynamic  Stress 

Computation  of  the  dynamic  stress  involves  the  frequency  and  sound  pressure  spectrum  level 
in  addition  to  certain  specimen  dimensions.  The  frequency  response  at  a  specific  tempera¬ 
ture  above  ambient  is,  in  turn,  dependent  on  the  room  temperature  fund  imental  frequency. 
Hence,  the  derivation  of  the  empirical  dynamic  parameters  will  be  disc  ^ised  in  the  order 
in  which  they  will  be  used. 

1,  Ambient  Temperature  Fundamental  Frequency 

Measured  room  temperature  fundamental  frequencies,  from  Tables  VII  and  VIII,  were  plotted 
versus  values  calculated  by  Equation  (30).  The  results  are  shown  in  Figure  71  for  all  alumi¬ 
num  and  titanium  specimens,  with  the  two  forms  of  regression  iine  indicated  .  Since  there 
was  very  little  difference  in  the  two  slopes,  the  approximation  was  selected,  and  the 
following  empirical  relation  resulted  for  room  temperature  fundamental  mode  frequency: 
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r  -  TEMPERATURE  RATIO 


- EMPIRICAL  EQUATION  (640) 

0  LOCAilON  1  -  CENTER  cr  PANEL 
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COMPARISON  Or-  EMPIRICAL  AND  MEASURED  THERMAL  STRAiNS 


THERMAL  STRAIN  -  M  IN/IN  -:x  ~  THERMAL  STRAIN  -  H  IN/lN 


-  EMPIRICAL  EQUATION  (64o) 

O  LOCATION  I  -  CENTER  OF  PANEL 


6  LOCATION  2  -  MIDPOINT  OF  LONG  SIDE 


THERMAL  STRAIN  -  /i IN/IN  O  "  THERMAL  STRAIN  -  f»IN/lN 


-  EMPIRICAL  EQUATION  (64b) 

A  LOCATION  3  -  MIDPOINT  OF  SHORT  SIDE 


FIGURE  70,  COMPARISON  OF  EMPIRICAL  AND  MEASURED  THERMAL  STRAINS 
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MEASURED  FREQUENCE  -  H 


FIGURE  71 .  FUNDAMENTAL  MODE  FREQUENCY  CORRELATION 
AMBIENT  TEMPERATURE 
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or,  alternatively, 


0.79 


luM- 

ab  [  v 


(1-v2) 


1 1/2 


(65b) 


This  equation  is  based  on  a  simply  supported  plate  analytical  model.  Use  of  more  realistic 
boundaries  for  the  analytical  model  is  not  justified  because,  like  the  buckling  temperature, 
a  certain  amount  of  scatter  is  present  in  the  measured  data  for  a  particular  specimen  con¬ 
figuration.  This  scatter  is  attributable  to  the  same  factors  that  create  the  scatter  in  the 
buckling  temperature,  as  discussed  previously. 


The  correlation  coefficient  for  the  plot  of  Figure  71  is  Rc  -  0.915,  hence  the  probability 
of  accurately  estimating  the  fundamental  mode  frequency  for  a  specific  design  configura¬ 
tion  is  considered  good.  For  instance,  the  error  in  the  estimated  frequency  should  be  It- ss 
than  i  18  Hz  for  a  confidence  level  of  60%. 


2.  Elevated  Temperature  Fundamental  Frequency 


The  fundamental  frequency  response  of  a  heated  panel  is  given  by  Equation  (30)  as  c 
function  of  temperature  ratio  and  the  room  temperature  frequency.  At  the  critical  buckling 
temperature,  r  =  1  ,  the  analytical  frequency  decreases  to  zero,  as  shown  in  Figure  4.  How¬ 
ever,  as  shown  in  Figure  33,  the  measured  data  did  not  exhibit  this  tendency.  Measured 
frequency  ratios,  f(r)/fQ,  were  plotted  against  the  measured  temperature  ratio,  r,  with  the 
data  from  all  specimens  plotted  on  the  same  graph.  The  data  distribution  is  indicated  by 
the  shaded  area  of  Figure  72;  the  data  points  were  distributed  uniformly  within  the  limits 
of  the  scatter  band.  The  data  were  separated  into  two  groups,  above  and  below  the  criti¬ 
cal  temperature,  and  me'-sared  frequency  ratios  plotted  against  calculated  frequency  ratios. 
A  least  squares  regression  line  was  computed  for  each  group  and  then  slightly  modified  to 
provide  identical  values  at  r  -  1  so  that  a  continuous  curve  would  result.  The  empirical 
equations  are  thus 


f(r)  =f 
e  o 


f(r)  =f  t 
e  o  u 


0.60  +  0.40  (l-r)’/2  1 

1  2  1 

0.60  +  0.44  (r-1)  '  ! 


(0  <  r  <  1) 


(r  S  1) 


(66) 


The  curve  produced  by  this  equation  is  also  shown  in  Figure  72.  No  means  were  available 
to  reduce  the  scatter  in  the  frequency  ratio  data  since,  for  identical  configurations,  the 
frequency  ratios  varied  from  the  minimum  to  the  maximum  shown  in  Figure  72.  Much  of 
this  scatter  is  directly  related  to  errors  accumulated  in  the  room  temperature  fundamental 
frequency  and  the  buckling  temperature  measurements . 
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3.  Dynamic  Strain 


Dynamic  stresses  were  calculated  using  the  following  equations  from  Reference  4  for 
clamped  edge  conditions: 


o  Rivet  row  at  midpoint  of  long  side  (x  -  0,  y  -  b/2) 

,/2 


~J7i 


o  Rivet  row  at  midpoint  of  short  side  (x  -  a/2,  y  -  0) 

M/2 

o  /  a  \  -nu  i  hi/  f 
Cy  =  ~J7i- 

The  aspect  ratio  parameter,  AR,  is  defined  as 


>\2  no 

f(0 

;/  ar 

_ 

side  (x 

-  a/1 

i  \2  no 

f(r) 

J  AR 

.  . 

ined  as 

I2  *  3  (a/b)2 

(67a) 


(67b) 


Dynamic  strain  analyses  were  made  several  times  during  each  test,  in  the  form  of  a  narrow- 
bond  analysis ,  as  discussed  previously.  The  re  ponse  frequencies,  sound  pressure  levels,  and 
overall  strain  levels  from  these  analyses  wer-  used  in  the  dynamic  strain  correlation  for  the 
two  locations  given  above.  These  locations  correspond  to  locations  2  and  3  for  the  thermal 
strain  correlation. 


The  calculated  strains  for  each  measurement  location  were  plotted  versus  the  corresponding 
measured  strains  as  shown  in  Figure  73,  and  the  least  squares  regression  lines  were  computed. 
As  noted  in  the  figures,  the  least  squares  curve  fit  has  a  very  gradual  slope,  with  residual 
values  of  176  and  137  ^in  /in  at  zero  calculated  strain.  The  approximate  zero-origin 
curve  fit  is  also  shown  for  each  case.  From  a  practical  viewpoint,  the  strain  without  noise 
excitation  should  be  zero,  which  is  not  the  case  for  the  least  squares  curve  fit.  Since  a 
large  variance  is  present  in  the  data  (correlation  coefficients  of  less  than  0.30  for  both 
plots)  either  line  can  be  used  without  loss  of  accuracy.  Because  of  this  lack  of  correlation, 
a  conservatism  factor  of  2  was  applied  to  the  approximate  relations  to  produce  the  following 
empirical  dynamic  stress  relations: 


o  Rivet  row  at  midpoint  of  long  side  (x  -0,  y  -  b/2) 

1 1/2 


;b\2  no 

f(r) 

\fj  AR 

i 

o  Rivet  row  at  midpoint  of  short  side  (x  -  a/2,  y  =  0) 

1 1/2 

=  1.30  {£1 

y. 


/cm2  no 

f(r) 

\h  /  AR 

' 

(68a) 


(68b) 


The  probability  of  accurately  predicting  dynamic  stresses  for  a  particular  application  is  low 
because  of  the  scatter  in  the  data;  hence,  these  relations  should  only  be  used  to  provide 
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FIGURE  73.  DYNAMIC  STRAIN  CORRELATION  -  ALL  TEMPERATURES 


= JL 


gross  estimates  of  the  stress  magnitude.  Assuming  a  60%  confidence  level,  errors  in  the 
estimated  stresses  can  be  expected  to  be  on  the  order  of  +95  and  +70  ^ir\/in  (rms)  from  the 
mean  value  for  locations  2  and  3,  respectively. 

The  above  results  have  been  derived  using  measured  strain  data,  regardless  of  the  test  tem¬ 
perature.  This  was  done  because  separation  of  the  data  by  test  temperature  produced  no 
significant  difference  in  empirical  stress.  Figure  74  shows  the  room  and  elevated  tempera¬ 
ture  data  separately  far  strain  location  2,  and  shows  approximately  4%  difference  in  the 
slopes  of  the  approximate  regression  line.  Such  a  small  difference,  in  comparison  with 
the  data  variation,  does  not  justify  establishment  of  separate  relations  for  room  and 
elevated  temperature  dynamic  strains. 
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V  -  DESIGN  METHODS 


A  useful  tool  for  the  design  engineer  is  the  design  nomogroph,  which  graphically  displays 
an  equation  for  rapid  solution.  The  empirical  equations  of  the  preceding  section  were 
formulated  into  such  nomographs  and  are  presented  in  the  following  subsections.  The 
results  of  the  empirical  derivation,  together  with  existing  room  temperature  criteria,  are 
summarized  here  to  clarify  application  of  the  design  technique.  The  empirical  subscript, 
e,  has  been  dropped  from  all  equations  presented  here  to  simplify  the  results. 

A.  Ambient  Temperature  Design  Cri terio 

The  design  criteria  ar  ambient  temperatures  are  unchanged  from  existing  criteria.  Only 
the  dynamic  response  of  the  structure  is  involved  in  the  design  os  long  os  the  ombient 
temperature  state  does  not  cause  buckling  of  the  skin. 


I  .  Skin  Design 

The  skin  design  criteria  of  AFFDL-TR -67-156^  are  valid  for  aluminum  structures  at  ambient 
temperatures.  Figure  75  depicts  the  nomenclature  for  i  simple  flat  panel  which  is  repre¬ 
sentative  of  a  single  bay  of  a  stiffened-skin  structure .  The  dynamic  stress  at  the  midpoint 
of  the  long  side  is  given  by  Equation  (36c),  Reference  4,  as 


a  =  I  .62  x  10 
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» j 


1/4  o’-25  *(f)  (b/a)'‘75~ks; 

h1.75  ..56(ar).84 
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Dynamic  stresses  at  this  po;nt  are  also  given  by  Equation  (68a),  or 
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where  fQ  is  given  by  Equation  (65),  or 


fo=0-79Fn  h 

ab 
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The  aspect  ratio  parameter  is  defined  as 


AR  =  3(b/a)2  +  3(a/b)^  +  2 


while  the  pressure  density  is  defined  by 

$(f)  =2.9]  x  lO^L  /20)  ~9^ 


ps. 
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Comparison  of  these  stress  equations  shows  that ,  for  identical  configurations,  Equation 
(69a)  gives  higher  stresses  than  (69b).  Since  both  ore  based  on  the  same  analytical  model, 
the  difference  is  in  the  data  on  which  these  empirical  relations  are  based.  These  empiri¬ 
cal  relations  can  then  be  considered  as  bounds  for  predicting  dynamic  stresses.  The  form 
of  the  latter  equation  lends  itself  to  much  easier  solution. 

Figure  76  presents  a  nomograph,  based  on  Equation  (69a)  from  Reference  4(  for  stiffened- 
skin  plating  design  .  This  nomograph  was  adapted  for  titanium  structures  by  using  the  dato 
of  Section  III  . 

EXAMPLE:  A  flat,  aluminum  alloy,  stiffened  structure  is  to  be  designed  to  withstand  on 
estimated  spectrum  level  of  120  dB  for  5x10®  cycles.  The  skin  design  is  determined  by 
the  fo  I  low  ing : 

Assume:  o  Damping  ratio:  '  -  0.012 

o  Stiffener  spacing:  a  -  4.75  inches 
o  Aspect  ratio;  b/a  u  1  .5 

Enter  the  nomograph,  Figure  76,  with  the  design  life  and  follow  through  the  nomograph, 
as  indicated  by  the  arrows,  to  obtain  o  skin  thickness  h  -  0.032  inch. 

The  fundamental  frequency  is  calculated  by  Equation  (65)  as  f  =  340  Hz.  At  this  fre¬ 
quency,  the  service  environment  spectrum  level  is  checked  with  the  spectrum  level  used 
above.  If  necessary,  an  iteration  can  be  made  to  obtain  agreement. 


2.  Stiffener  Flange  Design 

The  acoustic  loading  on  the  surface  of  a  stiffened  panel  is  transferred  to  the  substructure 
predominately  by  a  transverse  shear  loading,  causing  the  open  section  stiffeners  to  bend 
and  twist.  The  stiffener  loading  is  reacted  along  the  skin-stiffener  attachment  (rivet)  line 
and  at  the  clip  attachments  to  the  frames.  The  resulting  stresses  at  the  stiffener  flonge  ore 
given  by  Equation  (67),  Reference  8,  as 
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(70) 


This  relafion  is  valid  only  for  the  fundamental  mode  of  the  panel.  The  above  flange  stress 
is  used  in  conjunction  with  a  fatigue  curve  developed  for  flange  failures;  this  curve  is 
presented  in  Figure  77  (Reference  8,  Figure  44). 

EXAMPLE:  A  flat,  aluminum  o.-'oy,  stiffened  structure  is  to  be  designed  to  withstand  on 
estimated  spectrum  level  of  120  dB  for  5  x  10®  cycles.  The  stiffener  design  is  determined 
by  the  foil  owing: 
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NOTE:  ALUM  I  MUM  ALLOY 

DESIGN  NOMOGRAPH 
FROM  REFERENCE  £  . 


STIFFENED  PANEL  SKIN  DESIGN  NOMOGRAPH 
AMBIENT  TEMPERATURE 


FIGURE  77.  STIFFENED  PANEL  STIFFENER  FATIGUE  CURVE 
AMBIENT  TEMPERATURE  (FROM  REFERENCE  8) 


Assume:  o 


Damping  ratio:  £  =  0.012 
o  Stiffener  spacing :  a  -4.75  inches 
o  Aspect  ratio:  h/a  =  1.5  (F^  ^  —  2.17) 


From  the  previous  example, 


h  =  0 .032  inch 
f  =  340  Hz 


The  stiffener  is  a  zee. section  0.040  inch  thick  with  a  flange  width  of  0.75  inch  and  height 
of  1 .25  inch;  the  section  properties  give  a  volue  of  I*  -  0.01255  inch  . 

From  Equation  (70),  the  attachment  stress  is 


cf  =  0.9 


0.0121  (7. 1 25) 3  (1 .25)  (2.9  x  10"3)/  340  \  X'2 


(0.01255)  (2.17)' 


0.012 


1/5 


-  2 .25  ksi 


rms 


From  Figure  77,  the  life  is  estimated  to  be  N  -  9  x  10  cycles,  or  considerably  less  than 
the  design  requirement. 


The  above  procedure  should  then  be  repeated  using  a  thicker  or  deeper  zee  stiffener  until 
the  desired  life  is  achieved.  It  is  also  possible  to  reduce  the  stiffener  spacing,  thereby 
reducing  the  fundamental  frequency  and  perhaps  altering  the  excitation  (dependent  on  the 
spectrum  shape) . 


B,  Elevated  Temperature  Design  Criteria 

The  elevated  temperature  design  criteria  are  used  in  essentially  the  same  sequence  in 
which  they  were  discussed  in  Section  IV.  Since  more  fhon  one  method  of  application  is 
available,  all  the  criteria  will  be  summarized  and  then  followed  by  examples  of  usage. 


1.  Skin  Buckling  Temperature 

Tho  critical  buckling  temperature  of  a  single  panel,  such  as  that  shown  in  P:gure  75,  is 
given  by  Equation  (60),  or 

2 

T  3  5.25  h  F]|  ~  QF  obove  omblent 
a  ob(l  +  v) 
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The  temperature  ratio  is  then  defined  as 


r  =  T/Tc 

where  T  is  the  temperature  rise  of  the  structure  above  ambient.  Figure  78  is  a  nomograph 
of  the  above  equation  for  a  constant  value  of  Poisson's  ratio.  The  volue  of  v  =  0.32  was 
selected  os  representative  of  the  most  commonly  used  aircraft  alloys. 


2.  Skin  Buckling  Amplitude 

The  empirical  skin  buckling  amplitude  is  given  by  Equation  (61),  or 

~  inches 
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where  R  is  defined  by  Equation  (19a)  as 

2s,  2 


R  =  3  j  (5-v2)F]  j  -2(5  +  v)(1-v) 


A  nomograph  to  predict  buckling  amplitudes  is  shown  in  Figure  79;  this  nomograph  was 
also  developed  for  a  constant  value  of  v  =  0.32. 


3.  Thermal  Stress 


Thermal  stresses  due  to  in-plane  expansion  and  skin  buckling  are  given  by  Equations  (64), 
for  the  midpoint  of  each  side,  or: 


o  Midpoint  of  panel  long  side 


°x  = 

-  EoT  +0.82  EW02 
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o  Midpoint  of  panel  short  side 
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Thermal  stresses  must  be  computed  at  the  midpoint  of  both  sides,  since  the  short  side 
stress  is  greater  than  that  at  the  center  of  the  long  side.  This  is  opposite  to  the  magi- 
tudes  of  the  dynamic  stresses  at  the  two  locations. 

The  above  relations  were  simplified  to 
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PANEL  SKIN  TEMPERATURE  -  F  ABOVE  AMBIENT 
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for  development  of  design  nomographs.  These  respective  stresses  are  the  in-plane  expan¬ 
sion  and  buckling  stresses  as  inferred  from  ihe  preceding  equations.  Figure  80  is  a  nomo¬ 
graph  of  the  thermal  expansion  stress,  Oj  ,  while  Figures  81  and  82  represent  nomographs 

of  the  thermal  buckling  stresses  in  the  x  and  ^-directions,  respectively.  Again,  the  para 
meter  v  -  0.32  was  used  to  develop  these  nomographs. 


4.  Ambient  Temperature  Fundamental  Mode  Frequency 


The  fundamental  mode  frequency  for  a  single  bay  of  a  multi-bay  structure  is  given  by 
Equation  (65b),  or 
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The  nomograph  corresponding  to  this  equation  is  presented  in  Figure  83,  for  v  -  0.32. 

The  chart  was  simplified  by  taking  advantage  of  the  essentially  constant  ratic^of  E/y 
for  most  aircraft  structural  alloys.  An  average  ratio  of  E/y  =  3.98  x  10^  in  /sec^  was 
used;  this  is  an  average  of  the  ratios  for  aluminum,  titanium,  stainless  steel,  and  Inconel 
al  loys . 


5.  Elevated  Temperature  Frequency  Response 

The  fundamental  mode  frequency  at  a  temperature  increase,  T,  is  given  by  Equation  (66), 

‘1/2 1 
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Figure  84  is  a  nomograph  of  this  relationship  to  simplify  the  computation. 


6 .  Dynamic  Stress 

Dynamic  stresses  at  any  temperature  can  be  computed  by  Equofion  (68): 


o  Rivet  row  at  midpoint  of  long  side 
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FIGURE  80.  THERMAL  EXPANSION  STRESS  NOMOGRAPH 
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FIGURE  82.  THERMAL  BUCKLING  STRESS  NOMOGRAPH 
Y  -  DIRECTION  STRESS 
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FIGURE  83.  AMBIENT  TEMPERATURE  FUNDAMENTAL  FREQUENCY  NOMOGRAPH 


The  elevated  temperature  response  frequency,  f(r),  must  be  used  for  these  computations. 

The  stresses  at  both  locations  must  generally  be  calculated  for  elevated  temperature  appli¬ 
cations  because  of  the  interaction  of  the  thermal  and  dynamic  stresses. 

The  stiffened  panel  fatigue  tost  data  of  Tables  XI  and  XII  were  used  to  establish  the  design 
life  nomograph  of  Figure  85  for  elevated  temperature  applications.  This  nomograph  is 
applicable  for  7075-T6  aluminum  at  temperatures  of  150°  and  300°F,  and  6AI-4V  titanium 
at  temperatures  of  400°  and  600°F. 

7.  Application  of  the  Design  Procedure 

At  least  two  alternative  methods  of  application  are  possible  using  the  criteria  developed 
on  this  program.  These  alternatives  are  described  in  the  following  subsections  in  the  form 
of  sample  applications. 

a.  Design  Life  Nomograph  -  The  most  direct  method  involves  the  use  of  the  design 
nomograph  for  life.  Figure  85.  The  alloys  and  structural  temperatures  must  agree  with 
those  for  which  the  nomograph  was  developed.  This  nomograph  includes  thermal  mean 
stress  effects  in  the  data,  thereby  negating  the  need  to  compute  these  parameters. 

EXAMPLE:  A  flat  aluminum  structure  is  to  be  designed  for  a  service  life  of  100  hours  at  a 
sound  pressure  spectrum  level  of  1  20  dB  and  a  service  temperature  of  300°F . 

Assume:  o  Aspect  ratio:  b/a=3.0(F^  =3.33) 
o  Damping  ratio:  C  =  0.016 
o  Ambient  temperature:  80°F 

As  o  first  step,  assume  o  frequency,  at  the  service  temperature,  of  300  Hz.  Then,  the  life 
in  cycles  is 

N  =  (300  Hz)(1Q0  hr)(3600  sec /Hr)  =  1  .08  x  10®  cycles 

Enter  the  nomograph  of  Figure  85  with  this  life  and  follow  through  the  parameters  to  the 
skin  thickness  chart.  Several  spacing/skin  thickness  ratios  are  now  possible,  oil  of  which 
will  meet  the  design  life  goal .  Assuming  a  skin  thickness  of  0.050  inch,  the  panel  width 
is  found  to  be  5.0  inches. 

At  this  point  a  structural  configuration  is  defined;  however,  the  assumed  frequency  must 
be  checked.  Compute  the  fundamental  mode  frequency  at  the  ambient  temperature  using 
Figure  83,  (f0  =  390  Hz)  then  compute  the  skin  buckling  temperature  increase  using 
Figuro  77,  whero  a  =  13.5  x  10"$  ln/ln/>F  from  Figure  V-2,  Appendix  V.  (T  ■  33°F) 

The  temperature  ratio  Is  then  c 

r  -  T/T  *  220/33  =  6.7 
c 
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Figure  84  then  yields  the  frequency  ratio  f(r)/f0  -  1  .65,  and  the  elevated  temperature 
frequency  is 


f(r)  =  1  .65  f0  -  645  Hz 


Since  this  frequency  is  greater  than  the  assumed  frequency,  the  anticipated  life  will  be 
less  than  the  design  goal,  and  an  iteration  is  necessary.  The  above  procedure  is  repeated 
using  the  calculated  frequency  of  645  Hz.  One  or  more  iterations  may  be  necessary  to 
obtain  agreement  between  the  initial  and  final  frequency  (or  design  life). 

Several  spacing/skin  thickness  ratio  combinations  may  be  derived  using  this  method,  and 
the  weight' of  each  calculated  to  obtain  a  minimum  weight  design. 

This  method  may  also  be  used  for  structural  temperatures  different  than  those  indicated 
on  the  nomograph  by  assuming  a  linear  relationship  between  the  temperatures  shown  and 
interpolating. 


b.  Mean  Stress  Fatigue  Curves  -  An  alternate  design  method  involves  the  use  of  fatigue 
curves  where  the  mean  stress  effects  are  known  (i  ,e. ,  Figure  V-4,  Appendix  V).  This 
method  can  be  used  where  the  alloy  or  temperature  does  not  coincide  with  those  of  the 
Figure  8.5  nomograph. 

EXAMPLE:  A  flat  structure  is  to  be  designed  for  a  service  life  of  100  hours  at  a  spectrum 
sound  pressure  level  of  140  dB,  and  a  service  temperature  of  500° F .  Stainless  steel 
PH15-7Mo  is  selected  as  the  alloy  to  be  used  for  this  structure. 

Assume:  o  Aspect  ratio:  b/a=3.0(F^  -3.33) 
o  Damping  ratio:  '  =  0.016 
o  Ambient  temperature:  8Q°F 

1  7 

A  fatigue  curve  for  the  selected  alloy  was  obtained  from  Ml  L-HD8K-5B  to  give  the 
effects  of  mean  stress  on  fatigue  life.  This  axial  loading,  constant  amplitude  fatigue 
curve  was  converted  to  an  equivalent  random  amplitude  fatigue  curve  (Figure  86)  using 
the  method  of  Reference  1 . 

Assume  a  stiffener  spacing  of  a  =  6  inch  and  a  skin  thickness  of  h  -  0.050  inch.  From 
Ml  L-HDBK-5B, 

V  -  ( •  2 77/  386)  -  7 . 1 7  x  10  lb  -sec  i n^ 
a  =  6.1  x  10"6  irvW0F  (£  500°F 
E0  =  29.0  x  106  pSi  £  RT 
E  -  26.97  x  106  psi  &  500° F 
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REF.  17,  FIGURE  2.5.7.1.8(b) 

CONVERTED  TO  FGU1VALENT  RANDOM  AMPLITUDE 


FIGURE  86.  RANDG  A  LOADING  FATIGUE  CURVE  FOR  PH15-7MC  STAINLESS  STEEL  AT  500°F 


The  skin  critical  buckling  temperature  is  found  from  Figure  78  as  Tc  -  50  F.  Then 
r  =  420/50  -  8.4  and  the  buckling  amplitude  i-  ',V  =  0.240  inch,  from  Figure  79. 

The  thermal  stresses  at  the  midpoints  of  the  two  sides  are  found  from  Figures  80  through 
82,  or 

~  -ICO  ksi 

o  =75  ksi;  a  =  40  ksi 
*b  yb 

then 

o  =  o_  +  a  ~  -100  +  75  =  -25  ksi 

x  T  x, 

b 

a  =  a _  +  a  -  -iCO  +  40  =  -60  ksi 
7  T  yb 

The  ombier.t  temperature  fundamental  frequency  is  fQ  =  260  Hz  from  Figure  83.  From 
Figure  84  the  frequency  ratio  corresponding  to  a  temperature  ratio  of  8.4  is  1  .8.  Then 

f(r)  =  1  .Sf  =  468  Hz 
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The  dynamic  stresses  are  then 
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The  value  of  i(f)  =  2.9  x  10-2  is  the  acoustic  pressure  density  corresponding  to  140  dB, 
while  AR  =  29.33  is  the  aspect  ratio  parameter. 


Enter  the  fatigue  cu've  of  Figure  86  (Kf  -  4)  with  the  dynamic  stress  =  7.89  ksirrns  and 
thermal  mean  stress  ax  ~  -25  ksi .  This  combination  gives  a  life  of  approximately  4 .5  x 
10°  cycles.  Using  the  y-direction  stresses,  3y  =  3.17  ksirms  and  3y  =  -60  ksi,  gives  a 
life  greater  than  10^0  cycles. 

At  a  frequency  of  468  Hz  the  life  is 


....  N 

LlrE=33ooT 


4.5  x  10 

(3600X468) 


=  267  Hours 


which  Is  greater  than  the  100  hour  design  requirement.  The  design  can  be  optimized  by 
iteration  on  the  above  procedure  to  decrease  the  skin  goge  or  Increase  the  stiffener 
spacing  such  that  the  predicted  life  Is  equal  to  or  greater  than  the  100-hour  design  life. 


1 
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8.  Limitations  in  the  Design  Procedure 


Application  of  these  design  procedures  should  be  tempered  with  a  thorough  understanding 
of  their  limitations.  Certain  of  the  initial  assumptions  stated  during  the  analytical  devel¬ 
opment  were  negated  by  derivation  of  the  empirical  relations.  However,  the  limits  of  the 
physical  and  environmental  parameters  tested  in  the  experimental  program  then  opply  to 
these  design  criteria.  These  limitations  are  itemized  below. 

a.  Physical  Constraints  -  The  bounds  of  the  test  specimen  dimensions  were  used  to 
establish  these  limitations.  These  should  be  treated  only  as  a  guide  as  the  equations  and 
nomog.aphs  are  normolly  valid  beyond  these  limits.  Individual  judgment  must  be  applied 
in  unusual  coses  where  the  constraints  are  drastically  exceedea ,  part'culari  ly  in  the  cose  of 
the  design  charts.  The  guidelines  on  size  are: 

o  Panel  bay  width:  a  =  5  to  9  inches 
o  Panel  bay  aspect  rotio:  b/a  -  1  .5  to  3 
o  Panel  skin  thickness:  h  =  0 .024  to  0 .063  i  nch 

b.  Environmental  Constraints  -  The  acoustic  environment  generally  has  no  restrictions 

as  regards  applicability  of  the  design  criteria.  Spectrum  levels  below  120  dB  will  normally 
result  in  low  dynamic  stresses  and  a  long  fatigue  life.  The  higher  noise  levels  will  gen¬ 
erally  result  in  nonlinear  response,  but  these  effects  are  included  in  the  design 
criteria,  since  many  of  the  test  specimens  exhibited  a  high  degree  of  nonlinearity. 

The  thermal  environment  must  be  nearly  uniform  over  the  surface  of  a  panel  bay.  The  skin 
temperature  is  limited  to  the  maximum  temperature  for  which  the  alloy  retains  significant 
structural  properties.  The  limiting  temperatures  are,  for  the  alloys  considered: 

o  70  75-T6  aluminum:  30Q°F  maximum 
o  6AI-4V  annealed  titanium:  600°F  maximum. 

The  design  life  criteria  are  based  only  on  specific  nominol  temperatures,  requiring  the 
use  of  interpolation  for  intermediate  temperatures.  Extrapolation  beyond  the  temperature 
limits  may  be  permissible  to  some  extent  if  care  is  exercised  and  the  further  temperature 
degradation  effects  are  included. 

Probably  the  most  important  restriction  on  the  design  method  is  in  the  estimation  of  the 
ambient  temperature  and  the  state  of  the  structure  at  this  temperature.  All  thermal  res¬ 
ponse  relations  are  referenced  to  the  ambient  temperature  and  the  assumption  that  a  state 
of  stress  equilibrium  exists  (i  .e  . ,  r  ,  mean  s>resses).  It  is  impractical  at  this  stage  to  give 
guidelines  for  estimating  the  ambient  temperature  state,  because  it  is  dependent  on  the 
length  of  time  at  a  uniform  temperature,  external  constraints,  and  other  influences.  It 
will  be  noted  that  a  change  in  ambient  temperature  over  o  short  time  interval  constitutes  o 
temperature  change  as  far  as  the  analytical  relations  are  concerned. 
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c  ■  External  Constraint  -  The  external  constraints  imposed  on  the  test  panels  precluded 
significant  thermal  expansion  of  the  substructure.  This  is  considered  representative  of 
structural  applications  in  the  direct  flow  path  of  engine  exhausts  or  other  heat  sources, 
where  only  localized  areas  of  the  structure  are  heated.  The  criteria  can  also  be  applied 
to  design  applications  involving  gradual  heating  of  an  entire  structural  area,  where  all 
structure  expands  at  about  the  same  rate.  This  corresponds  to  a  relaxation  of  the  external 
constraints  from  those  considered  here.  In  this  cose,  the  thermal  buckling  amplitudes  and 
stresses  given  by  the  relations  delineated  herein  will  result  in  a  conservative  design. 

It  should  be  noted  that  the  empirical  results  presented  herein  are  applicable  only  for  the 
case  of  simultaneously  application  of  heat  and  noise.  Alternate  application  of  these 
environments,  wherein  significant  thermal  stress  cycles  are  incurred,  were  not  considered 
in  this  program . 


134 


VI  -  CONCLUSIONS 


Design  methods  were  developed  to  estimate  the  acoustic  fatigue  life  of  aircraft 
structure  exposed  to  simultaneous  noise  and  elevated  temperatures.  These  design 
methods  are  applicable  to  flat,  stiffened-skin  structures  fabricated  of  aluminum  or 
titanium  alloys.  The  methodology  wos  derived  through  a  combined  analytical/experimental 
program,  wherein  the  analysis  served  to  identify  parameters  important  to  the  experi¬ 
mental  effort.  Analytical  results  are  presented  for  a  single  panel  and  for  a  multi-boy 
panel  subjected  to  a  spatially  uniform  temperature  rise.  This  analysis  covers  the  pre- 
and  post-buckled  states  for  each  structural  model.  Coupon  fatigue  tests  were  conducted 
to  isolate  temperature  degradation  effects  on  the  fatigue  life  of  each  alloy.  A  total  of 
27  aluminum  and  21  titanium  stiffened  panels  were  subjected  to  a  thermal  environment 
and  tested  to  failure  under  high  intensity  random  noise.  The  data  from  the  test  program 
were  used  to  modify  the  analytical  results  and  to  provide  design  equations,  nomographs, 
and  a  computer  program  for  predicting  acoustic  fatigue  life.  The  following  conclusions 
were  drawn  from  the  results  of  this  investigation; 

a.  The  analytical  and  experimental  investigation  results  indicate  that  the  overall 
rms  stress  response  of  stiffened  structure  skins  to  acoustic  excitation,  at  ambient  or 
elevated  temperatures,  can  be  estimated  using  empirical  relations  based  on  thin  plate 
and  Miles'  response  theories.  However,  significant  voriance  between  estimated  and 
measured  stresses  can  be  expected  because  of  cumulative  errors  in  predicting  the 
excitation  and  response  characteristics. 

b.  Existing  substructure  design  criteria®  for  ambient  temperature  structures  are 
considered  valid  for  elevated  temperature  applications,  where  the  skin  gages  are  deter¬ 
mined  from  the  design  criteria  of  this  report. 

c.  The  assumption  that  the  substructure  was  thermally  isolated  from  the  skin  proved 
adequate  for  establishing  empirical  criteria  for  skin  thermal  buckling  effects.  However, 
the  analytical  development  should  be  extended  to  include  heating  of  the  substructure  to 
lefine  the  elevated  temperature  response  characteristics. 

d.  The  thermal  buckles  experienced  on  certain  structural  designs  may  prove  excessive 
for  operational  use;  hence  buckling  amplitude  should  be  a  significant  consideration  in 

the  design  of  elevated  temperature  structures. 

e.  Structural  buckling  temperatures,  amplitudes,  and  stresses  con  be  estimated  using  the 
empirical  results  of  this  program.  Since  each  of  these  parameters  is  based  on  a  temperature 
increase  relative  to  ambieit,  the  confidence  in  the  estimated  thermal  and  dynamic  response 
will  be  directly  related  to  the  accuracy  of  the  ambient  temperature  predictions. 

f.  The  use  of  faying  surface  sealant  between  the  skin  and  substructure  has  no  dis- 
cernable  effect  on  acoustic  fatigue  life,  at  ambient  or  at  elevated  temperatures. 

g.  The  analyses  presented  for  the  elevated  temperature  dynamic  response  of  box  and 
curved  structures  in  Appendices  I  and  II  are  preliminary  and  should  be  used  with  caution 
until  further  development  and  verification  by  experimental  data. 
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APPENDIX  I 


VIBRATION  ANALYSIS  OF  NINE-CELL  BOX  STRUCTURE  AT  ELEVATED  TEMPERATURES 


The  simple  panel  analysis  described  in  Section  II  is  based  on  the  assumption  that  the 
temperature  distribution  over  the  panel  is  uniform  and  in  a  steady  state.  Based  on  these 
assumptions,  the  application  of  the  simple  panel  theory  to  a  flat,  stiffened,  nine-boy 
panel  structure  was  presented.  For  the  nine-bay  stiffened  ponel,  the  assumption  of  a 
f  uniform  temperature  distribution  may  not  be  too  restrictive,  but  if  the  same  assumption  is 

applied  to  a  box  structure  configuration,  the  results  may  be  considered  unrealistic  . 

As  was  shown  in  the  nine -bay  panel  analysis,  this  assumption  of  a  uniform  temperature 
distribution  implies  thot  thermal  equilibrium  exists  at  each  temperature  considered .  That 
is,  either  the  stiffeners  experience  the  same  temperature  rise  as  the  cover  sheet,  or  the 
stiffeners  ore  insulated  from  the  cover  sheet  so  that  no  heat  conduction  can  occur.  Since 
this  analysis  did  not  consider  a  thermal  stress  in  the  stiffeners,  the  assumption  of  a  uniform 
temperature  distribution  in  the  cover  sheet  for  the  nine-bay  stiffened  panel,  as  applied  in 

■  j  Section  II,  implies  that  the  stiffeners  are  insulated  from  the  cover  sheet  and  that  thermal 

radiation  from  the  cover  sheet  to  the  stiffeners  con  be  ignored  in  the  analysis.  Since  the 
stiffeners  do  not  present  a  large  surface  area  when  compared  to  the  cover  sheet  area,  and 
,  since  the  temperature  range  considered  is  rather  low  (less  than  600°F),  neglecting  the 

■  thermal  radiation  effects  for  the  stiffened  nine-bay  panel  may  be  justified.  The  effect  of 
these  assumptions  os  applied  to  the  nine-cell  box  structure  will  nowbe  briefly  discussed  . 
This  analysis  is  presented  only  as  a  beginning  point  for  this  complex  structure;  no  experi- 

*  mental  effort  has  been  devoted  to  elevated  temperature  fatigue  testing  of  nine-cell  box 

structures. 

Consider  the  nine-cell  box  structure  configuration  os  illustrated  in  Figure  1-1  ,  For  this 
configuration,  the  two  cover  sheets  (located  in  the  planes  z  =  0  and  z  =  h|)  ore  of  equol 
surface  area,  and  the  surface  areas  of  the  ribs  are  of  the  some  order  of  magnitude  as  the 
areas  of  the  cover  sheets.  It  can  thus  be  supposed  that  thermal  radiation  from  one  element 
.  [  to  another  will  be  significant  in  determining  the  temperature  of  the  element .  Also,  it  is 

evident  that  if  the  two  cover  sheets  are  maintained  at  different  temperatures  the  ribs  will 
.  not  experience  a  uniform  temperature  distribution  unless  the  ribs  are  insulated  from  the 

[  cover  sheets. 

The  analysis  of  the  uniformly  heoted,  simply  supported  panel  will  now  be  used  to  obtain 
,  an  estimate  of  the  response  frequency  of  the  box  structure  when  subjected  to  a  temperature 

distribution  such  thot  the  cover  sheets  and  the  ribs  are  each  maintained  ot  a  specified 
,  uniform  temperature.  The  assumed  temperature  distribution  and  the  dimensions  of  the 

components  are  illustrated  in  Figure  I -I  . 

As  discussed  in  Reference  8,  three  mode  numbers  are  required  to  describe  the  vibration  of 
a  box  structure.  The  mode  numbers  are  denoted  by  the  nomenclature  (m,  n,  p),  where 
(m,n)  are  the  mode  numbers  for  the  cover  sheets  in  the  x  and  y  directions,  respectively, 
and  p  denotes  the  mode  number  of  the  rib  across  the  depth  of  the  structure. 
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FIGURE  1-1.  NINE-CELL  BOX  STRUCTURE  CONFIGURATION 


Using  the  results  of  Equation  (3),  imposing  the  slope  compatibility  relations  (Equations  43, 
Reference  8),  and  neglecting  the  detailed  algebra,  the  modal  mass  the  box  structure  for 
the  (1,1,1)  mode  is 
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Equation  (1-1)  represents  the  modal  mass  both  below  and  above  the  critical  temperature, 
since  the  mass  is  not  dependent  on  temperature. 


A.  Pre-Buckled  Frequency  Response 

The  modal  stiffness  of  the  nine-celi  structure  below  the  critical  temperature  is  obtained  by 
using  Equation  (34a)  and  the  slope  compatibility  relationships  between  the  cover  sheet 
bays  and  the  rib  bays.  It  is  first  assumed  that  all  components  are  manufactured  from  the 
same  material .  As  discussed  previously,  the  uniform  temperature  of  the  cover  sheet 
located  in  the  plane  z  =  0  is  denoted  by  Tj ,  the  uniform  temperature  of  all  ribs  is  denoted 
by  Tr,  and  the  uniform  temperature  of  the  cover  sheet  located  in  the  plane  z  =  h^  is 
denoted  by  Tj. 

For  these  assumptions,  the  modal  stiffness  for  the  (1,1,1)  mode  is  given  os 
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where  the  terms  are  defined  in  Table  l-l .  ^t  should  be  noted  that  the  temperature  rise,  r, 
as  defined  in  Equation  (1-2)  and  Table  l-l, is  the  ratio  of  the  temperature  rise  of  an  indi¬ 
vidual  component  to  the  critical  temperature  of  the  center  bay  of  the  cover  sheet  located 
in  the  plane  z  =  0.  That  is,  the  critical  temperature  for  the  entire  structure  is  defined  by 
the  buckling  temperature  of  the  center  bay  of  the  cover  sheet  located  in  the  plane  z  =  0, 
whether  or  not  this  happens  to  be  the  lowest  buckling  temperature  of  all  of  the  panel  bays 
of  the  structure . 
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TABLE  I -I 


TERMS  APPEARING  IN  THE  MODAL  STIFFNESS  EXPRESSION 
EQUATION  (1-2) 
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B.  Post-Buckled  Frequency  Response 


The  analysis  developed  in  Section  II  for  a  heated,  simply  supported  panel  above  the 
critical  temperature  will  now  be  used  to  derive  the  expression  for  the  modal  stiffness  of  a 
nine -cel  I  box  structure  above  the  critical  temperature.  As  for  the  nine-boy  panel,  the 
critical  temperature,  rQ,  is  defined  as  the  temperature  rise,  r,  necessary  to  cause  the 
modal  stiffness  expression,  Equation  (1-2),  to  vanish. 


For  a  rectangular  flat  panel  with  dimensions  a;  x  bj,  and  thickness,  t,  exposed  to  a 
uniform  temperature  distribution,  T,  the  strain  energy  of  the  panel  is  expressed  as 
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In  these  equations,  the  temperature  parameter,  r,  is  referenced  to  a  panel  with  dimensions 
02  x  b2  and  thickness  ]ts  which  is  exposed  to  a  uniform  temperature  increase  Tj  .  Also, 
;;W0  is  the  static  buckled  amplitude  of  the  panel  and  q; j(t)  is  the  dynamic  amplitude  of 
trie  panel  (jq;:|  <<:  ijW0).  In  the  derivation  that  follows,  the  analysis  is  limited  to  the 
box  structure  fundamental  mode  (m,n,p)  =(1,1,1). 

Following  the  analysis  developed  in  Reference  8  for  box  structure  at  ambient  temperature, 
it  is  assumed  that  all  components  of  the  box  structure  ore  manufactured  from  the  same 
material  and  that  the  structural  geometry  I*  as  Illustrated  in  Figure  I-) .  Using  the  above 
equations  and  again  Imposing  the  slope  compatibility  relations,  the  strain  energy  expres¬ 
sion  for  the  box  structure  for  the  (1,1,1)  mode  Is 
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where  K**  and  R**  are  defined  in  Table  I -II.  In  Equation  (1-8),  W0  is  the  static  buckled 
panel  amplitude  and  q(t)  is  the  dynamic  panel  amplitude  of  the  center  panel  bay  of  the 
cover  sheet  located  in  the  plane  z  =  0  . 

It  is  seen  from  Table  I -II  and  Equotion  (1-4)  that  the  stiffness  is  a  function  of  the  rotio 
(Wo/l  t,)2.  To  determine  this  parameter  as  a  function  of  the  temperature  increase  Tj,  or 
r  =  T,AC,  (where  Tc  is  the  critical  temperature  of  the  center  bay  of  the  cover  sheet 
located  in  the  plane  z  -  0)  set  q(t)  =  0  in  Equation  (1-8)  ond  obtain  the  strain  energy  as  a 
function  of  WQ .  Minimizing  this  result  with  respect  to  WQ  gives  the  result 
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Solving  Equation  1-9  for  (W  /|  t  )  as  a  function  of  r  =  Tj/Tc/  and  then  substituting  the 
result  into  the  expressions  for  K**  and  R**  gives  the  stroin  energy  expression  in  terms  of 
the  temperature  parameter,  r. 


The  modal  stiffness  is  then  derived  by  minimizing  Equation  1-8  with  respect  to  q(t)  so  that 
the  expression  for  the  modal  stiffness  of  the  fundamental  mode  is 
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where  the  terms  (Wq/^ appearing  in  K** 


have  been  determined  from  Equation  (1-9). 
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TABLE  I -I  I 


TERMS  APPEARING  IN  THE  STRAIN  ENERGY  EXPRESSION 
EQUATION  (1-8) 
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i  s 


2'  2'i's'  ~\  a„  /  '  '“2'“1  'i  s' 
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r  2' 


I  S  I 
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\a2b2< 


K*(b.  ,a . ,  .t  ,  T.) 

1  I  I  S  i 


iR;  =  R*<b2'ViW  +  2(^)R*<b2'arit$,Ti) 

t2(q)R,(b1'V1V1i>^(^)R*<br°r!YTi> 


lKr  ’  K*ll,rb2'lV  V  *  2(b^)K*(hrbri\'V 


2K*  =K*(h,,=2,2VTr,o(^)K-(h,.o1,2VTt) 


iR*  =  R'(brb2'iVV  *  2l  b^R'(VbriVV 


2R r  '■  R*(h1'°2'2,r/Tr)  +  2(  a2  r  (hl'V2W 
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APPENDIX  li 


APPROXIMATE  FREQUENCIES  Of  HEATED  CYLINDRICAL  PANELS 
1 8 

Szechenyi  piesenfs  on  approximate  technique  for  estimating  the  natural  frequencies  of 
curved  panels.  Szechenvi's  result  has  been  checked  against  analytical  and  experimental 
data  oresented  by  Petyt  with  good  agreement.  The  following  development  is  based  on 
tne  Reference  18  technique. 

For  either  simply  supported  or  clamped  boundaries  of  the  rectangular  cylindrical  shell 
illustrated  in  Figure  1 1  —  1 ,  Szechenyi  obtains  the  expression  for  the  frequency  of  the  (m,  n) 
mode  as 

f 

mn  v- 
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(11-1) 


where  m  is  the  mode  number  in  the  longitudinal  direction  and  n  is  the  mode  number  in  the 
circumferential  direction.  The  stresses  ax  and  are  the  iongitudinal  panel  edge  stress 
end  the  circumferential  panel  edge  stress,  respectively. 


The  terms  k^,  kn  and  G*  depend  upon  the  boundary  conditions  of  the  panel  and  are 
defined  for  simply  supported  edges  as 


k  =  Mn/a 
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k  -  nn/b 

r  i 


G*  =  0  - 


v?)  k4 
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m  n 


I  or  lamped  edgos,  the  factors  km  ond  kn  are 


(11-2) 


k  -  nm/6  a  ;  h  -  nrr/6  L 
mm  n  n 


(11-3) 


Vuiues  of  G*  can  be  obtained  from  Figure  M-2,  and  values  o  6  and  6n  con  be  obtained 
ft  n  Figure  l(-d.  m 


lr.  ter.TB  of  an  interna)  pressure,  p,  and  a  uniform  tempercture  rise,  T,  the  edge  stress 
c  -  .iponenls  are 
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FIGURE  11-2.  VARIATION  OF  PARAMETER  G*  WITH  RATIO  k  /k 

m  n 
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(H-5) 


substituting  Equations  (1-4)  into  Equations  (II  — 1)  gives  the  result 
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Setting  the  frequency  expression  given  by  Equation  (11-5)  to  zero,  the  expression  for  the 
critical  buckling  temperature  is  then 


T 

c 


1 2cr  ( 1  +  v) 

h2(k2  +k2) 
m  n 


(Pr  k2/D)  +(12  GVh2r2) 
c  n  c 

(k  2  +  k  2  )  2 
m  n 


(H-6) 


Hence,  Equations  (I I -6)  and  (1 1-5)  may  be  used  to  estimate  the  buckling  temperature  and 
frequency  response,  respectively,  for  curved  stiffened  structure.  It  should  be  emphasized 
that  these  results  are  tentative,  since  no  experimental  verification  has  yet  been  attempted. 


APPENDIX  III 


TEST  SPECIMEN  DETAILS  AND  TEST  INSTRUMENTATION 

The  stiffened  panels  used  for  the  second  phase  of  the  experimental  program  are  briefly 
described  in  Section  III.  That  section  also  describes  the  test  set-up  and  procedures 
used  for  each  test.  This  appendix  contains  a  more  detailed  description  of  the  test 
specimens,  together  with  schematic  diagrams  of  the  instrumentation  used  for  each  test. 

It  supplements  (rather  than  duplicates)  the  discussion  in  Section  III. 

A.  Stiffened  Panel  Specimens 

The  aluminum  and  titanium  test  panels  were  very  similar  in  design,the  only  significant  dif¬ 
ference  being  in  the  substructure.  The  channel  frames  and  zee  stiffeners  on  the  aluminum 
specimens  were  hot-rolled  to  a  bend  radius  of  approximately  three  times  the  thickness, 
whereas  the  titanium  substructure  was  welded  at  the  longitudinal  corners.  The  general  plan 
and  edge  view  of  the  stiffened  panels  is  presented  in  Figure  lll-l .  Figures  MI-2  and  IM-3 
show  details  of  the  frame  and  stiffener  members  for  the  aluminum  and  titanium  specimens, 
respectively.  The  edge  member  used  to  support  the  test  specimens  in  the  steel  test 
fixture  is  detailed  in  Figure  111-4.  The  edge  member  shown  is  for  the  titanium  specimens 
and  was  fabricated  of  6A1-4V;  the  aluminum  edge  member  was  similar  but  made  of 
7075-T6.  Only  four  edge  members  of  each  alloy  were  fabricated.  They  were  removed 
from  the  specimens  at  the  completion  of  testing  and  were  installed  on  the  next  set  of 
panels  to  be  tested.  Use  of  two  sets  of  two  edge  members  each  allowed  installation  and 
checkout  of  specimens  to  proceed  concurrent  with  testing. 


The  coupon  fatigue  testing  described  in  Section  III  was  accomplished  using  fixed  end, 
cantilever  beam  test  specimens.  The  specimen  support  was  provided  by  the  clamp 
blocks  shown  in  Figure  8  and  detailed  in  Figure  111-5.  Room  temperature  testing  was 
accomplished  using  phenolic  inserts  between  the  specimen  and  the  steel  blocks.  Elevated 
temperature  resting  used  stainless  steel  inserts  in  place  of  the  phenolic  with  asbestos 
insulation  between  the  steel  and  the  specimen  to  minimize  heat  flow  into  the  clamp 
blocks  from  the  test  specimen. 

The  instrumentation  set-up  used  for  the  coupon  fatigue  testing  is  shown  schematically  in 
Figure  1 1 1-6.  This  depicts  both  the  sinusoidal  frequency  input  as  well  as  the  random 
amplitude  input  used  for  fatigue  testing. 

C.  Stiffened  Panel  Tests 

The  test  instrumentation  discussed  in  the  following  subsections  parallel  the  discussion 
of  Section  III. 


150 


PLAN  VIEW  SECTION  A-A 


NOTE:  EDGE  MEMBER  OMITTED  FOR  CLARITY 

FIGURE  lll-l  STIFFENED  SKIN  SPECIMEN  CONFIGURATION 

bl 


FIGURE  111-2.  FRAME  AND  STIFFENER  DETAILS 
ALUMINUM  SPECIMENS 
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SECTION  A-A  TYPICAL  FRAME-STIFFENER  INTERSECTION 

FIGURE  1 1 1-3.  FRAME  AND  STIFFENER  DETAILS 
TITANIUM  SPECIMENS 


FIGURE  II 1—4.  STIFFENED  PANEL  EDGE  MEMBER 
TITANIUM  SPECIMENS 
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FIGURE  111-6.  SCHEMATIC  DIAGRAM  OF  COUPON  FATIGUE  TEST 
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1. _ Room  Temperature  Frequency 

The  room  temperature  fundamental  mode  frequency  test  set-up  was  shown  in  Figure  25. 
frequericies  were  determined  by  making  sinusoidal  frequency  sweeps,  at  constant 
input,  and  plotting  the  strain  response.  The  instrumentation  schematic  diagram  for  this 
is  shown  in  Figure  111-7. 

,2. _ Temperature  Effects  on  Frequency 

The  test  set-up  described  above  was  also  used  to  determine  the  elevated  temperature 
frequency  response.  Six  heat  lamps  w«"-e  positioned  above  the  specimens  as  shown  in 
Figure  26  and  the  frequency  was  the: i  rranuully  varied  to  maintain  a  fundamental  mode 
response.  The  schematic  diagram  of  the  set-uo  is  identical  to  that  of  Figure  lil-7, 
except  that  strain  response  plots  were  not  produced.  Skin  temperatures  were  manually 
read  from  a  single  channel  indicator. 

3.  Thermal  Strain  and  Deflection 


The  test  set-up  for  these  measurements  is  shown  in  Figures  27  and  28  and  shown 
schematically  in  Figure  1 1 1— 8 .  Temperatures  were  controlled  manually  and  read  from 
a  multi-channel  recorder.  Displacements  and  thermal  strains  corresponding  to  these 
temperatures  were  manually  read  and  recorded  on  data  sheets. 

4.  Excitation  Spectrum  Shaping  and  Fatigue  Tests 

The  specimens  were  subjected  to  high  intensity  noise,  and  strain  response  plots  were 
produced  at  the  desired  test  temperature.  The  schematic  diagram  for  this  test  is  shown 
in  Figure  111-9.  After  determining  the  significant  panel  response  modes,  the  desired 
random  spectrum  shape  was  obtained  using  the  equipment  shown  in  the  schematic  diagram  of 
Figure  111-10.  This  also  shows  the  instrumentation  used  to  record  noise  and  strain  on 
magnetic  tape  for  later  analysis. 

Three  types  of  data  analyses  were  accomplished:  (a)  narrow-band  analyses  of  the  noise 
spectrum,  (b)  narrow-band  analyses  of  the  strain  response,  end  (c)  probability  density 
analyses  of  the  noise  and  strain.  The  schematic  diagram  of  the  equipment  necessary  to 
accomplish  the  first  of  those  is  shown  in  Figure  111-11,  while  Figure  lil-12  shows  the  strain 
narrow-band  analysis  instrumentation.  The  only  difference  in  the  two  types  of  analysis 
was  in  the  log  converter  used  for  the  noise  analysis  to  obtain  a  dB  scale.  The  probability 
density  analysis  schematic  diagram  ;s  shown  in  Figure  111-13;  noise  and  strain  signals 
were  analyzed  with  the  same  instrumentation. 
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FIGURE  111-6.  SCHEMATIC  DIAGRAM  FOP  THERMAL  STRAIN 
AND  BUCKLING  AMPLITUDE 
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FIGURE  111-9.  SCHEMATIC  DIAGRAM  FOR  HIGH  INTENSITY  NOISE  FREQUENCY  S 
STRAIN  RESPONSE 


FIGURE  I ; : - 1  2 .  SCHEMATIC  DIAGRAM  FOR  NARROW-BAND  ANALYSIS  OF  STRAIN 


APPENDIX  IV 

CORRELATION  OF  EXPERIMENTAL  RESULTS  WITH 
NINE-BAY  ANALYTICAL  RESULTS 


The  analytical  section  presented  both  a  simple  panel  and  a  nine-bay  panel  analytical 
development.  The  results  of  the  simple  panel  analysis  were  correlated  with  the  test  data 
to  yield  the  empirical  equations  of  Section  IV  and  the  design  nomographs  of  Section  V. 

The  nine-bay  analytical  development  resulted  in  more  complex  equations  which  are  not 
readily  solved  without  the  aid  of  a  computer.  The  nine-bay  analytical  results  were  corre¬ 
lated  wi  th  the  ex  peri  menta  I  data,  and  a  digital  computer  progrcm  was  developed  to 
simplify  solution  for  routine  design  problems. 


A.  Room  Temperature  Fundamental  Frequency 


The  room  temperature  fundamental  mode  frequency  of  the  center  bay  is  given  by  Equation 
(46)  with  r  -  0,  or  . 

f  -  1  f  Kj(0)  1  '  Z  (IV-1) 


2v. 


M 


The  mass,  M| ,  is  the  total  combined  mass  of  the  skin  and  supporting  structure,  as  defined 
by  Equation  (43).  It  was  assumed  for  this  development  thot  the  analytical  mass  relation 
was  realistic,  since  only  alloy  density  and  specimen  configuration  ore  involved.  The 
structural  stiffness  is  more  difficult  to  represent  analytical  ly;  hence,  the  measured  and 
calculated  frequencies  were  correlated  to  provide  an  empirical  modification  of  the  total 
stiffness . 

The  stiffness  is  defined  by  Equation  (44)  which,  for  r  =  0,  becomes 


K  (0)  =  rr  D  F*  +  K  +■  K 
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^a2^2 


sx 


sy  I 


(IV-2) 


The  parameter  represents  the  skin  stiffness,  arid  and  Ks>/  represent  the  substructure 
stiffness . 


sy 


i 


fl.u  usu  of  these  equations,  without  modification,  to  calculate  narural  frequencies  for 
the  test  specimen  configurations  produced  frequencies  15  to  20  times  higher  than  the 
measured  vaiucs.  Comparison  of  the  skin  and  substructure  terms  revealed  that  the  sub¬ 
structure  teims  were  3f>  to  1435  tiroes  higher  Irian  ’he  skin  stiffness  terms. 

A  constant  C  |  was,  therefore,  introduced  into  the  stiffness  K  j (0)  as  follows: 


16'' 


K  (0)  — 

1  4°2b2 


F.  ,K  +  K 
r,  +  sx  sy 

C 


( I V — 3) 


The  constant  C  ]  was  varied  from  100  to  1200,  the  frequencies  calculated  and  plotted  versus 
measured  data,  and  standard  deviation  and  correlation  coefficient  computed  for  each  case. 
The  constant  C]  =  400  was  selected  as  the  value  which  minimized  the  standard  deviation 
(or  maximized  the  correlation  coefficient).  Figure  IV  —  1  shows  the  correlation  of  calcu¬ 
lated  and  measured  frequencies  for  Cj  =  400,  and  includes  the  least  squares  regression 
line  (A),  The  approximate  regression  line  (B) ,  plotted  through  the  origin  and  the  data 
centroid,  was  selected  for  the  empirical  equation,  since  the  differences  between  the  two 
lines  are  insignificant.  The  following  empirical  equation  for  the  natural  frequency  results: 


f  =  1.03 
e  2rt 


K,(0) 


=  0.164  [  Kl(0) 


(IV— 4) 


where  Mi  remains  as  defined  and 


Ft  +  .0025  (K  +  K  ) 
1  sx  sy 


(IV— 5) 


Comparison  of  Figure  IV  —  1  with  the  same  plot  for  the  simple  panel  (see  Figure  71, 
Section  IV)  shows  the  following  differences  in  statistical  properties: 

Standord  Correlation 

Deviation  Coefficient 


Simple  Panel  Equation 


20.69 


0.915 


Nine-Bay  Panel  Equation 


27.03 


0.850 


Although  there  is  a  slightly  greater  variance  in  the  nine-bay  results,  frequency  calcula¬ 
tions  using  either  method  will  provide  comparable  predictions. 


B.  Buckling  Amplitude 

J 

The  measured  buckling  amplitudes  were  correlated  with  analytical  buckling  amplitudes 
given  by  Equation  (51),  or 

W  _  [Fo  r  -  F.  -  (K  +K  )  l'/2  . 

o  -  2  1  sx  sy  (IV-6) 

T  - R* -  (r  *  r0> 


From  Equation  (57) 


F*  +K  +  K  =  F*  r 
I  sx  sy  2  o 
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MEASURED  FREQUENCY 


FIGURE  I V- 1  .  FUNDAMENTAL  t  ODE  FREQUENCY  CORRELATION 
NINE-BAY  ANALYSIS,  AMBIENT  TEMPERATURE 
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(IV— 7) 


and  the  buckling  amplitude  becomes,  upon  substitution, 


W 

o 
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As  defined,  the  term  rQ  is  the  temperature  ratio  at  which  the  strain  energy  becomes  zero. 
It  is  thus  a  parameter  by  which  the  simple  panel  buckling  temperature,  tc  ,  may  be 

multiplied  to  obtain  the  buckling  temperature,  Tc  ,  of  the  center  bay  of  a  multi-bay 

m 

panel  (where  the  sizes  of  the  two  panels  are  identical).  Hence,  the  parameter,  r,  in 
the  above  equation  is  the  temperature  ratio  of  an  equivalent  simple  panel.  For  instance, 
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where  r  is  the  temperature  ratio  for  the  center  bay  of  a  nine-bay  (or  multi-bay)  panel. 
S;nce  the  empirical  expression  for  buckling  temperature  is  based  on  multi-bay  panel  test 
data,  the  buckling  amplitude  can  be  redefined  as 


(r  ;>  0 
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( I V  -  8) 


where  r 

m 


T.  is  given  by  Equation  (60)  or  the  nomograph  of  Figure  78. 
m 


Calculated  displacement  ratios  were  plotted  versus  measured  displacement  ratios  for  each 
panel  configuration,  and  the  slope  of  the  regression  line  computed  for  each  plot.  The 
analytical  expression  dod  not  match  the  measured  data  for  all  aspect  ratios;  therefore,  the 
'■egression  line  slopes  were  plotted  versus  the  various  parameters  in  the  analytical  equation 
and  found  to  correlate  best  with  the  temperature  parameter,  r0.  This  result,  shown  in 
Figure  IV-2,  revealed  a  trend  of  decreasing  slope  with  increasing  values  of  the  parameter 
rQ.  Both  linear  and  exponential  curve  fits  were  plotted  through  the  data  points  to  deter¬ 
mine  the  best  analytical  representation.  Both  equations  were  used  to  calculate  the  buck¬ 
ling  amplitude,  and  the  resulting  curves  were  compared  with  the  measured  displacement 
ratios.  The  linear  equation  was  found  to  produce  the  best  agreement  with  the  test  data; 
hence,  the  displacement  ratio  for  the  multi-bay  panel  is  expressed  by  the  following 
empirical  relation: 
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T-— -  =  (3.37  -  0.20  r  ) 
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Unlike  the  simple  pane!  results,  the  buckling  amplitudes  given  by  this  equation  are  not 
constant  for  fixed  aspect  ratios  since  the  parameter  r0  varies  with  changes  in  substructure 
or  skin  stiffness. 


167 


-  SLOPE  OF  APPROXIMATE  REGRESSION  LINE 


C  .  Thermal  Stress 


The  thermal  stresses  for  a  multi-bay  panel  are  defined  by  the  same  relations  used  for  the 
simple  panel.  Since  the  only  parameter  that  differs  between  the  two  forms  of  analysis  is 
the  buckling  amplitude/  thermal  strains  due  to  skin  buckling  were  isolated  foi  comparison 
with  the  measured  data.  This  is  identical  to  the  correlation  describe  !  for  the  simple 
panel  in  Section  IV . 

Measured  buckling  strains  were  plotted  versus  calculated  buckling  strains  as  shown  in 
Figures  I V —3  and  IV-4.  The  data  for  the  x-direction  strains/  Locations  1  and  2,  were 
combined  as  shown  in  Figure  !V-3.  The  slope  of  the  approximate  curve  fit  was  used  to 
develop  the  following  thermal  stress  relations: 

o  Midpoint  of  panel  bay  long  side  (y  =  b/2) 

(IV- 10) 


(IV- 1 1) 


The  results  of  this  analysis  compare  favorably  with  the  simple  panel  thermal  strains 
presented  in  Section  IV. 
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D.  Elevated  Temperature  Frequency  Response 


The  elevated  temperature  frequency  response  is  given  by  Equations  (59)  os 


f(r)  =  f 

o 


(0  <  r  <  r  ) 
-  -  o 


(r  >r  ) 

o 


As  explained  previously/  r  is  the  temperature  ratio  for  on  equivalent  simple  panel; 
substituting  the  multi-bay  temperature  ratio,  rm,  for  r  gives 


f(r)  =  f 
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(  0  <  r  <  I)  (IV-12) 
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£c  -  ANALYTICAL  STRAIN  -  n  IN/IN 


THERMAL  BUCKLING  STRAIN  CORRELATION 
LOCATION  1  &  2  -  PANEL  LONG  LENGTH  MIDPOINT 
NINE  -  BAY  ANALYSIS 


FIGURE  IV -4 .  THERMA.  B'C'rOING  STRAIN  CORRELATION 

LOCATION!  o  PANEL  SHORT  LENGTH  MIDPOINT 
NINE  -  BAY  ANALYSIS 


These  equations  ore  identical  to  the  simple  panei  analytical  frequency  relations  given 
in  Section  II,  where  the  temperature  ratio  for  the  multi-boy  panei  is  now  u>ed  m  lieu 
of  the  temperature  ratio  for  the  simple  pu.nel,  hence,  the  empirical  equations  givon 
previously  in  Section  IV  are  valid  for  the  multi-bay  panel. 

£.  Computet  Program 

Since  the  nine-bay  panel  analysis  involves  extensive  computation,  a  digital  computer 
program  was  developed  to  simplify  solution.  The  program  was  developed  for  the  Univoc 
1  106  computer,  in  Fortran  V;  however,  the  program  can  readily  be  adopted  to  any 
d igi ta I  computer .  The  input  data  format  ir.  shown  in  Table  IV -I ,  whi  le  Table  IV -  I  I 
contains  a  definition  of  input  parameters.  The  computer  program  is  tabulated  in 
Table  IV  —  HI  and  Table  IV - IV  shows  a  sample  output,  using  specimen  AI-4  data. 

I  .  Subprograms 

The  following  subprograms  ore  required  to  run  rhe  dynomic  analysis  computer  program: 

ETEMP  (f,  IFF) 

ALPHA  (T,  IFF] 

SN  (SDYN,  STEMP,  TEMP,  CTF,  IFF) 

CTEMP  (TCALP,  TC,  IFF) 

PROP  (OPT,  B,  H,  T,  A,  RJ ,  GAMAT,  XIP) 


Subprograms  ETEMP,  ALPHA,  and  SN  ore  presented  in  Appendix  V.  Subprograms  CTEMP 
and  PROP  were  developed  to  compute  skin  buckling  temperatures  and  stiffener  properties, 
respectively,  as  discussed  in  the  following  subsections. 

a.  Subprogram  CTEMP  -  This  program  computes  the  skin  critical  buckling  temperature 
csing  the  individual  alloy  curve  for  coefficient  of  thermal  expansion  versus  temperature. 
The  product  of  critical  buckling  temperature,  T  ,  and  coefficient  of  thermal  expansion, 
o  ,  is  calculated  from  Equaiion  (60),  or 

T  =  5 . 25  h  FH 

a2^2^*  +  ^ 


The  program  then  uses  a  mathematical  representation  for  a  ,  as  defined  In  Appendix  V,  to 
iterate  for  the  actual  value  of  a  to  be  used  in  computing  Tc.  The  volue  of  Tc  is  then 
returned  to  the  coding  program. 


TABLE  IV-I 

DYNAMIC  ANA-  \  SIS  COMPUTER  PROGRAM  INPUT  FORMAT 


CARD  I 


name 

NCASE 

IFF 

C  OL  (FORMAT) 

1(12) 

3(12) 

CARD  2 


NAME 

OPTX 

BX 

HX 

TX 

COL(FORMAT) 

1(12) 

3(F8.4) 

1 1(F8.4) 

1 9(  F  8 . 4) 

CARD  3 


NAME 

OPTY 

BY 

HY 

TY 

C  OL  (FORMAT) 

1(12) 

3(F8.4) 

1 1(F8 .4) 

19(F8.4) 

CARD  4 


NAME 

A) 

A2 

B2 

B1 

COl.  (FORMAT) 

KF3.4) 

9(F8 .4) 

1 7(F8.4) 

25(F8.4) 

CARD  5 


NAME 

TS 

RHO 

RNU 

DAMP 

COL  (FORMAT 

KF8.4) 

9(F8 .4) 

1 7(F8.4) 

25(F8.4) 

CARD  6 


NAME 

PSL 

T 

COL(FORMAT) 

1(F8.4) 

9{F8.4) 
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TABLE  tv-l  I 


NCASE 

IFF 


OPTX 

BX 

HX 

TX 

OPTY 

BY 

HY 

TY 

Al 

A2 

Bl 

B2 

TS 

RHO 

RNU 

DAMP 

PSl 

T 


DYNAMIC  ANALYSIS  COMPUTER  PROGRAM 
INPUT  PARAMETER  DEFINITION 

Two-digit  identification  number 

Alloy  identification  code 

=  1  Titanium  Alloy  (6AI-4V  Annealed) 

-2  Aluminum  Alloy  (7075-T6) 

Input  parameters  defining  stiffening  member  parallel  to 
x-direction  -  see  Subprogram  PROP  for  definition. 


Input  parameters  defining  stiffening  member  parallel  to 
y-direcfion  -  see  Subprogram  PROP  for  definition. 


Panel  bay  dimensions 

Skin  thickness  -  inch 

Weight  density  of  skin  and  stiffening  member  alloy  *■  Ib/in 
Poisson's  ratio  for  structure  alloy. 

Damping  ratio  for  structure. 

Spectrum  sound  pressure  level  -  bB, 

Structure  temperature  rise  -  °F  above  ambient. 
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TABLE  IV-III 


COMPUTER  PROGRAM  FOR  ELEVATED  TEMPERATURE 
DYNAMIC  RESPONSE  OF  STIFFENED  STRUCTURE 


1 

c 

THIS  PROGRAM  CALCULATED  THE  DYNAMIC  RESPONSE  OP 

c 

A  NIUE-BAY  FLAT  STIFFENED  PANEL  EXPOSED  TO  A 

< 

-1 

c 

UNIFORM  ACOUSTIC  PRESSURE  AND  A  UNIFORM  TEMP¬ 

4 

k. 

c 

r 

ERATURE  RISE.  ROOM  TEMPERATURE  IS  BO  DEGREES  F. 

O 

V 

c 

T  IS  A  TEMPERATURE  RISE*  ABOVE  ROOM  TEMPERATURE 

7 

c 

b 

c 

subprograms  required:  ALPHA <T* IFF ) *  ETEMPIT.IF' 

9 

c 

SNISDYU.STEMP*T.CTF*IFF ) ,  C  TEMP (  TC  ALP  *  TC  >  1 F1  ) 

10 

c 

AND  PROP  (ORT  *B*H.T*A*2J*  aC  *PI-'l 

n 

c 

12 

c 

FUNCTION  DEPUTATION 

13 

c 

1“ 

P  <  B  *  A  )  =B/A ♦A/U 

lb 

R(B<A*PR)=3.*(  i.  b.-PR*  *2  >  * (B/A+A/B ) *  '2-2.  ♦  (b.+HP  ) 

\  0 

1  * ( 1 . -PR  )  ) 

17 

20u 

HEAD  '  b  *  3r)  1  )  1 (CASE  *  IFF 

lb 

READlb*302)0PTX*BX*HX*TX 

19 

READ(b*302)OPTY*6Y.HY#TY 

20 

READ ( b  *  30 1  )Kk 

21 

READ(b.303)Al.A2*B2*Ul 

2. ' 

RE AD ( b  <  303 ) TS  *  RMO  *RNU  *  DAMP 

2  3 

READ(b.3D3)PSL*T 

24 

c 

INPUT  DATA  FORMAT  STATEMENTS 

2b 

301 

FORMAT (2121 

2b 

302 

FORMAT  < I2.3FB.4) 

27 

303 

FORMAT (4FB. 4 j 

2b 

C 

CALCULATE  SUBSTRUCTURE  PROPERTIES 

29 

CALL  PROP (ORTX*BX*HX*TX*AX*yj*wCX*Xl  ) 

30 

CALL  PHOP(OPTY*BY*HY.TYiAY .YJ*wCY*Yl  ) 

31 

H=TS 

32 

GMSKHO 

33 

PRSRNU 

34 

c 

CALCULATE  STIFFENER  STIFFNESS  AND  MAS. 

3b 

RXl-0.0b06U*Al*Al*XJ/(wCX*( 1 • +PR  )  » 

36 

RX2=0. 0b06G*A2*A2*XU/  (  WCX*  ( 1 .  t-PR  )  ) 

37 

RY1  =  0,05060*B1*B1*YU/(WCY*< 1 . +PR )  i 

3H 

RY2=0.ObO6s*B2*b2»YJ7  ( VvC  Y  ♦  I  1 .  +PR  )  ) 

39 

SKXSWCX* ( 1.+RX2  +  2.* (A2/A1 1  * ( 1  ,+RXl  )  )/A2 

90 

5KY=WCY* ( 1 .THY2f2.*<B2/B1 ) *< 1 ,+RYl )  )/B2 

41 

H3SH*H*H 

42 

STR2473.740B* ( 1 .-PR*PR ) * (SKX+SKY ) / (H3*A2*B2 ) 

43 

STR=STR/40n,0 

4', 

A3=(A1/A2)*'3 

4  b 

B3= (01/92 1**3 
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TABLE  IV -111  (CONT'D) 


4o 

GM=GM/386. 

47 

SKM=0 ,2b*GW*H*A2*A2+ ( 1 . +2 . *A3  +  2 . *B 3+4 . *A3*B3 ) 

4b 

1  +9.8b9o*GM*(Xl*A2*(l.+2.*A3)/(B2*B2) 

44 

2  +  YI*B2* < l.*2.*B3)/(A2*A2)  ) 

50 

C 

CALCULATE  COVER  SHEET  STIFFNESS,  and  MASS 

bl 

F22=F (B2.A2! 

b  C. 

F21=F IB2.A1 ) 

53 

F12=F  ( B 1  *  A2  ) 

54 

F11=F(B1.A1) 

k,  . 

V  * 

F  1S=A27*F2-J*2.*  (A1/A2)  *F21*F21 

bo 

1  +2  .*  IB1/B2)  +F12+F12+4.*  (A1/A2  )  *  (B1/B2  )  *F1 1  *F  1 

57 

C 

R0  TEMP  STIFFNESS 

5H 

Sku  =  2.02937*ETEMP(80.*IFF ) *H3* (F 1S+STR  ) / (  ;  i,-pp*pp  j 

o  t  * 

1  *A2»B2 ) 

' )  _ 

C 

ROC;!-’  TE'-‘P  FREQUENCY 

ol 

Fu  =  o.  lo4*SGRT  (SK'j/SKM)  . 

ad 

C 

CALCULATE  ROD".  TEMP  MEAN  SQUARE  STRESS  RESPONSE 

o3 

AR  =  3.  *  (B2/A2 )  +  ‘2  +  3.»-tA2/B2)*'*2  +  2» 

CD  4 

c 

CONVERT  DB  TO  RSI 

ob 

SRL=2.91*10.** (PSL/20.-9.  ) 

n 

c 

CALCULATE  ROOM  TEMP  DYNAMIC  STRESS  AT  X=G.Y=B2/? 

„  7 

SXU=-j.  3b»B2*B2*SQRT  (FO/GAMP  !  *SPL/  (H*H*AR  ) 

4  H 

c 

calculate  room  temp  dynamic  stress  at  x=a2/2.y=o 

SY0=1 « 30*A£*A2  +SQRT ( F 0/DAMp ) *SPL/ ( H*H*AR ) 

/u 

c 

CONVERT  STRESc  from  PSl  TO  KSI 

71 

SX0=SXU/10UU  .0 

72 

SYo=sro/iouu,o 

73 

c 

CALCULATE  room  temperature  life 

74 

call  SN(SXO*(J.O*8U.O»CTF1*IFF  ) 

7b 

CALL  SN(SY0»U.0»80.U#CTF2f IFF ) 

-'b 

Xl=A2/2. 

7  t 

Y  1  =  0  t  0 

7b 

X2  =  0.0 

74 

Y2=B2/2. 

MU 

STEMPzO, 0 

Ml 

c 

PRINT  RO1" M  TEMPERATURE  RESPONSE 

M2 

WRITE <:  .  -os 

M3 

GO  TO  a  .  i‘j 2)  «IFF 

b  4 

201 

WRITE  o  b)  NCASE 

fib 

GO  10  2.13 

bo 

202 

WRITE (fa »40o)  NCASE 

«7 

203 

WRITE (0(41 0 )  PSL ( T 

H" 

WRITE (6(41b) 

59 

WRITE (6(416)  FQ 

176 


TABLE  IV -III  (CONT'D) 


40 

41 

92 

93 

94 
9b 

96 

97 
46 

* 

lOu 

101 

102 

I  u  3 
104 
10b 
106 

10  7 
108 

109 

110 

I I  ; 

1  1  2 
113 

1 1  4 
1  lb 
1  1  6 
1  1  7 
1  18 
1  19 
120 
121 
122 

123 

124 
12b 
126 
127 
126 

129 

130 

131 

132 
13.. 
134 
1  3b 


WHITE.  <6  >420) 

WK ITE ( 6  *  42b ) 

WRITE (6*43u )  X2*  Y2*SX0  *GTEM,P#C1  •-  1 

riR  ITE  t  6  ►  4  30  )  X  1  •  Y  1  *  S  YO  «  STEMP  <  CTR2 

THERMAL  STRESS  EFfECTS 

R22=R (B2*A2*PR ) 

R21 =R (B2»A1 *PN ) 


R32=RtBl*A2*Pu) 

R 1 1=R ( B 1 * A  1 »PR )  „ 

F2S=F22* t  F22+2 • ♦ (A1/A2 ) *+2*F21+2«* (Bl/b2) *  *  2+F 12 

1  44  .  * ( A  1 /A2 ) * *2* ( B1/B2 I **2*F 1 1 ) 

RST”R2 2+2 . +A3*R 12+2. *B3*R2l  +  4  «*A3  +  B3+Rli 
R0= (F1S+STR ) /F2S 

calculate  critical  temperature  rise*  TCP 

TCALP-b.2b*H*!l*F22  /  1  A2*B2*  (  1 .  +PR  )  < 

CALL  CTEMP ( TCALP * TCR *  IFF ) 

#  ‘NOTE* '  TCA  AMD  RS  ARE  BUCKL IUG  TEMPERATURE 

AMO  TEMPERATURE  PATIO  F uR  A'.'  EQUAL 
SI2E  SIMPLE  PANEL.  R9  IS  TEMP  RATIO 

FOR  f  ill  lE— BAY  PANEL 

TCA=TCR/RO 
RS=T/TCA 
R9=T/TCR 
TACT=T+8Q  .  0 

CALCULATE  MATERIAL  PROPERTIES  AT  TEMPERATURE 
£S=ETEMP  (TACT'.  IFF  ! 

ALP=ALPHA ( TACT  *  IFF  ) 

D=0 . 0833+ES+H3/ ! 1 . -PR+PR  ) 

°CALCULATE  RESPONSE  FREQUENCY  AT  TEMPER ATURC-.T 
SKT=D*F2S*R0/ 1 A2*82 ) 

FQT  =  0 .809  +  SQRT  <SK  f /SKM) 

*  HJOTE*  *  FGT=FO*  ROOM  TEMP  FREQUENCY 
STLlN=“ES+ALP*T/(  1  .-PR)/10O'J.0 
IF(RS-R0)20b*20b*210 

PRE-BUCKLEO  RESPONSE 
>  FTEMP-FCT+(0.60  +  U. 4  0  *SOR T ( 1 • -R9 )  ) 

SXT=STLlN 

SYT-STLlN 

W 0=0.0 


GO  TO  21b 

c  post -Buckled  response 

210  FTEUP=FOT* ( 0.6C+0.44*SQRT (R9-1 , ) ) 
r  CALCULATE  PLATE  BUCKLING  AMPLITUDE*  WO 

*0=(3.37-0.20*RG) +H  +  SORT 1 F  2S*R  0* ( R9-1 . ! /RST 1 
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TABLE  IV-III  (CONT'D) 


1 3o 
137 
136 
139 


C  CALCULATE.  THHERMAL  STRESSES 

Cl  =  i ./ I  A 2*82* ! 1 . — PR*PR )  ; 

SXT=STLlN  +  0.81*ES*C  1*  <  (  2  . -PR*PR  )  *B2/A2+A2*PR/B2  ) 
1  *v.  0*rt0/10u0 . 0 


14  0 
141 
1  42 

143 

14:. 

1  44 

1  40 
14  7 
146 

144 
1  56 
lb] 
162 
lb3 

lb4 

1  bv, 

I  bo 
3b7 
lbs 

A  u  •  > 
i.  .i  *? 

160 

lol 

1=2 

lc3 

lo4 

lob 

16') 

167 

lob 

lo9 

170 

17] 

172 

173 

174 

375 

17o 

I I  7 
1  7f; 
174 
160 


C 

21b 


C 

C 


C 

4  0  Li 


405 

406 
410 

41b 
4  16 
420 


425 

430 

435 

440 


445 


SYt=STL1N+ 1 .36*ES*C1* (PR *82/ A 2+  t 2 , -PR+bh ) +A2/B2 ) 
1  *W0*W0/1£)00.G 

CALCULATE  DYNAMIC  STRESS 
CONTINUE 

C2=SQRT ( FTEMP/F 0  3 

SX0=C2*SXO 

SY0=C2*SYu 


calculate  elevated  TEMPERATURE  life 

CALL  SN (SX0.SXT .TACT. CTF1 , IFF ) 


CALL  SN(S<0.SYT. TACT, CTF2, IFF) 

PRINT  ELEVATED  TEMPERATURE  RESPONSE 
wR  I  TE  (  6  *  4 3b  ) 

TA=TCR+80 . 0 
WRITE (o» 440 3  TCR 

W  R  I  T  E  (  6  >  4  4  5  )  .'.o 

.-.RITE  (6.416)  FTEMP 
/.RITE  (o.  420) 


.•.RITE  (6. 425) 

WRITE (6,430 )  X2.Y2.SX0.SXT.CTF1 
.-.'RITE  (6. 4  30  !  XI .  Y1  ,SY0  ,SYT  ,CTF2 

60  TO  200 


FORMAT  STATEMENTS  FOR  OUTPUT  DATA 
FORMAT <» 1  * ,25X. ’DYNAMIC  RESPONSE  OF  A«,/,19X» 

1 ’NINE -bay  stiffened  panel  exposed  TO  * . / »  2 IX  f 

2  *  ACOUSTIC  EXCITATION  AND  HEATING',/) 

FORMAT (29X, 'DATA  CASE ' . I  4 , // , 27X , 'MATER  I AL  l  TITANIUM') 
FORMAT (29X, 'DATA  C ASE ' • I  4 , // , 27x , 'MATER  I AL  t  ALUMINUM') 
FORMAT ( 5X » 'EXCITATION  SPECTRUM  LEVEL  =  »  ,  F4 . 0  ,  IX  .  •  DB ' . 
13X, 'TEMPERATURE  INCREASE  =  '  r F4 . 0 , IX »  'DEG  .  F',//) 

FORMAT (24X, 'ROOM  TEMPERATURE  RESPONSE',/) 

FORMAT ( 20X FUNDAMENTAL  FREQUENCY  =',F7.1,'  hZ'»/) 
FORMAT ( 5X , 'STRESS  AT  POINT ' f 3X , 'DYNAMIC  STRESS'. 3X, 

1 'THERMAL  STRESS  '  » 3X , 'CYCLES  TO  FAILURE') 

FORMAT  1 BX  » 'X ' ,7X ,  ' Y»  » 1 IX , 'KSI '  » 14X  » 'KSI ' »/) 
FORMAT(6X»F5.2»3X»F5.2»5X»F8.3,9x,F8.3»10X»1PE9,2»/) 
FORMAT  <////»  2  2  X . 'ELEVATED  TEMPERATURE  RESPONSE',/) 
FORMAT ( 10Xf 'BUCKLING  TEMPERATURE  =',F8.2,'  DEG.  F' 

1»*  ABOVE  ROOM  TEMPERATURE',/) 

FORMAT ( 18V , 'BUCKLING  AMPLITUDE  =«,F8.4, 

1'  INCHES1’./) 

END 
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TABLE  IV-IV 

OUTPUT  FORMAT  FOR  DYNAMIC  ANALYSIS  COMPUTER  PROGRAM 


DYNAMIC  RESPOND!  OF  A 
NINE-BAY  STIF  EN;  D  PANEL  EXPOSED  TO 
ACOUSTIC  E  XL  I  T  A i  I  ON  AND  HEATING 

DAT  C/.SL  4 

MATER  I  '.L  :  ALUMINUM 

EXCITATION  SPECTRUM  LEV  L  =  135.  DR  TEMP:  RATURE  lilCREAS-  =  2  0.  DEG.  F 


P0  M  TEMPERATURE  i  SPONSE 
FUNDAMENT.,!.  FREQU  I  ICY  =  17S.  .  M2 


STRES 

X 

AT  POINT 

Y 

DYNAMIC  STRES 
KSI 

THERMAL  STRES 
KSI 

CY  LES  TO  FAILURE 

.0 

6 .  L 

6.  0  7  S 

.0 

6  •  4  3  +  U  t  > 

3 .  U 

.1, 

5.485 

.0 

1.03+  7 

EL'  VATED  TEMPLRATURE  RESPONSE 

BUCKLING  TEMPERATURE  =  1  i  .  0 -  DEG.  F  ABOVE  RO.  M  TEMPERATURE 

BUCKLING  AMPLITUDE  =  .25:  INCHES 

FUNDAMENTAL  FREQUENCY  =  417.5  HZ 

STRESS  AT  POINT  DYNAMIC  STRESS  THERMAL  STRES  CY< LES  TO  FAILURE 

X  Y  KSI  KSI 

.01,  6.0-  9.371  -16.83",  1.32  +  05 

.0  L 


3.01: 


8.46U 


-24.342 


1.40+  .5 


The  input  parameters  are: 


TCALP  -  Product  of  critical  buckling  temperature  and  coefficient 
of  thermal  expansion. 

IFF  -  Al loy  Code 

=  I  for  Titanium  Alloy  (6AI-4V  Annealed) 

-  2  for  Aluminum  Alloy  (7075-T6) 

The  output  parameter  is: 

TC  -  Critical  buckling  temperature  -  °F  above  Ambient 
The  subprogram  is  listed  in  Table  IV-V. 

b.  Subprogram  PROF  -  This  program  computes  stiffening  member  properties  such  os  area 
and  moment  of  inertia.  The  basic  relations  are  from  Reference  8,  Appendix  I.  two 
different  sectional  shapes  are  available,  a  zee  or  a  channel  section,  with  the  parameters 
described  in  Figures  A— 1  —  1  and  A-l-2,  respectively,  of  the  referenced  report. 

The  input  parameters  are: 

OPT  -  Option  code  to  select  sectional  shape 

-  0  for  zee-section 

-  I  for  channel  section 

B  -  Flange  width  of  stiffening  member  -  in 

H  -  Height  of  stiffening  member  -  in 
T  -  Thick  ness  of  stiffening  member  -  in 

The  output  parameters  are: 

2 

A  -Cross-sectional  area  -  in 

RJ  -  St.  Venant's  Torsion  Constant  -  in  ^ 

GAMAT  -  Warping  constant  for  thin  walled  open  section  beam, 
with  the  pole  taken  at  the  shear  center  -  in  6 

XIP  -  Polar  moment  of  inertia,  referenced  to  rotation  about  the 
attachment  point  -in  ^ 

The  listing  of  this  program  is  presented  in  Table  IV-VI. 
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TABLE  tV-V 

COMPUTER  PROGRAM  FOR  CALCULATING  SKIN  BUCKLING  TEMPERATURE 


1 

SUBROUTINE  C TEMPI  TC  ALM>  TC  *  IK  5 

2 

C 

THIS  SUBROUTINE  CALCULATES  S*l!i  HUCklInG 

.5 

r 

TEMPFRATUPE  FOP  ALU'-*, I MUM  OR  TITANIUM  ALLOT 

u 

C 

STRUCTURAL  PAUFLr- . 

6 

c 

r> 

c 

TCALP  -  PRODUCT  OF  RuCkLJuG  TEMPERATURE 

7 

c 

AMD  ALPHA  FROv  CAL'  I  MG  PROGRAM 

H 

c 

TC  -  Pi  :CKL  IMG  TEMPER  AT' mE  -  0;"A  F  AMOVE 

Q 

c 

H  O  TE'-V-'ERAT'  iPE 

in 

c 

IF*-  -  /'Ll  GY  CODE 

1 1 

c 

=  1  titanium 

12 

c 

-  2  aluminum 

13 

c 

1*+ 

TC-0.0 

is 

GO  TO  ( 10  '  .20  U  .  IF: 

16 

c.  **♦ 

♦  tit. 

1? 

V- 

MATFRIAL  n A L “ 4 v  TITAtiliJM  AM!. BALED 

lfl 

r 

TEMPI.  RATurE  LIMITATION  100..  DEG.  F 

1  Q 

c 

2  0 

10-:- 

C  1  =  4 ,  4  5  F  —  0  o 

21 

C2=4 .a3F-0Q 

2-' 

1=1 

23 

1 

C  3  -  6  0  ,  *  C  2  +  C 1 

24 

TC= . 6 * SOR T  <  ! C3/C2 ) * • 2*4 . +TCALP/C2 

25 

T  =  TC  +  80  . 

26 

I F ( T - 2  6  0 .  )  50.50.2 

27 

2 

IF ( T -360  •  )  3.3.5 

28 

3 

0=4 . 9E-06 

?Q 

C2=2  < 5E-Oa 

30 

1  =  1  +  1 

31 

IF(I-,2)  1,1,50 

32 

5 

ClaS. HOE-06 

33 

TC= TCALP/ Cl 

34 

50 

TC  =  TC 

36 

f?  F  T  !  )R!  i 
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TABLE  IV-V  (CONT) 


36 

37 

38 

39 

40 

41 

42 
4  3 

4  4 
4*5 
Hh 
-*7 
u« 
4<J 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
6 1 


r  ****w*****> 

C  MATERIAL  7075-T6  ALUMINUM  ALLOY 

c  TEMPERATURE  LIMITATION  600  DEG.  F 


200 


201 


202 

203 


204 

205 


206 


81=12. 4E-06 
F2S5.0E-09 
1 1  =  1 

TC  =  0 • 50*SGRT  <  ( P3/F2 ) * *2  +  4 . 0 *TC ALP/F2 ) *0 . 50+F3/F2 

T=TO80 . 

IF ( T-100 « ) 500 •  50 0  » 202 
IF ( T-300 . )203  »203  »  204 
F1=19.625E~06 
F2=0.75E-O9 
11=11+1 

IF ( 1 1-2 ) 20 1 «  20 1 • 500 
IF!T-400. 1205 #205 *206 
Fl=13.0E-06 
F2=l . 5E-09 
11*11+1 

I F (  1 1-2 )201 » 201 *500 

Fl=13.6E-06 

tc=tcalp/fi 


500  TC=TC 
RETURN 
END 
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TABLE  IV-VI 

COMPUTER  PROGRAM  TUP  CALCULATING  SECTION  PROPERTIES 

SllHPO'-iT  IMF.  PROP!  OPT  .  4.  H.  T.  A.  P  J.  GAMAT  .  X  I  P) 

C 

C  SECTION  PPfjPfPT  IES 

C  I c  nnT  =  0  7  E  t  SECTION 

C  !  F  opt  t  l  Channel  cf^T  T  on 

C  H-  5TPIN0EB  HEIGHT.  CL  TO  CL 

r  »:  f|  AnGE  aIOTH 

C  t=  stpImgfp  thic*mESo 

C  REFERENCE !  AFEDL-TP-71-1  07 

C  wC  -  aARPInG  CONSTANT  ABOUT  SH^AR  CENTrp 

c  (jAvAT  -  .vABPI'iG  C^TTAiiT  ABOUT  ATTACK  P'"T‘iT 

r  A  —  C»ac.>  t  p  C  T  I  *  f ; A i  APE  A 

C  Rj-  ST.  VE  f  |A  I  ,Tr  TOP  a  I  O'  i  (■'••rTAi  .T  -  J 

c  y  ]  P  -  POLAR  Mr,urr-T  OE  I»iE«T  IA  -V  R  a.  i  T  ATTACH  POT :  IT 

T  F  (  OPT '  1  i|  ,5 

C  7ER  AT  I  =  f  F  •  ;f  ° 

1  . -T* ! H  +  ? , *B  > 

r.  ;h+*2*  (a  .  *h+hi 
f  l=T**2«'3.*h-|?.*H) 
c  y Is ( T» ( 0*01 '  ' /12 . 

02=2 ,*H+T 
03=2. *B-T 

y?  Ti-;  t*h«-,2«D3  !  75  . 

7  7  I  =  !  T  / 1  2  •  )  *  <  G  •  *  a  «  *  3  ♦  rt*  t  *  *  2  ) 

P j= ( T**3/3 . : * ( 2 . *R+H) 

wCiT»R»*3*H**2* (B*p .  *Ht / ( 12 .* i?.*R*hii 
Sx=B'2. 

E7=-a/2 . 

D4SGX**2+S2*  »2 
D5=D4*A 

x ips  tx 'i+7^ i+oa > 

GAVATSV.C+  (  57**2  )  *7  7  1-2  . *Gx*c7*x7  i+  (  EX*  »2  >  »x»  I 

DC  J  IP". 

C  CHANNEL  SECT  I  oil 

2  c-2  .  *a+H 
XHAO=R**2/c 

Fl=6  .  *B*H 

E  =  3.*R**2  7:rl 
CX=E+XBAP 
SxsE+ (fl/2 . 1 
S2=-h/2 . 

Ex=cx-sx 
AsT* ( h+2 . *B! 

F2=3.*h+2.*B 

XXl=T* (H**2*F1+T**2*F2)/12. 

E3=12.*H*XBAO**2+8.*R**3 

F4sB-Xc'AP 

E5sH-2 . *XBAP 

77  I  =  T*  t  E3-24  ,*XBAR*B*F4  +  1?  .  *R*FB*T+B  .+F4*T*x2+T*+3l/12. 
Pj=T**3*E/3. 

Fa--3.*B  +  2.*H 

wC'-  T*B*  *3*H*  «2  *EB/ (  12.*E1  1 

GAMAT  =  aC+  (  S7**2  >  *77  1+  (  c,x*  *2  )  I 

F7s ( EX* *2  + S’* *2 1  *  A 

X IP= ( X X T+77 I+F7 ) 

pETlipn 
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APPENDIX  V 


TEMPERATUR!  DEPENDENCE  OF  MATERIAL  PROPERTIES 

The  analyticol  development  presented  in  Section  II  focused  attention  upon  the  tempera¬ 
ture  effects  associated  wish  the  state  of  strain  of  the  structure.  The  results  obtained  from 
these  considerations  indicate  that  the  structural  stiffness  decreases  with  temperature  rise  up 
to  the  critical  temperature,  and  that  tor  a  temperature  rise  above  the  critical  temperature,  the 
structural  stiffness  increases. 

Beyond  these  considerations,  certain  alloy  properties  ore  temperature-dependent, 
and  this  dependency  must  be  included  to  predict  structural  life  accurately.  Algorithms 
were  developed  to  represent  the  temperature  dependence  of  each  property  considered, 
and  comouter  programs  were  written  to  expedite  their  use  with  the  dynamic  analysis 
computer  program  detailed  in  Appendix  IV. 

The  material  properties  considered  are  the  elastic  (Young's )  modulus  ,  E,  the  coefficient 
or  thermal  expansion,  .1  ,  and  trie  characteristics  of  the  fatigue  curves  for  each  alloy. 

Sic'rc  properties  were  obtained  from  Mi  L-HDBK-5B^,  and  generally  agree  with  those  of 
AFML-TR-6S-1 1:3^.  The  analytical  expression  for  fatigue  life  is  based  on  the  mean 
stress  variation  inferred  from  MIL-HDBK-5B  and  the  fctigue  characteristics,  with 
temperature  variation,  resuiting  from  the  coupon  fatigue  testing  discussed  in  Section  III. 3. 

A.  Static  Material  Properties 

This  discussion  will  concern  the  effect  of  temperature  on  the  static  material  properties 
of  7075-T6  aluminum  and  6A1-4V  annealed  titanium  sheet  material.  The  assumed  ambient 
temperature  for  this  development  is  80°F. 

1 .  Elastic  Modulus 

The  effects  of  temperature  on  the  elastic  modulus  of  7075-T6  aluminum  and  6A1-4V  annealed 
ti tan i urn  alloys  ore  shown  in  F igure  V-l .  These  curves  are  reproduced  from  Ml L-HDBK-5B, 
Figures  3. 2. 7. 1.4  and  5. 4.6. 2.  4,  respectively  . 

Tiie  elastic  modulus  curve  for  7075-T6  aluminum  decreases  steadi ly  with  increasing  tempera¬ 
ture.  Since  a  single  straight  line  representation  would  not  suffice,  the  temperature  range 
was  broken  into  three  regions.  For  each  of  these  temperature  regions,  the  elastic  modulus 
■was  represented  by  a  straight  line  element  as  follow: 

Room  temperature  to  20CPF  (80  1  T  200) 

F  =  1.020  -  0.0003  T  J  Ea 

bi  200°  to  400° F  (200  i  t  ■_  400) 

E  =  0.960  -  0.0007  (T  -  200)  E 

o 

c)  400°  to  600° F  (400  T  _  600) 

E  -  I  0.820  -  0.0016  (T  -  400)  I  EQ 
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PERCENT  OF  ROOM  TEMPERATURE  MODULUS  (E  ) 


where  (0  10,3  x  10^  p»i  ii  the  room  temperatirre  elastic  modulus, 

T hr?  clastic  modulus  for  titanium  hat  a  linear  relationship  with  temperature  below  800'F. 
Hence,  the  titanium  modulus  it  represented  by  the  following  tingle  relation: 

a)  Room  temperoluie  to  BOC^F  (80  T  800) 

£  Ml, 030  -  0,000375  T)  £*, 

where  the  room  temperature  modulus  it  L0  »  16.6  x  10^  pti. 

The  above  representations  of  the  elastic  modulus  were  developed  into  a  digital  computer 
program  for  vtr  by  the  analytical  program  discussed  In  Appendix  IV,  This  function, 
entitled  "LUMP",  it  presented  in  Table  V-l,  vine#  this  program  iso  function,  the  input  for¬ 
mal  is  timply 

MiMPu,  iff ) 

in  the  colling  program,  where 

T  -  Input  lempnratu'e  oi  which  the  elastic  modulus  it  desired  *  °f 
If  f  -  Alloy  code 

«■  I  Titaniom  Alloy  (6AI-4V  annealed) 

2  Aluminum  Alloy  (7075-T6) 

Data  points  calculated  from  this  program,  In  temperature  Increments  of  SO^T,  ore 
compared  with  the  MIL-HDBK-38  curves  In  Flgufe  V*l, 

jh_  Coelfic lent  of  Thermal  Expansion 

7 Is**  effects  ol  temperature  on  the  thermal  expansion  of  aluminum  and  titanium  alloys 
are  shown  in  Figure  V-2.  These  turves  ore  reproduced  from  MIL' HDBK-56,  Figures 
3,  2.7,0  and  5.4.6,  ! ,  respet lively. 

The  coefficient  o*  no* m<t I  expansion,  >»  ,  for  aluminum  alloy  Increases  exponentially 
with  increasing  temperature  as  shown  in  Figure  V-2,  Hence,  the  temperature  range 
wus  divided  Into  three  regions  and  the  following  representations  were  derived: 

'it  Room  temperature  to  lOC^F  (80  -  T  -  100) 

-  II?, 4  <  0.0050  1  |  I0'6 

L)  I  OOP  to  WPI  (100  -  1  -  300) 

|l?,V  r  0. 002/5(1  -100)|  10'6 

c)  300°  to  400Mt  (300  -  T  -  400) 

••  I  13,45  <  0.0015(1-300)  I  10*6 
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1 


TABLE  V-l 

COMPUTER  PROGRAM  FOR  CALCULATING  ELASTIC  MODULUS 


t 

FUNCTION  ETEMP  (  T »  IF!') 

2 

C 

3 

c 

THIS  FUNCTION  COMMUTES  ELASTIC  MODULUS  FOR 

.3 

4 

V_ 

aluminum  or  titanium  allot  as  a  function  of 

5 

C 

TEMPERATURE 

-rj 

ft 

c 

7 

c 

T  -  INPUT  TEMPERATURE  -  DEG.  F 

i 

8 

c 

IFF  -  ALLOY  CODE 

0 

c 

=  1  titanium 

10 

c 

r  2  ALUMINUM 

| 

11 

c 

I, 

12 

GO  TO  (100  *200) *  IFF 

13 

C  »♦*****♦  ******•************♦•*♦♦*  ***+****  +  •*♦♦**  *  M  +  1 

J 

14 

c 

MATERIAL  GAL-4V  TITANIUM  ANNEALED  SHEET 

IS 

c 

REFERENCE  MIL-HDBK-5B 

1ft 

c 

TEMPERATURE  LIMITATION  SOD  DEGREES  F 

| 

17 

c 

PT<T <ft00  F 

1ft 

10n 

IF (Q00-T ) lftO *  190  *  150 

19 

180 

ETEMpr ( 1,030-0. 0O0375*T)*lft.ftE+0ft 

-jn 

RETURN 

1 

21 

C 

T>800  F 

T 

27 

1B0 

ETEMP=12. lE+Oft 

-j 

23 

WR  ITE ( ft  *  333 ) 

24 

RETURN 

i 

25 

C  **  ***♦**'♦  ****************  *****************  *  ******  * 

- 

2ft 

c 

MATER  I Al  7075-Tft  SHEET 

2  7 

c 

REFERENCE  mIl-HDBK-513 

2ft 

c 

TEMPERATURE  LIMITATION  ftOO  DEGREES  F 

29 

c 

PT<T<20O  F 

30 

20t; 

IP  (20O-T)  2.30, 21 0*210 

. 

31 

210 

ETEMPr ( 1.020-0.00H 30 *T) *10. 3E+ 0ft 

32 

RETURN 

33 

c 

200<T<400  F 

34 

220 

IF(40'J-T)240*230*230 

35 

230 

ETEMP:  (0 ,9ft-0 .00 '-70*  (  T-20n  )  )  *10 . 3E+06 

3ft 

RETURN 

37 

C 

4  0 11  <  T <ftOO  F 

3ft 

24  0 

IP (60n-T )2ft0 ,250 *250 

39 

250 

ETEMP: (0.82-0. On  1ft* ( T-40O )  )  *10.3E  +  0ft 

40 

RETURN 

41 

T  >#5  0 1  ■  F 

42 

2ft0 

ETEmp=0. 50*10. 3E+0h 

43 

WR  I T  E,(  ft  *  353  ) 

4  4 

3  S3 

FORMAT!/* 5X«  »UP”EP  TEMP  LIMIT  ON  ELAST  MODULUS'* 

45 

1  '  E  X  C  E 1  OED' */ ) 

4  5 

RETURN 

47 

END 
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-  COEFFICIENT  OF  THERMAL  EXPANSION  -  10‘6  IN/IN  /°F 


FIGURE 


a)  ALUMINUM  ALLOY  7075-T6 


b)  TITANIUM  ALLOY  6AI-4V 

-2.  TEMPERATURE  EFFECTS  ON  COEFFICIENT  OF  THERMAL  EXPANSION 

I  FROM  i-'E  FI: PENCE  I  7) 
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The  coefficient  of  thermol  expansion  for  titanium  increases  from  200°  to  40CPF  and  is 
a  constant  thereafter.  As  no  data  were  available  for  temperatures  less  than  20Q°F, 
the  slope  of  the  segment  ar  200°F  was  projected  to  room  temperature  to  develop  a 
relationship  for  the  computer  program.  The  resulting  analytical  representations  are: 

a)  Room  temperature  to  200°F  (80  T  200) 

»  =  (  4.45  <  0.00425  T  )  10'6 

b)  200°  »o  400° F  (200  T  -  400) 

«  -  (4.80  +  0.0025T)!0"6 

c)  400°  to  1 000°F  (400  T  -  1000) 

“  =  5.80  x  10*6 

The  preceding  representations  of  the  coefficient  of  thermal  expansion  were  formulated 
into  a  digital  computer  program  similar  to  that  for  the  elastic  modulus.  This  "ALPHA" 
function  is  presented  in  Table  V-ll.  The  input  format  is 

AlPHA  (T,  IFF) 

in  the  calling  program,  where  T  and  IFF  are  as  defined  in  the  preceding  subsection. 

B.  Fatigue  Characteristics 

The  analytical  description  of  the  fatigue  curves  described  here  resulted  from  a  review  of 
fatigue  data  presented  in  MIL-HDBK-5B  and  the  data  from  the  coupon  fatigue  tests 
ot  roam  and  elevated  temperatures,  ft  was  assumed  that  for  all  values  of  mean  stress,  cm, 
dynamic  stress,  a,  and  temperature,  T,  the  resulting  fatigue  curve  would  be  linear 
when  plotted  on  log-log  scales. 

1 .  Mean  Stress  Effects 

The  effect  of  increasing  mean  stress  is  basically  a  decrease  in  fatigue  life  at  constant 
dynamic  stress.  Figures  V-3  and  V-4  present  fatigue  curves  for  axially  loaded  aluminum  at 
room  temperature,  and  titanium  at  room  temperature  and  600°F.  These  curves  were  obtained 
from  MIL-HDBK-5B  in  the  form  of  constant  amplitude  test  data  and  were  converted  to  an 
equivalent  random  amplitude  fatigue  curve  by  the  method  of  Reference  1  .  Based  on  these 
room  ond  elevated  temperature  curves,  it  was  determined  thot  the  effect  of  increasing 
meon  stress  was  to  lower  the  RMS  stress  level  by  a  constant  amount  from  the  zero  mean 
stress  curve.  This  decrease  in  dynamic  stress  was  approximately  1  ksi  for  every  10  ksi 
increase  in  meon  stress.  Then,  if  only  a  zero  mean  stress  fotigue  curve  were  available, 
as  is  the  case  for  the  coupon  fatigue  data  of  Section  III,  the  dynamic  stress  can  be  corrected 
by  subtracting  0.1  <7m  from  the  dynamic  stress,  where  Om  is  the  mean  stress.  The  equation 
for  the  zero  mean  stress  fatigue  curve  is  of  the  form 

log  N  =  A  +  B  log  or  tv-) ) 

Then  the  equation,  including  mean  stress  effects,  becomes 

log  N  =  A  +  B  log  (  o  -0.1  an()  (V-2) 
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TABLE  V-ll 

COMPUTER  PROGRAM  POR  CALCULATING  COEFFICIENT  OF  THLRMAL  EXPANSION 


1 

FUNCTION  ALPHA ( T  # IFF) 

2 

C 

3 

C 

THIS  FUNCTION  COMPUTES  COEFFICIENT  OF  THERMA 

4 

c 

EXPANSION  FOR  ALUMINUM  OR  TITANIUM  ALLOY  AS  , 

5 

c 

FUNCTION  OF  TEMPERATURE 

6 

c 

7 

c 

T  *  INPUT  TEMPERATURE  -  DES.  F 

8 

c 

iff  -  alloy  cooe 

9 

c 

=1  titanium 

10 

c 

=2  ALUMINUM 

11 

c 

12 

GO  TO  (100*200) . IFF 

13 

c  ***************** ***  ******************  ******* t* 

14 

c 

MATERIAL  6AL-4V  TITANIUM  SHEET 

15 

c 

annealed 

18 

c 

REFtRENCE  MIL-HDBK-5B 

17 

c 

TEMPERATURE  LIMITATION  1000  DEGREES  F 

18 

c 

RT<T<200  F 

19 

100 

IF(200-T>180. 350*150 

20 

*50 

ALPHAS  <4, 45+0. 00425 *T) *1.0 E-06 

21 

RETURN 

22 

c 

200<T<400  F 

23 

180 

IF(400-T> 185*190# 190 

24 

190 

ALPHAS) 4, 80+0. 0025*T>* 1.0 E-06 

25 

RETURN 

26 

C 

400<T<1000  F 

27 

185 

IP(1000-T)195*198#198 

28 

*98 

WRITE (6*500 ) 

29 

500 

FORMAT l A  *  5X  # ' UPPER  TEMP  LIMIT  ON  COEFF  OF  EXPAN 

30 

1 ’EXCEEDED*  */> 

31 

196 

ALPHAS5.8E-06 

32 

RETURN 

33 

c  *********************************** ************* 

34 

c 

MATERIAL  7075-T6  SHEET 

35 

c 

REFERENCE  MIL-HDBK-5B 

36 

r 

TEMPERATURI  LIMITATION  600  DEGREES  F 

37 

C 

RT<T<100  F 

38 

200 

I F  (  100-T)2£*0*210*210 

39 

210 

ALPHAS (12 ,4+0. 0050*T> *1.0 E-06 

40 

RETURN 

41 

C 

100<T<300  F 

42 

220 

IF(300-T)240*230  #230 

43 

230 

ALPHAS (12.9+0 . 00275* ! T-l 00 ) ! *1 .0E-06 

44 

RETURN 

45 

C 

300<T<4QO  F 

46 

240 

IF(400-T)260#250*250 

«*? 

250 

ALPHAS (13. 45+0, 00150* (T-300) >*t.0E-06 

*• 

RETURN 

** 

C 

T>400  F 

■0 

260 

ALPHAS13.6E-06 

51 

IF(6P0-T)280#270*27O 

52 

260 

WRITE ( 6 i 500 ) 

153 

270 

RETURN 

■4 

€HD 
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AXIAL  LOADING 
DATA  COVERTED 
TO  EQUIVALENT 
RANDOM  AMPLITUDE 


RANDOM  AMPLITUDE  FATIGUE  CURVE 
6AI-4V  TITANIUM 


This  relation  was  found  to  be  valid  for  both  the  aluminum  and  titanium  data  of  Figures 
V-3  and  V-4, 


From  the  riveted  coupon  fatigue  tests  of  Section  III,  the  equotions  for  the  regression  lines 
of  the  room  and  300°F  temperature  fatigue  curves  were  as  follows  for  the  aluminum  data: 

log  N  =  10.42  -  4.61  log  O  (80oF) 

(V-3) 

log  N  =  9.13  -  3.85  log  a  (300°F) 

Since  only  two  temperature  points  were  available,  it  wos  assumed  that,  if  the  regression 
line  was  linear  on  a  log-log  scale,  then  the  slope  and  the  intercept  of  each  equation 
varied  linearly  with  temperature.  Therefore,  the  intercept,  A,  and  slope,  B,  con  be 
expressed  as  functions  of  temperature  by 


C1.+C,2T 

Using  the  constants  in  Equations  (V-3)  gives 
Cn  +  80  C  =  10.42 
CM  +  300  Z  -  9.13 
Solving  Equations  (V-4)  gives  the  constants 
Cn  -  10.89 
C)2  = -0.00584 

and  the  intercept  and  slope  become 
A  =  10.29  -  0.00584  T 


and  +  C2  T  ~~  ® 


C21  +  80  C22  =  -4.61 
C2)  +  300  C22  r  -3.85 


C2]  ”  -4.89 
C22  =  0.00347 


B  =  -4.89  +  0.00347  T 


Then,  the  general  expression  for  the  fatigue  curve,  at  zero  mean  stress,  is 

log  N  =  (10.89  -  0.00584T)  -  (4.  89  -  0. 00347T)  log  a 

Introducing  the  mean  stress  variation  as  assumed  for  Equation  (V-2),  the  fatigue  curve  is 
given  by 

log  N  =  00.89  -  0.00584  T)  -  (4.89  -  0.00347  1;  log  v  a  -  0. 1  o  )  (V-5) 

where  N  =  life  in  cycles  to  failure 
a  -  dynamic  stress,  ksi  rms 
a  =  mean  stress,  ksi 
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T  -  temperature  in  degrees  Fahrenheit. 

The  equations  of  tiie  regression  lines  for  the  titanium  alloy  riveted  coupon  fr.tigue  data 
o  f  Section  1 1 1  are 

log  N  -  12.28  -  5.26  log  a  (80°F) 

(V-6) 

log  N  -  10.33  -  4.34  log  a  (600°F) 

Repeating  the  procedure  given  for  the  aluminum  alloy  gives  the  following  constants: 


CM  -  12.58 


C| 2  =  -0.00376 


C2r-5.40 

C22  =  0.00176 


and  then 


A  -  12.58  -  0.00376  T 
B  =  -5.40  *  0.00176  T 

The  fatigue  curve  for  titanium  alloy,  at  any  mean  stress  or  temperature,  is  then  given  by 
log  N  =  (12.58  -  0.00376  T)  -  (5.40  -  0.00176  T)  log  (  o  -  0. 1  a  )  (V-7) 

where  all  parameters  are  as  defined  for  the  aluminum  alloy. 

These  relations  were  formulated  into  a  computer  program  for  use  with  the  general  analysis 
program  of  Appendix  IV.  This  program,  in  the  form  of  a  Subroutine,  is  presented  in 
Table  V-lll  and  is  entitled  "SN."  The  input  parameters  are 


SDYN 

STEMP 


Dynamic  Stress  -  ksi  rms 
Thermal  (or  Mean)  Stress  -  ksi 
Temperature  -  °F 
A I  loy  Code 

1  1  i tanium  (6A1 -4V  annealed) 

2  Aluminum  (7075-T6) 


The  output  to  the  calling  program  is 


Life  in  cycles  tofailuie 
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TABLE  V-lll 

COMPUTER  PROGRAM  FOR  CALCULATING  FATIGUE  LIFE 


1 

2 

3 

4 

s 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 
18 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 


subroutine  sn<sdyn*stemp»temp#ctf.ipf) 


c 

c 

c 

c 

c 

c 

c 

c 


THIS  SUBROUTINE  CALCULATES  FATIGUE  LIFE  FOR 
ai  i  iu  t  mi  im  OR  TITANIUM  ALLOY  AS  A  FUNCTION  OF 
?EmpeSa?ur?  and  mean  STRESS.  THIS  SUBROUTINE 
IS  BASED  ON  COUPON  FATIGUE  TEST  DATA  AT  ROOM 
awn  flEVATED  TEMPERATURE. 


ROOM  TEMPERATURE  IS  BO  DEG.  F 


C 

c 

c 

c 

C 

c 

c 

c 

c 

c 

c 

c 


SOYN  -  DYNAMIC  STRES  -  KSI  (RMS> 

STEMP  -  THERMAL  (OR  MEAN)  STRESS  *  *SI 
TEMP  -  TEMPERATURE  -  DEG,  F 
CTF  -  LIFE  IN  CYCLES  TO  FAILURE 
IFF  -  ALLOY  CODE 
=  1  TITANIUM 

=  g  aluminum 


GO  TO ( 100  *  200 ) *  IFF 

************************** *********** 
MATERIAL  6AL_4V  TITANIUM  SHEET 
TEMPERATURE  LIMITATION  600  DEG. 


*************** 

ANNEALED 

F 


100  C 1= 12 • 58-0 . 00376* TEMP 
C2=-5.40+0.00*76*TEMP 
ARF=Cl  +  e2*ALOG10  f SDYN-0 . 1 *STEMP > 

ctf=io.**arf 


C 

C 

C 

C 


*  **!*!*  **** ***** ******* ************************* ****** 
MATERIAL  7075-T6  ALUMINUM  SHEET 
TEMPERATURE  LIMITATION  300  DEG.  F 


200  C1=10.89-0,00584*TEMP 

C2=-4.89+C.O 0347* TEMP 
ARF=C1+C2*ALOG10 (SDYN-0. 1*STEMP) 


CTF=10 . **arf 


RETURN 

END 
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